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PREFACE 


The reports in this volume represent two years of work under a 
contract with the Air Force Cambridge Research Center. These reports 
are presented as separate investigations, and will be submitted 
as such to scientific journals. 

They form a coherent group of studies, however, much more 
so than would appear from reading them separately. This is due 
partly to the overlap of subject matter but more importantly to the 
fact that they illustrate or depend upon the same general principle*. 

Our main interest is to bring out matters of general importance 
although these may appear to the reader to have become lost in a 
mass of details. Consequently the new results will be summarized 
in this preface. 

The research has been (1) a systematic study of the optical spectra 
of transition metal ions in corundum (2) a systematic study of 
transition metal ionB in wurtzite type lattices (3) an interpretation 
of the absorption spectrum of nickel fluosilicate (4) a study of the 
exchange coupled ion pairs iri ZnSsMnS and (5) completely organizing 
all of crystal field theory in a form sufficiently general for application 
to any transition metal or rare earth ion in any crystal field. 

The main naw result of this work is that we now have a clear under¬ 
standing of the principles which must be applied in order to interpret 
optical spectra of ions in crystals. At the beginning of our work it 
appeared that crystal field theory could be used to calculate the spectra 
with the only complication being the occurrence of lattice vibration 
structure superimposed upon the electronic structure. We tacitly assumed 
that when a foreign cation is introduced into a crystal lattice it occupies 
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a substitutional cation site and that the surrounding anions have the 
same symmetry as in the pure crystal. This assumption has now been 
abandoned, and evidence for three types of exceptions to it has been 
found: 

1. In a non-cubic crystal having a cation site which lacks an 
inversion center, a foreign cation may move away from the substitutional 
site without altering the symmetry. The example we have investigated 

is corundum in which the foreign cation probably moves about 0,1 X 
away from the substitutional site along the threefold axis, 

2. At a site of high symmetry, such as the tetrahedral sites in 
ZnS, the foreign ion may destroy the local symmetry if its electronic 
structure causes large quartic terms to appear In the potential function. 
This can be interpreted as a second order Jahn-Teller effect. We have 

a little spectroscopic evidence for this type of distortion in the 
spectrum of Ni ++ in ZnO, ZnS, etc. and it fits in with magnetic evidence. 

3. For ions whose excited states are degenerate in a symmetrical 
configuration there exists the possibility of the Jahn-Teller effect 
and the consequent displacement of the equilibrium configuration toward 
a less symmetrical form. We have found numerous examples in which the 
electronically excited ion and its immediate surroundings do not have 
the symmetrical configuration of the ground state. It has become 

a general rule that the optical spectra do not reflect the symmetry one 
expects from the ground state except in a few cases. The most probable 
type of distortion appears to be tetragonal. 
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These results do not mean that crystal field theory is inapplicable; 
in fact without the comprehensive theoretical work which has been carried 
out in crystal field theory it would be impossible to interpret the 
spectra. The theoretical calculations must be carried out for the 
stable equilibrium configuration of each electronic state, and this 
may mean several different calculations for the same ion-host crystal 
combination. Since we usually do not know this configuration in advance 
the calculations must be programmed with great flexibility to permit 
making a number of trials until the spectrum can be fitted. We have 
only carried out a few rather crude approximations to this procedure, 
but have reached the point where it is now possible and profitable to 
do the more complete calculations. 

Further, We find that crystal field theory applied t6 ions in 
tetrahedral siteB accurately describes the observed spectra. This means 
that configuration mixing between <3d) n and (4j>)(3d) n " 1 , which arises 
from the lack of a center of symmetry, gives a negligible contribution 
to the relative energies. The cubic crystal field parameters Dq which 
are found are about half as large as those found for the same transition 
metal ion and ligands in octahedral coordination, which is a result predicted 
by the simple point charge model for the potential. We have found it necessary 
to go to a high order of perturbation theory in order to get an accurate 
fit of theory with experiment since the spin-orbit, crystal field, and 
electrostatic repulsion energies are all roughly comparable. The parameter B, 
which is related to the difference in the interlectronic repulsion energy 
between the angular momenta states P and F, is found to be depressed from 
the free-ion value. We have obtained an empirical correlation between 
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the parameter B for the ion in the crystal and the polarizability of 
i.mi livomm, v?e have also related the t <-*odve depression of rhe B values 
to a rough measure of the covalent character of the bonding and find 
that the covalency is in the order S* > Br” > Cl” > 0~ and iB the same 
for both cobalt and nickel bonded to the same ligands. 

The foundation of all the interpretations of spectra are the 
theoretical calculations. These have been carried out systematically 
so that any crystal field problem involving _d or f_ electrons (except 
for a few states of some £ configurations) can be set up on a computing 
machine and solved completely. It is no longer necessary to work in 
either of the extremes of coupling. This work is presented here in the 
form of tables of transformation coefficients and reduced matrix elements. 
The Wigner coefficients which have been calculated are not given here 
because they form a very bulky table, and they have already been published 
by another group. 

The persons principally involved were Dr. Santiago Polo, who carried 
out the bulk of the theoretical work, Dr. Herbert Weakliem who did 
spectroscopic, preparative and theoretical work. Dr. Donald McClure 
who carried out experimental and theoretical work and was responsible 
for the direction of the project, and Mr. Clyde Neil who capably assisted 
in all of the experimental work. Important assistance was rendered by 
Mr. Sherwood Skillman who programmed several lattice sums for the IBM 650, 
and Mr. Phillip Baltzer who made available his program for the energy 
2 

levels of the d system in tetragonal crystal fields. 






Optical Spectra of Transition Metal Ions in Corundum 
Donald S. McClure 
RCA Laboratories, Princeton, N. J. 

ABSTRACT 

The polarized optical spectra of the ions Ti +++ , V 1M , Cr 1 11 , Mh +++ , 
Co' 11 and Ni 111 in corundum single crystals have been studied at temper¬ 
atures from 4.2°X to 1200°K. A theory of the intensity based on the 
point charge model and p-d mixing has been developed and applied to the 
data with results in fair agreement with experiment. The effects of 
temperature show that the vibrational-electronic contribution to intensity 
is quite small at low temperature but may be appreciable at high tempera¬ 
tures. The crystal field parameters have been calculated as convergent 
lattice sums. The observed trigonal field parameter has the opposite 
sign from that calculated by the point charge model if the impurity ion 
is assumed to occupy an Al" t ~ H " ion position in the lattice, but has the 
same sign as calculated for an ion 0.1 8 displaced along the axis 
toward the empty octahedral site. Details of the spectra have been 
Interpreted as showing that the surroundings of an ion are distorted In 


some electronic states. 
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Optical Spectra of Transition-Metal Ions in Corundum 
Donald S. McClure 
RCA Laboratories, Princeton, N. J. 


1. Introduction 

A foreign cation at a substutional site in a crystal will be subject 
to a crystal field which is approximately that of the undisturbed crystal. 
In a cubic crystal the foreign ion either resides exactly on the cation 
site or the symmetry of the crystal is reduced in its vicinity. Electronic 
transitions at a cubic site are forbidden in transition-metal ions and 
their spectra are complicated by the presence of vibrational-electronic 
interaction. The crystal field must, however, have cubic symmetry except 
in cases where the radius ratios are quite different for host and guest 
cations, or when the Jahn-Teller effect occurs, In a crystal of the 
corundum type, there is no center of symmetry at the cation site and 
electronic transitions occurring in transition-metal ions are permitted 
without the intercession of vibrations. The spectra may be considerably 
simplified as a result. On the other hand, the impurity cation need not 
change the symmetry of the lattice if it moves away from the cation site. 
The normal adjustments of the positions of the iona in the locale of the 
impurity following the change in ionic radii will cause a shift of the 
impurity along the three-fold axis. The crystal field is therefore not 
as well defined as in cubic crystals. 

In this study of the spectra of the first transition-group ions in 
corundum, we have attempted first of all to understand the spectra in 
detail: this includes relative intensities and polarizations, band 
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energies and details of band features. The interpretation has been made 
by using the point charge model in various ways, but never to rely on its 
predictions of absolute magnitudes. By studying the entire series of ions, 
regularities are found and the level of confidence in specific results 
is made clearer. Some deductions relevant to the position of the ion 
in the lattice are then made. 

II. Experimental 

The crystals used were grown by the flame fusion method. The impurity 
content added to the starting material was on the order of 0.2 mole percent. 
The crystals were oriented by x-rays and were then cut to the required 
dimensions and polished with a diamond abrasive. The cut moat often 
used was one in which the c-axis lay in the plane perpendicular to the 
direction of propagation of the light. In this way, one could observe 
spectra in both polarizations of light with the same crystal, 

The spectra were taken with a Cary Model 14 spectrophotometer, and 
a Jarrell-Ash 3-meter grating (15000 lines/in) in the Ebert mounting. 

In each case, the spectra were polarized; by means of a Poucault prism 
mounted between the final slit and the chopper compartment in the Cary I4j 
and by means of a Wollaston double image prism abunted far enough in front 
of the slit of the grating spectrograph to produce separated Images of 
the crystal on the plate. 

All spectra were observed at temperatures of 77 and 298°K. In some 
cases noted later, it was necessary to use helium temperatures, and some¬ 
times to go to high temperatures, on the order of 1300°K. The high 
temperature spectra were done with a nichrome wound furnace which fits 
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into the Cary 14 cell compartment, A chromel-alumel thermocouple mounted 
in contact with the sample was used to measure temperatures, The tempera¬ 
ture was not controlled, but never fluctuated more than 10 degrees, which 
was not serious for the experiments being done. The infra-red optics of 
the Cary 14 were used so as to remove most of the black body radiation 
of the sample by the monochrometer. Since most of the measurements were 
in the visible or ultraviolet, the sample radiation was not a serious 
factor as high as 1000°C with this arrangement. 

Wavelengths were converted to wavenumbers by using Kaysers tables 
111. Theoretical Framework 

The crystal structure of corundum^'^ is illustrated in figs, la 
and lb, showing respectively a view parallel to c and perpendicular to c . 
Pairs of A? -atoms are stacked along the c-direction. These are surrounded 
by trigonally distorted octahedra of oxygen. The two A t -atoms of the 
pair share the equilateral triangle of 0-atoms lying between them. The 
AX -0 distance in this case is 1.969 5. The other triangle of each 
octahedron is larger and lies closer to the A jL -atom than does the first. 
This AH -0 distance is 1.856 Thus the site symmetry of the A J? atom 
would be Cj v if no other distortions occurred. However, the upper and 
lowar triangles are rotated away from the <T~ y planes by 2°8.5' in 
opposite directions, thus destroying the (f~ planes. This loss of 
symmetry is so slight however, that the site group C^ v may be retained 
for many purposes. The relations between the relevant symmetry groups 
are given in Table I. 
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The potential function of crystal field theory for An impurity ion 
in the crystal lattice may be calculated assuming that the ion occupies 
the A1 1,1 position. The symmetry of the surroundings of the ion are not 
al'ered if it moves along the c-axis away from the A1 ' 1 1 position. 

Potential constants were therefore calculated for two displaced positions: 

+ 0.1 8, toward empty octahedral site on c-axis and -0.1 8, the opposite 
direction. These calculations were made assuming no other atoms moved. 

The computations were made by summing directly over the lattice. 

This procedure is adequate for the second and higher powers, but not for 
the first, where the Ewald method should be used. The zeroth term is 
simply the Madelung energy per Al +++ ion. Since the calculations were 
made on a small machine (IBM 650) the summations over the lattice were 
not carried out to convergence of the second power term, but for this 
case we could compare the results to Bersohn's^. 

The lattice parameters and atomic positions used are given in Table II, 
The impurity atom was placed at OOw plus 000, OOz, or 00-z, where z 
corresponds to the 0.1 8 displacement. 

The potential function for an electron on the impurity ion may be 
written: 

V * B° + Bj r cos & + B° r 2 (3 cos 2 0 -1) + r 3 (5 cos 3 0-3 co«<9) 

+ B 3 r 3 sin 3 & cos 3 0 + C 3 r 3 sin 3 & sin 3 0 
+ r^ (35 cos 4 0 -30 cos 2 0 +3) + B^ r^ sin 3 & cos & cos 3 0 
+ C 3 r^ aln 3 & cos & sin3 P + higher power terms. 

Here, r, O and 0 are coordinates of an electron. The B constants are 
the ones which appear when the site symmetry is assumed to be C^ v , while 


(1) 
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the C constants are brought in when the true symmetry, C 3 , is used. 

Terms higher than the fourth power are not needed when we are only 
dealing with problems involving d-electrons. 

The potential constants are the following lattice sums over functions 
of the atomic coordinates when the point charge model is assumed: 
n° ne 

B o = r 

B? = cos a 

1 ir 


,0 


J 2 



—« <3 cos 2 & -1) 

4R J 

(5 C 0 S 3 & -3 cos e ) 
4 It 

~~T~ sin 3 & cos 3 $ 

8R 4 

| sin 3 <9 sin 3 ^ 


64R- 


16R J 
7 One 


(35 cos 4 ^ -30 cob 2 & + 3) 
sin 3 0 cos & cos 3 $ 
sin 3 & cos & sin 3 ^ 


( 2 ) 


Here, R, & , 0 are coordinates of a lattice site* n is either 1 or 3 
depending on whether or not the representative atom lies on the threefold 
axis passing through the impurity atom. Thus each term of B or C refers 
to a set of equivalent atoms. The orientation of the coordinate axes is 
important. Take the c-axis to be positive in the vertical direction. 

The impurity is located in the A1 site lying immediately below the empty 
octahedral site. The angle 9 is measured from the positive c-axis and 
the angle jt is measured counterclockwise around c from one of the upper 
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TABLE II 

Atomic positions and lattice parameters for corundum. Hexagonal coordinates. Space group 
D^j. (parameters and lattice constants from A. Newnham, ref. 3). 


a 

b 

c \ 



t 

Aluminum: 0 

0 

W 


0 

0 

0 


0 

0 

-W 

+ 

1/2 

2/3 

2/3 

+ hk <£ 

0 

0 

1/2+ w 


2/3 

1/3 

1/3 


0 

0 

1/2-wy 






Oxygen: 

V 

0 

0 > 







V 

V 

0 

| | 

0 

0 

0 



0 

V 

0 


1/3 

2/3 

2/3 

| + hlt'i 


V 

0 

1/2 

|*| 

2/3 

1/3 

1/3 



V 

V 

1/2 i 

1 ' 






0 

* 

1/2 / 






w 

- 0.102 

V 


0.306 







oxygens in the nearest neighbor shell. The quantity e was given the value 
14.3923 ev, R is expressed in 8. In the summations over the lattice, the 
oxygen contributions were multiplied by +1 and the aluminum contributions 
by -3/2, and these were then added. The summations were carried out to 
several values of R maj£ : 3.86, 5.81, 7.59 and 9.52 8. The convergence 
behavior of the sums could be Judged from these results. The third 
and fourth power terms changed very little after 5.81 8. The term Bj 
kept the same sign as for nearest neighbors but changed considerably 
between 7.59 and 9.52 Bersohn carried out the sum to about 20 & and 

showed that the sum had ceased to fluctuate appreciably. He obtained 
-3 3 

10.4209 ea Q while we obtained 13.335 ea‘ . Bersohn’s result has some 
experimental verification since it gives an antishielding factor for the 
A1 nucleus having the correct sign and approximately the correct magnitude. 
His value will be used In this paper. The values of the second power 
term in displaced positions were obtained from our results at 9.52 X 
multiplied by 10.4209/13.3350. The potential constants are given in 
Table III. 

The basis functions which have been used in this paper are the trans¬ 
forms of cubic field baslB functions shown in Table IV. Also shown in 
Table IV are their functional forms In a cartesian coordinate system whose 
z axiB is the 111 axis of a cubic coordinate system. The x-axls lies in 
the approximate reflection plane when the surroundings of the ion are 
supposed to have the symmetry Cj v . Therefore one partner E + of an E 
function is invariant under tf” v and the other partner, E”, changes sign. 
The same transformation matrix converts both the basis functions and the 
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coordinate system when the correlations of Table IV are followed. Four 
of the basis functions fall into two E representations of C^ v and one falls 
into A l of C 3y . 

The matrix elements of the potential between d-electron states only 
involve the second and fourth power potential terms. The energy of a d- 
electron is affected mainly by the cubic field, a fourth power term, and 
very little by the trigonal field which arises from a second and a fourth 
power term. The actual field is conveniently expressed as the sum of a 
cubic and axial field. In a strictly cubic field we would have the 
relation (B 3 ) c = - 20/"2(B°) c . The trigonal component of the fourth power 
term can thus be expressed as the deviation from this equality, and we 
may write: 

V. = [B° - (B?) ] (35 cos 4 6 - 30 cos 2 0 + 3)r 4 

4 4 4 c ( 3 ) 

+ B? [sin 3 & cos & cos 3(2^ - * (35 cos 4 & -30 cos 2 ^ + 3)] r 4 

4 207? 

where the first part is the trigonal (axial) field, the second the cubic 

B? 

field. For B? -(B?) we can substitute B? + - - . 

4 4 c 4 201/2 

The energy matrix of this potential has in it different diagonal 
elements of the trigonal field in each of the states , E and E. The t 
states are split apart by the trigonal field but the e state is only 
displaced. The separation of the e state from the center of gravity of 
the t states will be called 10 Dq even though a small part of this splitting 
id sue to the trigonal field. This traditional constant must therefore 
be defined in the following way: 
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TABLE IV 


Transformation from cubic to trigonal one-electron functions, and from cubic to trigonal coordinates. 
The t 2 functions and coordinates are transferred by the same matrix, T. The cubic functions trans- 
form like: t iy - xy, t„ - xz, t y2 - yz, e b » 2 ^(2z 2 - x 2 - y 2 ), e, - £(x 2 - y 2 ). 


cubic trigonal 



2 /s/6 - Us/6 

- 


z 


X 

T ■ 

0 - 1/vlT 

1 a/2 

; T • 

y 

- 

y 


1/VT 1/v/T 

l/v/T 


X 


z 


T 


cubic trigonal 






-~r(x 2 - y 2 - \/2 xz) 

A 

*** 

■ 

•> 

.M 

/ 4. 

~5=-(/s’ *y + y 2 ) 

t„ 


t s 


— ( Z 2-X x 2 --L y 2 ) 

r» 




VT 2 2 


Rop’n in 

c 3v 

E + 


cubic trigonal 


1 0 




e, 


-L (x 2 -y 2 + 2 N /5‘xz) 




= 


-M 


0 1 


e # 


e r 


-~L (xy-V^ yz) 

S 
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3 VT , 


' 27 


20 fT 


(4) 


The splitting between the t(A.) and t(E) levels will be given by 


6.0 1 , 160 


< B 4 + -X > r 
4 20f? 


t 2 + 6r 

and the small coupling term between the t(E) and e(E) states by 

v . = . iR s° 7 


20 VT 

The energy matrix then has the following form: 


a l (t) 

S (t) 
E (e) 


3 v - 4 Dq I 
■ j v - 4 Dq 
v 1 


v’ 

6 Dq 


(?) 


The three parameters in these equations are in principal directly measurable 
spectroscopically. Their formulas in terms of the point charge model are 
given in Table V. 

In this work we have avoided using the values of r n by taking ratios 

or by considering them as adjustable parameters. The values were actually 

( 6 ) 

computed from Watson's S.G.7. functions } but were used only as rough 
guides as to magnitudes of the potential constants. Their variation with 
atomic number is much too strong to result in realistic comparisons with 
experiment. For example, r^ changes from 7.07 a.u. in Ti* 11 to 2.38 e.ti. 
in Co 1 11 , while the values of Dq (rw?) change from 1980 to 1830 cm'*. 
Clearly some other very important factor is not taken into account by 
the point charge model. We may be able to believe the types of predictions 
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of this model which do not depend on v“, or we may be able to determine 

T "T 

empirically such ratios as r /r', and use them for valid predictions. 

In fact the results obtained by these procedures do seem to be believable. 

The absence of a center of symmetry in the corundum lattice makes it. 

possible for electric dipole d -■> d transitions to occur. This can 

come about because odd power terms may mix the even d-functions with odd 
functions. The latter may be either atomic p-functions or odd molecular 
orbitals on the ligands. We will be concerned mainly with one-electron 

transitions of the type t -> e. Then whatever the type of perturbation, 

there will be certain symmetry relations between the various possible 
transitions. These are expressed by the following relations for the 
symmetry C^. 

(t x' X 1 e x ) " ' (t x I y 1 V = ' (t y 1 y 1 V = ■ (t y I x I V * A 

(«,i »i. K )-(t.iny-« (t) 

(t »i' 1 V ‘ V' 1 V ■ 0 

The direction of the transition moment is given by x,y, or z, where these 
are in the trigonal coordinate system, x and y are perpendicular to the 
c-axis while z is parallel to it. The constants A, B and C represent 
the magnitude of the three types of transition moment. Their values in 
terms of the point charge model, and assuming that p-states are mixed 
with the d-orbitals, are given in Table VI. 

It was found that the selection rules for ere not accurately 
followed and it 1 b necessary to use the true site symmetry, C^. Then 
eight additional one-electron transitions are possible. Note that there 
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would be 18 when no symmetry Is present, and 16 are possible in Cy Three 
additional constants appear on reduction of symmetry from C^ v to Cy How¬ 
ever, since the departure from is slight, only one additional constant 
will be introduced, the other two being related to the first by means of 
the point charge model. Any error introduced by this procedure will only 
be an error in a rather small correction term. The constant introduced 
is the magnitude D of the one-electron transition (e y | x | t^). This is 
representative of transition moment integrals which change sign under 
reflection in the pseudo symmetry plane. The one-electron transitions are 
summarized in Table VI. 

The four constants which determine the intensity may, in principle, 

be obtained empirically from the spectra. They are expressed in terms of 
-J 

only two constants, r and r , when the point charge model is used. 

The optical spectra are expected to be determined mainly by selection 
rules for C 3v , shown in Table VII, but strictly by Cy also shown in this 
Table. The states are mainly determined by the cubic field but split into 
trigonal components as shown in Table I. 

The intensity formulas for several-electron atoms are derived from 
the results of Table V with the help of wave functions for the many electron 
case. The wave functions chosen are shown in Table VIII, and the intensity 

formulas are given in Table IX for one, two and three electrons. The high- 

4 1 

spin states of d have the same intensity formulae as d and the singlets 

of t 3 e in d^ have the same formulae as d 3 when the proper correlations are 

made, namely: ^(t 6 ) -> ^A 2 (t 3 )j 1 T 1 (t 5 e) -> S 2 (t Z e); 

1 T 2 <t 5 e) --> \ <t 2 e). 
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v.l».. ot canataat. to, p. d mixing ia poia, ah.,go apodal, 
A ■ j- Bf p - A B» p~ 

B ■ ■& efr * a Ab>,T 

c B«T a i£ Bf r T * | g3 f T 

D - A c/ pT 
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TABLE VII 

Selection and polarization rules for electric dipoie transitions in ions in corundum. 
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The notation for electronic states which will be used in this paper 
is illustrated in Table VIII. The cubic field parent such as 3 T. is 
given first, showing the multiplicity and the cubic field representation. 

2 

Following this is the strong field configuration such as t , and following 

this is the trigonal field representation, such as giving a complete 
3 2 

three part symbol Tj t A 2 - 

IV. Results: Intensities and Energies 

The questions to which we have obtained some answers are these: 

1) How well does the point charge model give the energies and 
intensities of the spectra? 

2) What is the explanation of the band shapes? 

3) Do ca tions which substitute in the corundum lattice occupy the 
Al +++ positions fairly closely, or do they move away appreciably 
in order to accomodate the different ionic radius? 

The data which has been collected are; integrated band intensities 
as a function of temperature, band positions as a function of temperature, 
trigonal splitting measured from spectra or from high temperature experi¬ 
ments, and details of band structure. Figures 2 through 7 give a general 
comparison of the spectra from Ti* 11 to Ni 111 , to d 7 . Table X gives 
the centroid positions and integrated intensities of the observed bands 
at 77°K, and the ground states. The bands are labeled according to the 
cubic field state from which they arise and their polarization. Chemical 
analyses were carried out on some of the samples so that absolute intensities 
have been given. Although the relative band intensity values are correct 
for a given ion, there are uncertainties in the absolute values vhlch will 
be discussed at the end of this section. In section V the spectral details 
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TADLE VIII 

Wave functions for two and three electrons for high-spin strong-field states of ions in 
corundum. A term such as e.r.i, means e s (I)t,(2)t s (3). The antisymmetrizing operator 
muse be applied to the functions given* 
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TABLE IX 

Intensity formulas for one, two and three d-eleccron ions in corundum. 


No. of 

Ground 

Excited 

Intensity 


d*electrons 

State 

State 

1 -L c 

lie | 

1 

*r 2 tE 

2 E e E 

2A 2 + 16D 2 

2C 2 + 18D 2 


J T J ‘ A 

2 E e E 

b 2 + d 2 

0 

2 

3 T lt 2 A 2 

3 T 2 et A 

0 

I8D 2 



3 T 2 et E 

A 2 + 8D 2 

0 



3 Tj ct A 

0 

2C 2 



3 r t et E 

A 2 + 8D 2 

0 


3 T! t 2 E 

3 T 2 et A 

b 2 + d 2 

0 



3 T 2 et E 

(A-BA/2) 2 + 9/2 D 2 

C 2 + 9/2 D 2 



3 T ie t A 

B 2 + 0 2 

0 



3 Tj et E 

(A + B//?) 2 + 25/2 D 2 

C 2 + 9/2 D 2 

3 

4 A 2 t } A 

*r 2 et 2 A 

0 

18D 2 



<T 2 et 2 E 

(A-B/y6) 2 + 9/2 D 2 

0 



% et 2 A 

0 

2C 2 



% et 2 E 

(A + B/y/5) 2 + 25/2 D 2 

0 
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observed at low temperatures are interpreted) in section VI the high 
temperature results are given and in section VIII the nature of the trigonal 
field is considered. Further discussion must b* carried out for individual 
ions. 



The polarized spectrum at 77° is shown in Fig. 2, The diffuse 

double humped band due to -•> is stronger when observed with 

light polarized parallel to the axis than with light polarized perpendicular 
to it. The polarization ratio was determined by integrating the absorption 
over frequency for each component. The ratio found at 77°K is « 3.2. 

Both the total intensity and the polarization ratio change by less than 
ten percent when the temperature is raised to 300°K. 

The observed polarization is compatible only with the assumption 
v > 0. In this case the E component of the 2 T 2 state lies below the 
component and the transition may occur in x,y or z directions. If were 
to lie lowest, the z-polarized absorption would be forbidden, in either 
Cj v or C^, according to Table VII, whereas it is observed to appear strongly. 
The theoretical polarization ratio in Cj v symmetry as obtained from Table 
IX is: 

III _ c 2 

— ’ 7 <»> 

Using a point charge model of the surrounding anions, the intensity ratio 

is calculated to be 8 when only the dipolar part of the anion field is 

used, and 17 when only the octupolar part is used. The calculation is thus 
qualitatively correct. When the symmetry is reduced to C^, the intensity 
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ratio becomes 4.4 and therefore is in nearly quantitative agreement. The 
addition of the dipole term as an adjustable parameter never results in 
a value lower than 4.4. 

The actual magnitude of v can be found by measuring the absorption 
spectrum at a series of temperatures high enough to populate the 
component of ^Tj. The polarization ratio would change as the absorption 
from became significant, and presumably could be calculated in terms 
of the intensity parameters. No such effect was observed up to 800°K, 
however; the value of v must therefore be greater than 500 cm”*. 

The excited state does not split in a trigonal field, yet a double 
peak is observed. This can be taken as evidence for a distorted state^*. 
The Jahn-Teller theorem^ shows that an E state is unstable when the 
three-fold axis is present. The potential is double valued, and the two 
peakB represent transitions ending on the two sheets of the potential 
surface. The separation of the peaks is 1850 cm"* at 77°K. This is larger 
than the trigonal splitting parameter measured for V +++ , Cr M1 or Co 1 ' 1 . 

In spite of the presence of a strong distortion of the upper state, the 
polarization seems to be unaffected. 

The centroid of the absorption band lies at 19380 cm”*, giving 
Dq = 1938 cm’*. 



The polarized spectrum of V 11 ' in corundum at 77°K is shown in Pig. 3. 

( 8 ) 

The spectrum has been discussed by Low' * and in great detail by Pryce 

(9) 3 2 

and Runciman . The ground state in an octahedral field is T^(t ), and 

3 3 

in corundum, the trigonal field forces the A^ component below E by an 





amount v. Second order spin-orbit coupling, mainly between the two 
3 -1 

trigonal components, causes the A 2 state to split by 8 cm , leaving 
the = 0 state lowest. At low resolution both components are observed 
at once and the selection rules are based on the orbital part of the 
3 A 2 state. 

There are two excited triplets which arise from the et configuration, 

3 T 2 an<1 3 T^. In the C^ v field they split into 3 A^ + 3 E and 3 A 2 + 3 E 
3 3 

respectively. The A 2 -> A^ transition is forbidden, but in the 

3 3 

transition becomes the permitted A -> A. These two states of et 

should give rise to absorption bands polarized |l C and _L C whose maxlmc 
are displaced from each other by v/2 by the trigonal field. The A component 
of each should be the lower one when v > 0 as in Ti'' 1 . 

The remaining triplet state is 3 A 2 , arising from the e 2 configuration. 
The intensity of this transition must arise from configuration interaction 

M a t 2 -■> e 2 transition is a two-electron change. Ballhausen^^ 

has suggested that the main source of this interaction is the coupling 
between ^(Aj) and 3 A 2 (A 2 ) in the trigonal field. The cubic field is 
not sufficient. 

Weak sharp lines from the doublets of the t 2 configuration are also 
observed. 

The experimental situation is as follows. The ^ state behaves as 
expected. The splitting between the absorption maxima in the two polariza¬ 
tions is 380 cm -1 apparently giving v = 760 cm -1 (but see later). The sign 
of v is therefore the same as found for Ti* * *. The magnitude is less 
than assumed by Pryce and Runciman, 1200 cm" 1 . The intensity ratio of 
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the oppositely polarized components of ^ Is 6, the largest anisotropy 
found in any of these spectra. The theoretical formula for this intensity 
ratio is 



( 10 ) 


The ratio using the point charge model and only the dipole term (so that 
r cancels out) is 16} and using only the octupole term it is 8. If we 
could ignore the interference between these terms it would be possible to 
fit the observed result. For the moment, we shall consider ourselves 
fortunate that the simple model gives qualitatively correct results, as 
was true for the case of Ti’ H " 1 \ We observe however that if the symmetry 
is used, eliminating D in (10), the polarization ratio is calculated 
to be much too high. 

3 

The T 2 band has some unexpected characteristics} namely, it is Just 
about as strong in the parallel spectrum as in the perpendicular} and there 
is no appreciable trigonal splitting. Furthermore there is very pronounced 
vibrational structure, and a band which appears clearly to be an electronic 
origin. 

The parallel spectrum can be explained by again using the correct 
site symmetry rather than C 3v , which would require zero intensity. 

The theoretical formula for the intensity ratio is: 


r 


18D‘ 

T7 


(ii) 
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The experimental value is 0.8 at 77°K. This value is given approximately 
by using the point charge model: it is 1.7 using octupole terms only and 
would be reduced by inclusion of the dipole term. 

The absence of trigonal splitting calls for a npw feature in the 
interpretation; considering the possibility that the upper state may be 
distorted into a configuration of low symmetry. This will be discussed 
in section V. It is related to the appearance of the vibrational structure. 
The polarization does not appear to be affected by the distortion, as we 
will see in section V. 

3 

The absorption in the region of the A 2 state also has unexpected 
features. It should be a weak transition appearing only in the parallel 
spectrum. What is observed are unresolved shoulders on a strongly rising 
background at 31,000 cm" 1 _L.C, 33,000 cm’ 1 |l C. Another shoulder occurs 
at 39,500 cm -1 and is more distinct X C than || C. These bands must not 
arise from pure d —•> d transitions. 

An important prediction of the intensity formulae is that in C 3v or 

the E components of and T^ should have the same Intensity of absorption. 

Assuming that I(T 2 E) is the intensity of perpendicular absorption in the 

region of the 3 T, state, and that I{T,E) is that in the region of 3 T,, 

2 I (T,E) 1 1 

their apparent ratio is t = %-j 

This discrepancy can only be explained by a special perturbation. The 
predicted equality of T^(E) and T^E) band intensities does not even depend 
on the use of the point charge model, since it is independent of the 
constants. It is a symmetry condition which depends for its validity upon 
the purity of the configurations. Configuration mixing cannot explain the 
large difference of intensity here. 
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The perturbation which seems most likely to operate here is the 
vibrational transfer of intensity. The Tj^A) band is predicted to be 8 
or 9 times stronger than T^E), so that a small fraction (about 15 percent) 
of the intensity of the A component appearing in the perpendicular spectrum 
would add an appreciable intensity increment. This can come about because 
the absolute value of the intensity Induced by the asymmetry of the 
environment is not very much higher than that capable of being induced 
by vibrations. Vibrations belonging to the E representation of would 
be responsible for transferring intensity from the parallel spectrum of 
TjA to the perpendicular spectrum. 

These observations on the spectrum of V* ** in corundum suggest that 
the most efficient way to compare the theoretical and observed intensities 
is first to use ratios which are independent of polarization then to 
compare the calculations to the polarization ratios which are least sensitive 
to vibrational perturbations. 

The two bands of the et configuration provide a comparison with the 
intensity formulae which is Independent of polarization. The total absorp¬ 
tion strength to either band is measured by A || + 2 A J_ , where A is 
the area under the absorption band component. The ratio of this quantity 
for the two bands is 8.8 and the theoretical formula is 

A || +2 AKTj) C 2 + A 2 + 8D 2 

A || + £ A^(T 2 ) = 9D 2 + A 2 + 8D 2 (12 > 

In this instance the octupole terms alone give a value of about 2.7, but 
the dipole term added in makes agreement with experiment possible. 
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I' 


A systematic comparison of all the intensity data with theory is given by 
the following list. The results using octupole terms alone and dipole 
terms alone, are given in the second and third columns. The dipole terms 

always spoil the agreement when the value of B° from Table III, and the 
_ “3 

S.C.F. value of r/r are added to the octupole terms. Therefore the 
0 _ "3 

quantity x = B^ r/r is considered as a parameter and its most favorable 
value is sought. This value and the resulting relative intensities are 
given in the third column. The unadjusted value of x is -3.04 ev/A°^. 


V 

t l e 

V T 2 

The agreement of the values for x = -j with experiment is nearly perfect 
if intensity transfer from T^A to T^E is permitted. Subtracting 5.5 from 
21 to give the observed 15.5 for IjA, we must add 5.5/2 to the TjE 
intensity giving 3.2, in good agreement with observation. The value of 
x is only small relative to the octupole term, and its difference from the 
unadjusted value is due to our uncertainty about the value of B° and to 
the failings of the point charge model. This model probably is much less 
reliable when comparisons of terms differing in the power of r are made 
than for terms with the same power. 

We can state some important conclusions at this point. The absorption 
we see in a given polarization is not entirely that due to a single electronic 


Exp. 

0.77 

1 

15.5 

4.4 

8.8 


Oct. Dipole 
Only Only 

1.7 0 

1 1 

8.1 16 


—2 
0.81 





29. 


state. This means that the separation of absorption maxima in two polariza 
tions may not be an accurate measure of the trigonal field. The effect 

3 3 

of the T^ A state on the perpendicular absorption in the region of T^ E 

is to shift the maximum toward lower energies by an amount equal to 

v/2 - hv g , where v g is the frequency of the depolarizing vibration. The 

apparent trigonal splitting could be less than the real splitting unless 

v/2 = hv e « The conclusion of this problem is discussed in section VI. 

A second conclusion is that we must use the correct site symmetry, 

Cj, in order to explain the intensity of the transitions forbidden in C^. 

3 3 

The transition A 2 -> T 2 A^^ forbidden in C^, is too strong relative 

3 3 

to A 2 -■> T 2 E to be explained by depolarizing vibrations, and only 

the use of the intensity formulas gives an adequate explanation. 

Cr +3 

The polarized spectrum at 77°K is shown in Fig. 4. The two bands of 
the transition from ^A 2 (t 3 ) to the quartets of the t 2 e configuration 
dominate the spectrum. The transition to the next higher configuration, 

2 

te , is forbidden in the approximation neglecting configuration interaction 
as it is a two electron change. It probably appears, however, as the weak 
band at 39,000 cm"*, partially hidden by stronger bands to the violet. 

2 

The splitting of the t e bands is easily measured, as the individual 
peaks can be located to within 50 cm”*, and are not overlapped by other 
absorption bands. The bands are separated by 800 cm”*; the T 2 by 
-1 2 

450 cm . The splitting of the te band is not accurately measurable, but 
appears to be around 500 cm"*. These splittings are somewhat higher than 
in the case of V 11 *, and do not match the theoretical expectations for 
first order splittings. These would be v/2, v/2 and v, respectively. 
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3 2 

Tanabe and Sugano showed that the bands of the t -■> t e 

transitions have approximately the correct polarizations when the dipole 
term of the crystal field potential is used^ 1 ^. The relative intensity 
of the two bands was not explainable in this way however. 

The two bands Of the et^ rmflonmUnn nf ry 1 1 1 »!•« *•« thccc 

of the et configuration of V 1 1 1 , but the intensity formulas are somewhat 
different. The systematic comparison of relative intensities as carried 
out for V 144 yields for Cr 111 the following table: 

Exp. 

T 2 A 0.27 

t 2 e 1 

TjA 2.3 

T^E 1.2 

t x n 2 2.i 

The additional adjustment in the last column is to reduce the value of 
3 

to one half of its value from Table III. Its only effect is to reduce 
the value of TjA, while adjusting x only affects TjA and TjE. The con¬ 
clusion again seems to be that the dipolar perturbation is unimportant 
relative to the octupole, and that it is much smaller than the computed 
value of Table Ill. 

The temperature effects discussed in section VI show that TjA 

3 

probably gains much of its intensity from vibrations. Thus may be 
even smaller than supposed. However, the other octupole terms are 
essential to explain the intensity ratios. 


Oct. Dipole * 

Only Only x= -.1 
0.55 0 0.30 

1 0 1 

2.6 4 4.5 

.52 1 0.57 

1.4 2.9 
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The spectrum of this ion shown in Fig. 5 consists of one strong band 

centered at 19,650 cm -1 . This is the spectrum expected for the ion in a 

5 3 5 2 2 

E(t e) ground state, since there is only one excited quintet, T£(t e ). 

The polarized spectrum should correspond to transitions from the degenerate 

^E state to componenents of ^T^« The state does not split in a trigonal 

field, nor is it split in first order by spin orbit interaction. The 

5 

components should be dtermined mainly by the trigonal field, giving 
below 5 E, separated by the amount v which is expected to be about 1000 cm" 1 . 

The polarized spectrum under these conditions should consist of a perpendicular 

band ( 5 E -■> 5 Aj) at about 19,000 cm' 1 , and a band of mixed polarization 

at about 20,000 cm -1 (^E -> 5 E). 

The spectrum actually observed is in fair agreement with these 
expectations. Two peaks are observed in opposite polarizations, but they 
are separated by 1900 cm" 1 , about twice the value expected. The lower 
energy peak at 18700 cm" 1 is predominantly in the perpendicular polarization, 
while the higher peak at 20,600 cm" 1 has mixed polarization. 

The intensity ratios which can be measured are the ratios of the two 
peaks summing over both polarizations, and the polarization ratios. If 
the first absorption region is identified as the Aj component and the 
second as the E component, the theoretical intensity ratio which does not 
depend on polarization is 

I (A) B 2 + P 2 

1 2A 2 + C 2 + 25D 2 


(13) 
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where the absorption is calculated for both partners of the ground state 
and both of the upper E state. A reliable experimental value of this 
quantity is difficult to obtain. Probably the best procedure is to choose 
an energy which appears to separate the absorption to the E component 
from that due to the A component. The value 19100 cm" 1 appears to be 
correctj and the quantity: 


A(<,L1) + 2A«, _L ) 

= °- 5 ' 


(14) 


where A (<,11) m eans the area under the parallel absorption band less 
than 19,100 cm" 1 . 

The calculated value of the ratio (14) is 0.42 using octupole terms 
alone, and is 0.2 using dipole terms alone. The octupole term gives 
nearly correct results and these are not improved upon by using the dipole 
term as an adjustable parameter. 

The polarization ratio for the transition to the E component of 


m 


1.25 


This ratio is given by the intensity theory as 


r 

i 


ii; 

x! 


2 

2 


(15) 


(16) 


Its experimental measurement is quite inaccurate because of the uncertainty 
associated with the cutoff position of 19100 cm" 1 . This is the same ratio 
which was evaluated to explain the polarization of the Ti 1 * * spectrum, and 
it was successful there. In this case the observed polarization ratio is 
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much less than the 4.4 given by the octupole terms. This result must be 
regarded as an even more drastic failure of the unmodified intensity theory 
to give correct polarization ratios than in the cases of vibrational 
transfer in V* * 1 and Cr 11 1 . 

The polarization ratio of the A component of is quite high and 
therefore in agreement with the selection rule forbidding parallel absorption 
to this component. 

Since the Mn* 11 ion has a ground state arising from the configuration 
t^e, it may well be unstable in the nearly cubic configuration of the 
A1 111 site in corundum. The distortion expected to result would be the 
tetragonal elongation of the octahedron observed in Mn 111 salts. If the 
ground state is so distorted, the upper state may be split under the 
resulting potential plus that due to the natural trigonal field. The 
resulting symmetry would be G fl if C 3v were the original site symmetry, 
and if the true site symmetry is considered. The energy matrix for 
a single d-electron in C fi symmetry is shown in Table XI. The matrix 

A 

elements for d quintets are one-electron matrix elements with signs 

4 

changed. For a tetragons], elongation of the octahedron in d , Dt < 0, 

Ds < 0, so that for this distortion e g falls below e b and t + < t” < t^ 
if Dt Ds > v/3. The transitions must therefore begin on e fl and end on 
t + , t~, t xy , with the intensity ratios: 

I II _9D^_ C 2 _9D^_ 

1 S 6D 2 + j (A + /2B) 2 ’ A 2 + 8D 2 ’ y- D Z + | (A - \/2B) 2 

(17) 

respectively. The agreement with experiment is only moderately good here. 
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The formulas predict a lower degree of polarization for the first band 
than the unmodified ones, and predict almost the same high polarization 
of the rest of the spectrum, which we previously found did not agree 
with observations. The agreement with a purely trigonal field seems 
therefore to be much better than the assumption that the Jahn-Teller 
distortion of the ground state persists in the upper state. The only 
way to obtain agreement with the intensity formulas is to have t fc the 
lowest component of the upper state, in which case the zeros of Table VI 
give the perpendicular polarization which is observed. Only a trigonal 
field can do this. The conclusion here is that the upper state of Mn* 11 
is subject to the unmodified trigonal field of corundum. The parentage 
2 3 2 3 

of the upper state is c ( T^)e ( h ^), and the absence of distortion may 
be connected with the fact that the e-electrons are coupled into a shell 
having nearly full cubic symmetry. In Sec. V we will present evidence 
that the ground state actually is distorted, however. The low polarization 
ratio of the E component might he explainable if we abandon the first 
order calculation and use perturbed wave functions appropriate to the 
distorted ground state. This has not been done however. 



The spectrum of this ion is due to transitions from the ^A^(t^e^) 
ground state to excited quartets, and is therefore very weak. The optical 
spectrum observed in crystals containing about 0.1 percent of Fe was of 

-1 4 

poor quality. The broad band expected to appear at 8000 cm ( T^) 

(as in aqueous solutions) could not be seen, but the bands at 25,600 
4 4 4 -1 

( A^ + E) and 18200 ( T 2 ) cm were observed. The upper states of these 
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TABLE XI 


One-electron matrices for c^symmetry; trigonal plus tetragonal field, 
(Same for C g within the d-configuration) 


A* 


-4Dq ♦ 2Ds - Dt ^2v/3 

^ 2 v /3 -4Dq ♦ | - Ds ♦ 4Dt 

2vVf6 -v'/4T 


2v’/ ^6 

-v'/fs 

6Dq - 2Ds - 6Dt 


-4Dq - v/3 - Ds + 4Dt -v’ 

.V* 6Dq + 2Ds - Dt 


cubic field basis functions;t « a> « b< 1 = 


nT? 
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bands arise principally from the state. The value of Dq obtained from 
these measurements is about 1650 cm 1 , but the bands are so broad that this 
value is quite uncertain. 

Co W 

The spectrum of this ion Fig. 6 shows that its ground state is the 
non-magnetic (t^) state. The other possible ground state is "^(t^e^), 
and the corresponding spectrum would then consist of one transition to "*E, 
probably a double band similar to that observed in Ti 111 . The observed 
spectrum, however, consists of two well separated bands plus a number of 
very weak ones. The strong bands are the singlets of the configuration 
t 5 e, 1 T 1 (15,560 cm'*) and ^ (22,980 cm' 1 ). They are split by the 
trigonal field by 360 and 370 cm" 1 (A > E) respectively, at 77°K. The 
first order splitting of each band is v/2, so that v = 720 cm' 1 , and is 

11 | 

therefore of the same sign and of about the same magnitude as for V 
There is a suggestion of another strong band at 30,000 cm' 1 near 
the sharply rising absorption (probably charge transfer absorption). 

This is in the right position to be the transition to (t^e^). 

Two weak bands, undoubtedly due to transitions to triplet states are 

seen at 18,880 cm' 1 and 19,800 cm' 1 . These are probably the lowest ^Tg 

4 2 3 3 4 2 

of t e and the accidentally degenerate pair T^ + Tj of t e , 

The ^ (t^e) band is doubled. It is believed that the low energy 

component in the parallel spectrum is a triplet state intensified by 

coupling to the singlet. 

2 

There is an exact correspondence between the states of the et 
configuration of Cr 11 1 and of the et^ configuration of Co H1 , The intensity 



formulas for Co + *”* are identical to those for Cr 1 1 1 when the correspondence 
4 T x (Cr 3-f ) - X T 2 (Co 3 ^) 4 T 2 (Cr 3+ ) - ^ (Co 3+ ) is established. 

The lower band of Co * * * is forbidden in the parallel polarization, as in 
V '** and Cr +++ , and is actually found to be very weak. The same type of 
intensity comparison as was made for Cr 111 can be made for Co 111 : 


Oct. Dipole 

Exp. Only Only X= 0.1 X= 0.05 


V 

0.31 

0.55 

0 

0.55 

0.55 

TjB 

1 

1 

0 

1 

1 

V 

1.2 

2.6 

4 

1.1 

1.9 

T 2 S 

0.69 

0.52 

1 

0.36 

0.44 

V T 1 

1.12 

1.5 

“ 

0.74 

1.11 

The results in this case also 

seem to be improved by adding in a 

. small dipole 

contribution x * 

0.05, but it 

is of the opposite sign from V 111 

and Cr +++ . 

Intensity transfer from T,,A to TjE amounting 

to 

35 percent is needed to 


explain these relative intensities. As for Cr +++ , the value of C 3 seems 
too large. 

Ni 444 " 

The spectrum of this ion, Fig. 7, consists of a weak sharp band at 
12300 cm two broad bands in the visible and some bands in the near 
ultraviolet which are too strong to be crystal field bands. The absolute 
intensity of the visible bands is about right for crystal field bands, 
and the infrared band could be an intercombination. It is difficult to 
assign the ground state of this ion from its spectrum. The choices are 
t 4 e 3 4 T X as in Co" 44- *" salts, or t^e a spin paired state which would be 
similar to the ground state of Mn'' 1 , 3 E. 
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In Che spin-unpaired state the lowest multiplet component of 
would be ( cu bic field) and the next multiplet component would be a 

G state about 500 cm * higher. At all temperatures below room temperature, 
only the transitions from E^ 2 would appear. The trigonal field is not 

expected to change the symmetry of the ground state. 

2 

In the spin paired state E, the spin-orbit representation in a 
cubic field is G, which splits in a trigonal field into E^ 2 + E 3/2’ The 
splitting between these components could only be on the order of 10 cm"* 


as it arises from the product of trigonal and spin-orbit perturbations. 
There is the additional possibility, as was the case for the Mn 1 * 1 
ground state, that the Jahn-Teller effect removes the degeneracy. 


Thqs..relevant. eel^tion^rule,? fsr-Scansifriofis .fry .?jre given in the 

lower right hand corner of Table VII. In the spin-unpaired case, transitions 
have to be either perpendicular, or unpolarizedj in the paired case, they 
will follow either E^ 2 or E 3 ^ 2 rules, or be a combination of the two, 
depending on the magnitude of the Jahn-Teller splitting. 


The observed polarizations are as follows: 
Band , 


Position, cm 

12.300 

16.300 
16,350 
16,700 
17,750 
17,850 
19,600 


Description 

line 

broad band max 

origin of group of lines 

broad band max 

line 

line 


Polarization 


broad band 



39. 


The occurrence of features polarized exclusively along c is compatible 
only with a ground state Such a level can arise only from the spin- 

paired ground state, but in order for the polarization to remain strongly 
marked, the ground state must be split, with lowest end not 

contributing to the absorption. If there is a distortion, however, the 
symmetry is no longer correct and the selection rules have to be re¬ 
examined. Since the new symmetry is probably C^, only Intensity calcula¬ 
tions can be used, but it is probable that these would give results 
closely resembling the symmetry rules of (as was found for the similar 
case of Mn 111 ). The assignments of these absorption bands is not at all 
uleckiT, however. 

The trigonal splitting parameter v might be measurable from the first 
band, which has its maximum at 16300 _L C, 16800 ||C. This splitting is 
probably v/2. Thus v = + 1000 cm"*, and has the same sign and about the 
same magnitude as in the other ions. 

In this section we have been concerned with spectral assignments, 
energy parameters and intensities. The results of the comparison of the 
intensities to the theory we have been using are perhaps rather unclear 
at this point/ so we will summarize them. We have made little use of the 
absolute intensities because their a N priori calculation involves too many 
numbers which an difficult to evaluate. 

The variations in total intensity from one ion to the next are also 
not well understood. A possible source of error in these numbers is the 
chance that the ion may be present in more than one valence state. The 
divalent ions are very weakly absorbing relative to the trivalent, and they 
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would not be detected in spectral absorbance measurements. A block of 
corundum 3 x 3 x 10 mm containing Ni*** was heated in a vacuum at 1600°C 
for 15 min. The yellow color of Ni 1 1 1 was entirely removed by this 
treatment, and the expected green or blue of Ni** did not appear. The 
spectrum of this colorless crystal showed some weak ultraviolet absorption, 
but none which could be ascribed to Ni**. Upon reheating this sample in 
an oxygen rich oxyhydrogen flame at 1600°C for 30 min, the yellow color 
of Ni 1 ' 1 was restored. In order to determine the spectrum of Ni** in 
corundum, a corundum slab was embedded in NiO and heated in a corundum crucible 

at 1700°C in a vacuum for 1 hour. The slab turned sky blue, and showed the 
++ -1 

optical spectrum of Ni , with Dq = 930 cm . The colored region of the 
crystal was, however, very close to the surface of the corundum slab, and 
the exact phase in which Ni** was present is in some doubt. In the case 
of Ti***, the samples all Contained what appeared to be a cloudy precipitate. 
This may be Ti 0^,, and of course the chemical analysis would not distinguish 
this from Ti***. This may account for the low apparent intensity. These 
few experiments show that considerable additional work will have to be done 
if reliable absolute intensities are to be obtained. 

The relative intensity values of bands of the same ion are experimentally 
in very little doubt. The interpretation by the point charge theory is 
qualitatively correct in every case so far observed. The degree of 
polarization is always less than the theory gives, and it appears that 
vibrational transfer of intensity can explain this. The ratio of intensity 
of two bands where measured by the quantity A || + 2 A _L is Independent of 
polarization, and a reasonable fit, using only octupole terms is found for 




V , Cr 


and Co' ‘ By adjusting one parameter, x = B“ r/r~, good 


agreement of all intensity ratios with experiment is possible for these 
three ions. The value of this parameter changes from -0.5 for V* f1 to 
-0.1 for Cr 44 * to +0.05 for Co""^ but the significance of this change, 
if any is not known. It might be interpreted as a change in ion position, 
in which case the ion moves from a region of negative to positive electric 
field as it is changed from V’ +++ to Co 1 1 1 . The intensity data all seem 
to require that the octupole perturbation of the ion by the lattice is 
greater than the dipole perturbation. 

V. Analysis of Spectral Details 

Up to this point we have discussed only the gross features of the 
spectra. At 77°K and below, fine structure is observed in all spectra 
except that of Ti 11 1 . This is generally identifiable as vibrational structure 

the vibrational intervals being characteristic of the corundum lattice. 

(121 

The best available assignment of these vibrations, due to Krishnan' ', is 
shown in Table XII. 

The introduction of an impurity atom reduces the symmetry from D^ d to 
Cj, therefore vibrations in the non-degenerate classes of D 3<J may all interact 
with each other, and vibrations within the two degenerate classes may 
interact with each other. There is therefore a possibility that the 
vibration frequencies around the impurity may be different from those of 
pure corundum. The extent of the change will depend on changes in force 
constants, which probably will not be great, and on the difference in 
mass betweel kjL and the impurity atom. The atomic weight of kJL is 
27 and of Ti through Cu 47.9 to 63.6. The frequencies involving motion 
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of the metal atom can change at most by the square root of the mass ratio, 
assuming the same force constants, or by a factor of from 1.3 to 1.53. 

If the motion involves the oxygens for the most part, the frequency change 
will be much less. 



The spectrum of Co shows a great deal of vibrational structure 
on the *Tj t“*e E peak, and many corundum vibrations may be identified. 
Most of these frequencies appear to be nearly the same as corundum 
frequencies. The correlation is shown in Fig. 8. Spectra were taken at 
4.2°K under high resolution, but no features other than those visible 
with the Cary 14 were found. 

Among the most definite features are the strong bands at 666 and 
742 cm” 1 from the origin. They correlate well with 642 e g and 751 a^, 
although the latter could be an overtone of 393. The 393 peak appears 
strongly and in association with two other peaks. It could be related 
to 375 6^ of pure corundum. A peak at 235 cm -1 could be the 244 a^ u 
corundum vibration. The equally strong 259 is not obviously related 
to any corundum frequency, although it might be one of the normally 
inactive ones. 

Although the bands are broad there are enough near coincidences 
so it seems safe to conclude that the local vibration frequencies around 
the impurity are not very different from those of the pure crystal. 

According to the identifications made in Fig. 8, not only non¬ 
degenerate vibrations, which become totally symmetric under C^, but also 
degenerate vibrations, are excited. The degenerate modes are non-totally 
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symmetric in C^, and their appearance in the spectrum can be ascribed 
either to vibrational-electronic, interaction or to a distorted equilibrium 
configuration in the upper state. There is no progression of eg modes as 
in the case of Cr 1 1 1 and V* H , so there is no obvious way to distinguish 
these two possibilities. The effect of temperature is the same as for 
Cr 1 11 j the Tj E band increases by a factor of 1.36 from room temperature 
to 1300°K (see Sec. VI). It appears therefore that vibrational-electronic 
coupling is no more important here than for Cr' 1 1 . If the excited state 
is distorted, the distortion must be smaller than in Cr 111 because only 
single quanta of the eg modes appear; and it must be different from Cr 11 * 
because many eg modes appear rather than one. In particular, the 194 eg 
mode does not appear in the Co* 11 spectrum, so the distortion must be 
orthogonal to this particular displacement. The fact that the 0-0 band 
of Co 111 is much stronger relative to the rest of the band than in the 
case of Cr 111 is another indication that its excited state is geometrically 
more similar to the ground state than for Cr* 11 . 

Figure 8 shows that there is a weak band in the parallel spectrum 
53 cm" 1 to the red of the strong 0-0 band. This could be a "hot band", 
but the helium temperature spectrum has not been sufficiently investigated 
to determine whether or not this is true. There is also a band +80 cm" 1 
from the 0-0 band. What these intervals are is not known as they are less 
than any of the reported vibration frequencies of corundum. It may be 
that -53 represents the origin of the parallel transition and that +80 is 
a 133 cm" 1 interval. An interval of about this magnitude (110 cm’ 1 ) is 
observed in Ni +++ . It is possible that it represents one of the inactive 
modes, the twisting motion of the octahedron about the 0., axis. This 
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mode might be easily excited because the transition is made allowed by 
3 

the parameter C 3 which measures the twist of the octahedron from the 
ideal C 3v geometry. It is a totally symmetric mode in C^, and is 
permitted to appear by the selection modes. 


These two ions have similar Tj bands as shown in Figs. 9 and 10. 

The most significant difference is that the sharp line structure a 

-1 | I I 

progression of 200 cm intervals, in the Cr band is highly polarized 
_L C, while in V* 11 the structure appears with about the same intensity 
in both polarizations. 

The origin band in Cr * 1 ' would appear in the perpendicular Bpectrum 

4 

according to selection rules for C^, because the lower component of T 2 

la *B, and the ^A 2 -•> transition is polarized x,y. But in the V 1 11 

3 3 

case the A component of is lower in the.trigonal field, and the 
3 3 ° 

polarization for A -■> A is z. The x,y absorption is Just as strong 

as the z, however, This cannot be explained on the basis of selection 
rules, and no vibrational perturbation is possible if this is the origin 
band. It was therefore important to establish whether or not these first 
bands (in the low temperature spectra) of V* *'* and Cr 1 1 are really 0-0 
bands. 

A series of spectra at various temperatures was obtained in order 
to locate the first "hot bands". If the apparent origins are actually 
0-1 bands of an enabling vibration, this vibration would appear as a 1-0 
band at temperatures high enough to populate its first quantum level in 
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the ground state. On heating the crystal from 77°K, the first feature 
to appear is a broad weak band 200 cm 1 from the apparent origin. It is 
most distinct at 120°K and can still be distinguished at 175°K. The 
spectrum becomes quite broad at high temperatures but there seems to be 
no evidence for the hot band which should appear at an interval roughly 
twice the 0-1 energy from the apparent origin. The apparent origin is 
probably a real origin, therefore, and the one hot band at -200 cm" 1 is 
one quantum of the progress ion-forming frequency in the ground state. 

| | | m 

The same experiment performed with Cr again showed the 200 cm 
hot band but it was less distinct. There was no evidence of any other 
band. Although this evidence is not conclusive it seems most likely that 
the apparent origins are really 0-0 bands in both cases. 

According to Table XII the totally symmetric frequencies of corundum 
are 578 and 751 cm" 1 . In the locale of an impurity the inactive vibrations 
and the two A 2u vibrations, 244 and 847, become totally symmetric. The 
progression frequency is not as close to any of these as it is to the 
194 cm" 1 Eu mode, and this is its probable assignment. If this is the 
correct assignment the appearance of such a progression is evidence that 
the excited state is distorted so as to destroy Cy 

The appearance of the unpolarized origin in V +++ is also evidence 
that the upper state has lost C^. Since the same progression appears 
in both Cr*** and V +11 they must both have the same distortion. If we 
knew the normal coordinate corresponding to the progression mode we could 
calculate the equilibrium configuration of the excited state from the 
intensity distribution on the progression. The distortion must be quite 
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large as the maximum intensity occurs in the eight or ninth member in V 1 '' 
and in the seventh in Cr f 1 1 . A crude application of theFranck-Condon 
principle gives the order of magnitude of the oxygen displacements as 
0.3 - 0.4 X. This is based on a harmonic oscillator model, however, and 
we must be cautious about applying such calculations to a system in which 
there are undoubtedly important linear, cubic and quartic contributions 
in the potential function from Jahn-Telier effects and Renner effects. 

These could make a vibrational progression stretch out without there 
being large displacements in the upper state. 

The origin of the distortion of the T^ state is undoubtedly electronic 
in nature and could be considered as a static Jahn-Teller effect. In the 
case of ions in corundum there is already a large trigonal potential which 
splits T 2 . The commonest type of Jahn-Teller splitting results in a 
tetragonal distortion of the octahedron. The axis of this distortion 
lies at 54.7° to the trigonal axis of the octahedron and may be directed 
along any of the three equivalent cubic axes. In corundum the octahedron 
is distorted, but thsre will still be three equivalent pseudo-cubic axes. 

In order to see if the observed spectra are consistent with an upper 
state having a tetragonal distortion, a calculation of the spectral 
energy and intensity under these conditions was made. 

The basis functions were chosen as cubic functions combined so as to 
transform properly under the symmetry plane containing the axis and 
the distortion axis, C^. This is only a pseudo-syranetry plane in the 
true site group of the impurity, but it should have considerable validity 
in energy calculations. These basis functions are those of Table XI and 
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the cubic functions are designated as in Table IV. Two- and three- 
electron functions were constructed and the matrices of trigonal plus 
tetragonal plus spin-orbit were then constructed. Some of these results 
are shown in Table XIII. The one-electron matrices shown in Table XI 
define the perturbations used. The eigenvectors of the matrix in Table 
XIII are essentially the basis functions when Dt > v, as is likely to be 
the case. These were used to calculate intensities. The intensity 
formulas and energy diagram are shewn in Fig. 11 for Dt > v. 

The sum of the parallel intensities is the same as is given in Table 
IX for the A -* A transition, and similarly for the perpendicular 
intensity- The tetragonal field merely redistributes the intensity. The 
discussion of Sec. IV is therefore not affected by these new considerations. 

2 2 

The intensities would agree with experiment if A ~ 6D , a reasonable 
requirement. The spin-orbit splitting of the lowest component must be 
small since we see only very closely spaced lines in the origin band. 

Fig. 9 shows that the origin band consists of at least two lines, 
differently polarised. The calculrted spin-orbit splitting of the lowest 
level, 0^, is a second order effect given approximately by: 

* 3 < “> ■ e 3 (0 > s -usr- 


(18) 
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TABLE XIII 


Wave functions for et configuration in C 2h (or C s ) symmetry, and matrix for tetragonal trigonal + 
spin orbit energy of 3 T 2 (et) (cross terms between trig, and tetrag. potentials are negligible). The 
roots of the analogous matrix for 4 T 2 (et 2 ) are very similar to the roots of this matrix. The cubic 
basis functions of Table IV are used here except that t + =1 ifl (t^ + t y!6 ), t“ -l/V?(t 1B ~ t yB ). 


3 T 2 (et) 


<!>\ - 1/2(V3 e 4 t + + e b t“) 
<f> 2 - l/2(\/Te 4 t~ + e b t + ) 

*3 ’ e b‘*y 


A" 3 T,(et) <j >4 * 1/2(vT e b t + - e 4 t”) A' 
A' - l/2( v /Te b f-e ll t + ) A" 

A A 


(The reflection plane is x - y, A'symmetric, A "anti-symmetric under reflection) 



a - (I +i)£/8 v * trigonal parameter 

b = i(/4 

$T a ^j(S o -1), etc. 


Dt = 4th power tetragonal parameter 
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Using 35Dt/4 = 1500 cm , £ = 140 cm a splitting 0.8 cm result 
This pattern would not be resolved by the Cary 14, nor was it cleai 
(91 

resolved by Pryce and Runciman in their higher resolution work, 
although they state that there appear to be several components at l 


The assumption that the level of V is tetragonally distorted 
is therefore consistent with experiment: the spin-orbit splitting is 
quenched and the intensity distribution is reasonable. The calculation 
for Cr 111 is the same as for V +1 * except that the trigonal; tetragonal and 
spin-orbit matrix elements must be given the opposite signs. The state ^ 
is the lowest one in this case. The intensity formulas are not the same 
(see Fig. 11), but the essential feature of is unchanged from the 
case of d^: it has zero moment jj C. This is a necessary requirement 
for agreement with experiment. Thus the similarities and differences of 
the two spectra are all explained by the same theory, and we conclude 
that the T 2 states of Cr +I ' and V* 4 " 4 are tetragonally distorted. 

There is one point which must be examined. The state ^ in d^ has 
appreciable spin-orbit splitting. The bands in the progression in Cr* 11 
are broader than those in V +++ , and this may be the reason why spin orbit 
structure is not resolved, or the spin orbit structure may be the reason 
for the broadening. 


Vibrational structure appears exclusively in the perpendicular 
ipectrum of Mn +++ near the beginning of the absorption band. (Fig. 5) 
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The first line, 16760 cm -1 probably marks the. electronic origin, although 
a more concentrated sample is needed in order to be sure of this. The 
intervals do not form a nearly regular progression as in V 1 * 1 and Cr* 1 1 . 

Each interval measured from the origin is close to one found in Co +++ ; 
but there are some remarkable omissions such as 666 and 742. 

j | | 

One feature of the Mr spectrum is not observed in any of the 
previous spectra of the series. A pair of narrow bands at 19130 || C, 

19340 _L C appear, Just barely resolved in the middle of the broad 
absorption band. These must mark the origin of the second component of 
the state, ^E. The origins of the A and E components are thus separated 
by about 2500 cm'*. The separation between the two lines, 210 cm'*, 
could be spin-orbit splitting, but the polarization and energy calculation 
should be carried out to check on the reasonableness of this suggestion. 

The large separation of the origins of the two components would be 
understandable if the upper and lower states have quite different con¬ 
figurations. The long buildup of vibrational progressions in the first 
band is evidence that the configurations are different. The origin band 
16760 cm'* is relatively weaker in Mn H1 than it is in V 111 , Cr +++ or 
C0 HI . The failure to observe an effective tetragonal splitting in the 
upper state of Mn +++ (sec. IV), due to the suspected Jahn-Teller effect 

in the ground state is now understandable. The upper state, belonging 
2 2 

to the configuration t e does not have a great tendency toward distortion 
as the e-shell is half filled (the parentage of this state is t^T^e^Ag)). 
In addition, the trigonal field forces the non-degenerate component lowest, 
reducing even more the tendency toward distortion. This lower component, 
which appears in the perpendicular spectrum must be a nearly pure "t^ hole" 
with the trigonal field around it quite undisturbed. 
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The long vibrational progression can only be accounted for if the 

ground state equilibrium position is appreciably displaced away from the 

trigonal geometry. Again, if we could identify all of the vibrations in 

the spectrum and knew their normal coordinates, we could determine what 

the ground state is like. It is possible to correlate most of the intervals 

with an e or e vibration of the corundum lattice: this of course would 
8 u 

establish the distortion as one which destroys the three-fold axis. This 
correlation is as follows; 

Observed interval: 224 380 492 694 885 1021 

Corundum vibration; 224(a u )375e u 450 g u 642e y 909eg ? 

Within the rather large limits which must be allowed, this correlation 
could be carried out in other ways. It is at least suggestive that the 
distortion is present and of the form required for a Jahn-Teller static 
distortion of the ground state. 

VI. Absorption Strength as a Function of Temperature 

By studying the strength of the optical absorption as a function of 
temperature, one can determine the contribution of vibrations to the 
absorption process, and can determine the excitation energy of low lying 
components of the ground multiplet if there are any. 



The spectrum of ruby at various temperatures is shown in Fig. 12. 

The total absorption strength to the t^e configuration is measured by 
A || + 2A J_ , where A || and A _L are the areas under the absorption bands 

observed parallel and perpendicular to the axis. These areas are taken 
4 4 

as the sum of those of the T 2 and Tj bands. Absorption strengths of 
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individual bands are measured similarly. The total absorption strength 
does not increase by more than two percent (the error of measurement) 
between 77 and 300°K. Above 400°K, the increase can be represented by 
the formula 

A (T) . . - ^D/T (19) 

A ( 6 ) - 1 + e 

where @p = 1050 + 30°K. This is nearly equal to the Debye temperature 
of the corundum lattice, 980°K. At 1000°K, the absorption strength has 
increased to about 1,35 times its value at 77°K. 

Individual bands have temperature dependences differing somewhat from 
each other. The strongest bands increase least with rising temperature. 

4 4 

The region of the forbidden transition, k^ -> Tj^), increases far 

more rapidly than the others. These effects are illustrated in Fig. 13, 
where the function 19 is also plotted. 

Two effects of temperature are to permit transitions from vibrations 
of the ground state, and to change the lattice constants. The excitation 
of certain vibrations will permit allowed transitions to appear in the 
opposite polarization from that given by the electronic selection rules 
and will cause violation of electronic selection rules of forbidden transi¬ 
tions. 

In the case of Cr 111 , we must distinguish between two possibilities; 
does the increase of absorption in the region of the partially forbidden 
4 

transition to T£ A.^ indicate, (1) vibrational intensification of this 

4 

transition, or (2) depolarization of the "allowed" transition to Tj E ? 

If the process is vibrational-electronic coupling to "enabling" vibrations, 
the forbidden transition will gain intensity from excited vibrations of the 
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ground state and will ultimately shift to the red with rising temperature 
by an amount roughly equal to the frequency of the enabling vibration. 

If it is the depolarization of the transition to the E state the parallel 
and perpendicular absorption should eventually have their maxima at the 
same wavelength, and the trigonal splitting should disappear. Super¬ 
imposed on these shifts would be the one due to expansion of the lattice. 

4 

Figure 12 shows that the band maxima remain about the same distance 
apart from 300°K to 1200°K. At the latter temperature all of the corundum 
vibrations are fully excited to at least the first quantum level, and more 
transitions will be taking place from vibrational levels than from the 

vibrationless level. It is clear therefore that the depolarizing process 

4 

is less important here than vibration-induced transitions to the T 2 
state. The forbidden transition almost triples its strength between 300 
and 1200°K and its strength becomes nearly as great as that of the allowed 
transitions. Thus the vibrational-electronic contribution to the intensity 
in the corundum lattice can become comparable to the contribution of the 
static field. 

The shift of the band maxima of the allowed transitions should 
reflect the changes of the crystal field with temperature. Theoretically, 
the band positions depend on Dq.In the cubic approximation, Dq varies 
with R - "* where R is the metal-oxygen bond length and changes in R, in 
turn, are related to the thermal expansion coefficient. In the equation 

Dq = _a 


(20) 
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where £ is a constant determined by the metal ion and the ligands , and R 
is the ion-ligand distance, the exponent n can be found from the temperature 
dependence of the band maxima. The energy of an electronic state is given 
by 

E = b Dq + C (21) 


where b and C depend on the state. We can take the logarithmic 
derivative of this equation and obtain 



Iff...: 9 .) 

is R 


n 


( 22 ) 


The temperature dependence of the band maxima give ^ 
and the thermal expansion coefficient gives i > and their ratio is 

-n according to Eqn. 22. The data of Fig. 13 for the band of Cr*++ 
give -gpj; = -1.1 cm’^/°C and since C = 0 this can be used directly. The 
linear thermal expansion coefficient of corundum is approximately 

- 8 x 10"^ °C. Assuming that the bond distance R is proportional 
to macroscopic dimensions, H , we find from Eq. 25, n = 7. This is in 
good agreement with expectations based on a point charge model (n = 5) or 
a point dipole model (n = 6). We cannot take the exact value seriously 
until more examples become available. 



The intensity of the two principal absorption bands of V +++ in 
corundum as a function of temperature are shown in Fig. 14. The points 
give the integrated band areas at each temperature divided by the area at 
77°K. Some of the data on which the figure is based is given in Table XIV. 
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3 

The integral of the ' band is somewhat uncertain at high temperatures 
because of overlapping by strong ultraviolet bands. It is relatively certain 
however, that the quantity A || +2A _[_ for this absorption does not change 
significantly with temperature. This is in spite of the fact that the 
parallel absorption decreases rapidly in maximum extinction coefficient 
as the temperature rises. The increase of band width and transfer of 
intensity to the perpendicular component makes up for the loss of maximum 

3 

absorption strength. The area ascribed to the T 2 band is actually the 
area measured from t’ . peak to longer wavelengths. This area is free of 
overlapping by the J T^ band. If the asymmetry of the band remains 
constant with temperature, the portion we have integrated will remain 
proportional to the entire band area. This assumption appears to be 
justified, as the band shape does not change below temperatures at which 
overlapping begins to become important, i.e. below 500°C. 

The total absorption to the first excited state of V 11 1 increases 
by a factor of 2.15 from 300° to 1000°K, while in Cr +++ this band increased 
by 1.36. Furthermore, the so-called forbidden component, 3 T 2 (Aj) in the 
z-polarized spectrum, which is 75 percent as strong as the E component 
at 77°K, becomes stronger rapidly and at about 700°K surpasses the intensity 
of the allowed component. At 1300°K it is about 18 percent more intense. 

This would not be expected if the mechanism of intensification were 
vibrational-electronic coupling or the depolarization process described 
previously. The two components could at most become equal in this way. 

The additional fact which must be considered in order to explain 
3 3 

these observations is that the E component of the ground T^ state lies 
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The integrated areas, A ||+ 2A j^> for the T^ and T, 


(expressed as optical density X wave number). I(T)/I(77°K) is the ratio of 
A| |+ 2Aj^ for T°K and 77°K. 
used for study of the bane 


-1 


1 1 | and vi are band maxima in cm . The sample 
id h ‘ 


had the same concentration as the one for 


but 7.5 times the path length. 



T°K 

I* 

2 II 

2 A 1 

I(T)/I(77) 

V | i 



77 

333 

429 

1.000 

17450 



300 

435 

508 

1.218 

17450 



375 

556 

590 

1.455 

17360 



425 

566 

570 

1.430 

17280 



475 

586 

615 

1.520 

17300 



530 

642 

651 

1.630 

17270 



575 

684 

723 

1.785 

17220 



675 

722 

726 

1.825 

17100 



780 

867 

774 

2.020 

16990 



875 

877 

875 

2.200 

16890 



975 

972 

867 

2.270 

16750 



1080 

1010 

908 

2.370 

16690 



1175 

1099 

940 

2.500 

16660 



1275 

1169 

991 

2.640 

16610 


3t i 

T°K 

A || 

A 1 

1(T)/I<77) 

V M 

V 1 


77 

2547 

725 

1.000 

24960 

25340 


300 

2550 

754 

1.01 

24880 



525 

' 2117 

1100 

1.08 

24560 



775 

2110 

1081 

1.07 

24440 
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only a little way above the A 2 ground component and becomes populated at 

moderate temperatures. Determination of this splitting provides a measure 

3 2 

of the parameter v. The intensity of the transition from Tj(E) t to 
3 

the T 2 state is theoretically quite high in the -3-polarization, and 
moderately high in the other. The intensity formulas for transitions 

3 

from the components of the ground state to those of the T^ upper state 

3 

are shown in Table IX. The ratio of the absorption strength to T 2 from 
the E state compared to that from the A 2 state is 

1(E) _ 3 B 2 + 2A 2 - 2 yY AB + C 2 + 31/2 D 2 


1(E) 

: i7a7 “ 


Using the dipolar potential of the lattice this ratio becomes 6 and the 
octupole part gives 20. Thus the observed increase of Intensity with 
temperature is qualitatively explained. 

3 2 

A value of the splitting between the E and A components of T(t ) may 
be obtained from a detailed analysis of the data. The populations of the 
two states at a temperature t are; 

N a (T) = N/(l + 2e" ^ /T ) 

N e (T) = 2N c‘ @/T /(l + 2e" ^ /T ) (24) 

where the factor 2 is the ratio of statistical weight of the E to the 
A state, 0 is the energy of E above A and N = + N^. Combining this 

information with the intensity formulas above, the intensity as a 
function of temperature may be derived: 


I (T) _ 1 + 2a e~ d/ T 

^ = 777F W7T ~ 


= p 


(25) 


where a is the limiting ratio of absorption strengths given in (23). 
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This formula must, however, be corrected for the vibrational intensifi¬ 
cation which was shown to operate for Cr* * *. We will not be much in error 
if we assume that the same functional dependence on temperature applies 
to both Gr 1 1 ' and V* 11 . Therefore the theoretical absorption strength of 
transitions beginning with t 2 of V 111 should contain the factor 

1 + e' , and that for transitions beginning with 3 T^ t 2 E should 

-(B + IB wt 

contain 1 + e ' ' * The ratio of these two factors will be 

denoted g: 


1 + e ' ( S + >)/T 



and in place of a in eq. 25, we will have ag. The value of & will be 
close to e p , so we simply substituted 9 = in g. Then the value 

of & can be obtained from the equation 

9 = 2.31 T log (27) 

The values of 9 for a = 8 and 8=6 are shown for a series of tempera¬ 
tures in Fig. 15. The values are moderately constant for a = 6 and give 
6 = 1380 + 100°K or 960 cm"*. This is also the value of v. 

The value of v obtained this way is somewhat larger than that 

3 

obtained from the splitting of the Tj^ band. The 960 value is probably 
more reliable, as the other valve may be affected by the vibrational 

transfer discussed in Sec. IV. The value obtained by Pryce and Runciman 

3 2 

from an analysis of the spin-orbit splitting of the T^ t, A 2 state is 
1200 cm’^. This value is affected by uncertainty in the value of the 
spin-orbit coupling factor and the spin-spin coupling correction. At this 
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stage, however, the differences between the values of v obtained in 
different ways seem minor, and it is gratifying that they all have the 
same sign and approximately the same magnitude. 

The value of a obtained from the temperature effect is equal to that 
calculated from the dipo e terms of the crystal field perturbation. This 
is in slight disagreement with the low temperature results, where a 
mixture of dipole and octupole terms was necessary. A value of a = 12 

is calculated for the case x = -i used in the low temperature.case. 

3 

The quantity corresponding to a calculated for the T^ band does 
not differ much from unity for any value of x. This is in agreement with 
the results shown in Fig. 14, where Che total absorption to this state is 
shown to be virtually independent of temperature. For this reason, it 
would not be possible to obtain a value of v from the temperature dependence 

3 

of the band, and because of the intensity transfer phenomenon, one 
cannot use the temperature dependence of the polarization ratio. 



The high temperature spectra of this ion are shown in Fig. 16. There 
5 2 4 

was a possibility that the (e t ) state could become populated 

sufficiently at 1000°K to give an absorption band due to the ^T^ -->^E 

transition. Although there is a strong absorption at about 25,000 cm"* 

this is probably a charge transfer band, as the 5 T 2 -> transition 

should appear at 18,000 - 19,000 cm" 1 . The temperature dependence of the 
*1^ band (6500 X) is the same as for the corresponding band of Cr 1 H \ 

VII. The Trigonal Field 

The study of six ions in the corundum lattice has made possible a 
comparison of the effects of the trigonal field on each. The separation 
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of band maxima of the components of a T^ or cubic field parent have been 
taken as the measure of the trigonal field. In Sec. IV it was shown that 
intensity transfer from a strong oppositely polarized component may shift 

j j | ^ 

the band maximum. In V the absorption where T^ etE is supposed to 
appear probably is largely the component which could absorb in the 
perpendicular direction in its e vibronic states. Since the e-vibration 
energies are not too different from the trigonal field energy and the 
component lies below the E, the depolarization may not have much effect 
on the position of the band maximum. The value of v obtained from the 
temperature dependance certainly agrees well with the directly measured 
band splitting. 

In the first broad band of V 111 the parallel and perpendicular 
absorption are mixed by the Jahn-Teller effect, as shown in Sec. V, The 
splitting of -90 cm 1 has nothing to do with the trigonal field. In the 
analogous band of Cr , however, it can be seen by reference to Fig. 11 
that the observed splitting may well be the trigonal splitting. The first 
component is polarized purely in the perpendicular direction and has the 
greatest intensity in this direction of any of the component states. The 
second component has the greatest parallel intensity, and is separated by 
v/3 from the lower one. The splitting of 450 cm’* - is therefore v/3, unless 
the depolarizing processes are appreciable. The value of v = 1350 cm"* 
is itv fairly good agreement with .the value obtained from the 800 cm"* 
splitting of the band, 1600 cm’*. 

The analysis of the Mn 111 spectrum given in Sec. V requires v = 1950 cm 
for this ion. Co' 11 and Ni 1 1 1 seem to have smaller values 720 and 1000 cm"* 
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The e:,ccited states of these two ions probably do not undergo as large 
Jahn-Teller distortions as V* ' 1 and Cr^*, so that the band splitting 
measurements give the correct values of v. 

The second order corrections to the splitting do not amount to more 
than ten percent, so they have been neglected. The best values of v for 
the different ions are given in Table XV. The values are not all the 
same, but they all have the same sign. 

TABLE XV 

Trigonal field parameters for ions in corundum. In cm"*. 

Ti V Cr Mn Co Ni 

v >0 960 1425 1950 720 1000 

The optical determination of the values of v has some uncertainties 
connected with it, as noted earlier. It does not seem possible that there 
is any error in the sign of v for all the ions, and it seems highly probable 
that the values of v are not the same for all the ions. 

The sign of v is given incorrectly by the point charge model if the 
S.C.F. values of r n are used. Another method of calculation is to find 

■x 

r from the value of Dq. For the constants given in Table V and the 
formulas for the constants in the normal and displaced positions, the 
results are as follows: 


Position 

7 . 

v (ev.) 

+0.1 

2.95 

—2 

-0.200 re + 0.494 

0.0 

2.58 

-0.119 r^ + 0.067 

-0.1 

2.39 

~ 

-0.083 re- 0.062 
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The value of Dq Is taken to be 0.223 ev or 1800 cm \ which is almost the 
normal value for all the ions in corundum, e is the charge assigned to 

the 0 atom. Since v= 0.1 ev or greater, the correct sign and magnitude 

~~2 

cannot be achieved in the normal position for any value of r e, but they 
can be achieved in the 0.1 position. The validity of this conclusion 
depends upon the assumption that the point charge model gives correct 
ratios of terms of the same power. The actual value of the displacement 
~~2 

would depend on the value taken for r e, but there are too many un¬ 
certainties connected with this choice to make further calculations profitable. 

The variations of v are seen to be quite possible theoretically if 
the ion position in the lattice changes a little from one ion to another. 

The fourth power term in v is particularly sensitive to displacement. Thus 
the positive sign of v may be caused by a displacement of the ion along 
the threefold axis, and the variations of its magnitude could be caused 
by differing araountf! of this displacement. 

The reason why a displacement of the impurity ion occurs has been 
suggested in the introduction. The increase of the ionic radius of the 
impurity over that of Al”^’* 7 *, 0.57 X is about 0.08 X. The increase in 
the force of repulsion due to this increase in radius can not be the same 
from the upper triangle as from the lower (Fig. 1), because these '•wo 
arrays of 0" ions are at different distances and they form triangles of 
differing size. In order to calculate in a simple way the magnitude and 
direction of the force acting on the impurity we will first assume that 
0 ions are not displaced and that they interact with the impurity 
ion by way of a repulsion potential of the form 
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V 


Ae 


% - 


( 28 ) 


where is the bond length, r n is the radius of rhe impurity and f* 
is a scale factor like the Rorn-Mayer constant. The increment of the 
potential upon changing an A1 * 11 ion into an impurity ion can be 
written down easily. When we do this for the corundum geometry we have: 

V' = A [exp (- p ( R u " r n " z cos 0 „)) + ex P <" ”■ ( R /. ‘ r n + z COB % ^ 
- exp (- ^ ( R u - r A ’ z co * ^„)> - ex P (- ^( R / - r A + z cos ))] 

(29) 


where z cos 0 u is the shortening of the bond length to an upper oxygen, 
and z is the displacement along C^> measured positive toward the empty 
octahedral site) r A is the radius of A1 1 1 1 . The force due to thiB potential 
term may be calculated directly: the resulting displacement would have to 
be calculated from the complete potential, and for this, other terms would 
have to be added, perhaps a harmonic restoring force, for example. In 
theundisplaced Al +++ posi.ion the force is: 

F z = - -37 - (e x P (^) - ex P (-cos © u exp (- ^±) + cos Q exp <- —■) 

(30) 

The force is positive (toward empty octahedral site) when r n > r A and when 
p > 0.238. The larger p the "softer" the force, and it is not unreasonabl 
to have 0 = ions with a fairly large value of p . 

This calculation makes it reasonable to expect ion displacement:, but 
can it account for the observed differences between the ions studied? The 
ionic radii of the trivalent ions are all very similar and it does not seem 
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possible Co explain Che large differences in Che values of v in these 

of Che metal ions as well as that of the 0~ ions. It is clear that the 
simple models do not explain these fine details and a more sophisticated 
approach will be necessary. 
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FIGURjs (JAfTlGNS 


1. Position of neighbors of Al-ion in corundum crystal. The lines shown 
outline an ideal octahedral site to illustrate the departures in the 
real crystal. 

2. Absorption spectrum of Ti* *' in corundum at 77°K. 

3. Absorption spectrum of V* 1 * in corundum at 77°K. 

4. Absorption spectrum of Cr 111 in corundum at 77°K. 

5. Absorption spectrum of Mn 111 in corundum at 77°K. 

6. Absorption spectrum of Co' ' 1 in corundum at 77°K 

7. Absorption spectrum of Ni 111 in corundum at 35°K. 

The high energy band of this spectrum may be a charge transfer band. 

8. Vibrational structure of the first strong band of Co 1 * 1 In corundum 
at 77°K. The vibrations of pure corundum and their assignments 
(ref. 12j sec text) are boxed in. 

9. Vibrational structure of the first strong absorption band of V +++ 
at about 5°K. 

10. Vibrational structure of the first strong absorption band of Cr 11 * 
at about 5°IC. The trace of structure in the || -spectrum is thought 
to be due to crystal imperfections or polarizer imperfections. 

11. Energy diagram and intensity formulae for 3 T 2 (et) of V 111 and S? 2 (et^) 
of Cr+++ in trigonal plus tetragonal field. The A’ and A" designate 
the representation in Cs. The spin splitting of 3A' Is given in the 
text. The diagrams given here apply to the case Dt » v »J“ . 

12. Absorption spectra of Cr * 11 in corundum at various temperatures (°K) 
Upper set X C, lower set || C. 

13. The increase in the area of absorption bands of Cr 1 '* with temperature; 
given as the ratio at temperature T°K to 298°K. 

14. Change of area of absorption bands of V* 11 with temperature. 

15. Plot of 0 / 2.31 vs T for a = 8 and a = 6 (see text). 

16. Absorption spectrum of Co 111 in corundum at various temperatures. 

The 1300°K spectrum is from a separate run. 
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Optical Spectra of Ni ++ ., Co + + , and Cu’^in Tetrahedral Sites in Crystals 
Herbert A. Weakliem 

RCA Laboratories, Princeton, New Jersey 
ABSTRACT 

The polarized absorption spectra of Ni* + and Co ++ in crystals of ZnO, 

ZnS, and CdS: Ni ++ in crystals of Ca 2 ZnCl^ and Cs 2 ZnBr^; and Cu ++ in ZnO 
have been measured at 4°K, 77°K, and room temperature. The spectra have 
been interpreted by the use of crystal field theory for the states 
of the (3d) n configuration acted on by a potential of predominately Tj 
symmetry. Certain details of the spectra are accounted for by smaller 
contributions from fields of lower symmetry, notably a C^ v potential 
contribution for the transition metal ions in ZnO. Crystal field, electro¬ 
static repulsion, and spin-orbit parameters have been obtained for all these 
cases. An empirical correlation between the electrostatic repulsion parameter,B, 
for the ions in the crystals and the ligand polarizibility has been obtained. 
Although the configuration mixing between the states of the configurations 
(3d) n and (3d) n "*' (4]>) has been found to give a negligible contribution 
to the calculated relative energies of the levels, it does partially 
explain the large values of the observed intensities of the transitions. 

The relative intensities of transitions between the spin-orbit componentr 
have been calculated by using models for both £ - d mixing and o-bondlng 
with the ligands. Neither model gives a quantitatively good account of 
the observed relative intensities and there is evidence, at least for the 
Ni ++ spectra, that the observed low energy, intense charge transfer absorptions 
play a predominant role in the intensity mechanisms of the d-d transitions. 

*The research reported in this paper has been sponsored by the Electronics Research 
Directorate of the Air Force Cambridge Research Laboratories, Office of Aerospace 
Research under Contract AF19 (604)-5541. 






INTRODUCTION 


A study of the absorption spectra of tetrahedralLy coordinated 
transition metal ions in crystals has been undertaken. The spectra of Ni ++ 
in ZnO, ZnS, CdS, CSgZnCl^ and Cs^^Br^; Co ++ in the former three crystals; 
and Cu ++ in ZnO are discussed in this paper. The room temperature spectra 
of powdered solid solutions of NiO and CoO in ZnO 1- ; the visible region of the 
spectrum of Ca 2 (Zn ^Ni ^Cl^ 2 ; and the 77°K axial spectrum of a crystal 
of ZnOiCu have been measured previously. In each of these cases, the spectrum 
has been interpreted as being due to the transition metal ion in tetrahedral 

coordination. Our results on the crystal spectra are in agreement with this 
45 

conclusion. Crystal field theory has been used successfully for detailed 
interpretations of the results. Many of the foraal aspects of the theory for 

it 

ions in octahedral coordination, for which an extensive literature exists , 

are directly applicable to the case of tetrahedral coordination since the point 

groups 0 and T d are isomorphous. One non-trivial difference in the theory 

applied to these two cases is that the potential of a tetrahedral array of 

charges has an odd power term proportional to xys. This term causes no 

further splitting of the degenerate d n levels, but may be effective in 

removing the parity classification of the eigenfunctions, for example by 

causing the ground state (3d) n even parity configuration to be mixed 

with the higher energy (3d) n “ 1 (4£) odd parity configuration. The effect 

of such mixing on the relative energies of the (3d) n levels and on the 

intensities of transitions between the energy levels has been examined. 

Corrections to the relative energies are shown to be negligibly small, which 

agrees with recent corrected estimates by Low and Weger^. 

The odd parity contributions to the wave functions due to £-d 

mixing qualitatively account for the observed intensities of the bands, but 

fails to give a good account of the relative strengths of individual spin-orbit 

*The review articles cited in References 4 and 5 are devoted to crystal field 
theory and contain most of the pertinent literature references. 
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components. Other contributions to the intensity of the transitions, including 
a molecular orbital (MO) model and mixing with charge transfer states, are 
discussed, but a detailed model which gives agreement with experiment has not 
been developed. 

Spin-orbit, electron repulsion, and crystal field parameters have been 
obtained and are correlated with certain physical properties of the host lattice. 

It has been found that the electron repulsion parameter B depends linearly on 
the ligand polarizability. The values of X obtained for Co** are larger 
than those found from spin resonance^’® and magnetic susceptibility^’^ 
studies on similar materials. 

Each section of the paper is usually subdivided into several parts. 

The first part, which treats the general case, is followed by three sections 
specifically concerned with Cu ++ , with Ni*"* , and with Co** and the last 
part deals with general conclusions. 

EXPERIMENTAL 

The ZnO crystals were grown from a lead fluoride melt by slow cooling 11 . 

The transition metal ions were incorporated in the crystals by adding about 1% 
of the oxide or carbonate to the initial mixture of zinc oxide and lead fluoride. 

The crystals usually grow in the form of thin hexagonal plates with the hexagonal 
c: axis normal to the plate face. Other crystal habits were sometimes obtained, 
but the crystals were of too poor a quality for optical studies. Spectra 
obtained with the plate face of the crystals oriented normal to the beam are 
polarized with Ei and perpendicular to ^c, The other polarization was obtained 
by mounting a section of one of the more dilute crystals so that c was perpendicular 
to the direction of propagation of the light, In this orientation, the crystals 
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were at most 0.2 mm in one direction, but by carefully masking with a slotted 
copper disc and placing the crystal at the exact focus of the spectrophotometer, 
satisfactory spectra were obtained. 

ZnS crystals grown from the vapor phase which were predominately hexagonal, 
but which had numerous stacking faults, were doped with the transition metal ions 
by vapor phase diffusion of the metal sulfide in the presence of a sulfur 
atmosphere. The same method was used to obtain crystals of CdS:Ni and CdS:Co. 

Pure CdS crystals obtained from the Eagle-Picher Co, and crystals grown from 
the vapor phase** in these laboratories were both used. Nickel diffused quite 
rapidly into these crystals at 800°C, as much as 0.1 percent being incorporated 
in a few hours. The incorporation of cobalt proceeded more slowly at this 
temperature at^d heating for several days was required in order to obtain concentrat: 
as high as 0.1%, 

Cs 0 ZnX )i (where X - = Cl" or Br") were prepared 2 by crystallization from 
an aqueous solution containing stoichiometric amounts of CsX and ZnX,, 
acidified by the addition of HX. To the molten salt under an atmosphere of 
HX was added the desired amount of either NiXj or CoXg. The melt contained in a 
sealed, evacuated Vycor tube was then crystallized by use of the Bridgeman 
technique. The resulting large single crystal fragments were readily cleaved 
to give suitable samples for the measurement of spectra. 

The crystal spectra in the range 0.2 p to 2.5 p were obtained with a 
Cary model 14M spectrophotometer. Low temperature spectra were measured 
with the crystals attached with silver paste to the cold finger of a brass, 

Grown by C. W. Struck, RCA Princeton. ————— 

** Grown by L. A. Barton, RCA Princeton. 
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two chamber cryostat equipped with quartz windows. A Foucault prism which was 
placed in the optical path between the beam divider and the entrace slit was 
used to obtain the polarized spectra. The degree of polarization was better 
than 95% in the visible region of the spectrum. The resolution generally 
obtained was 10 - 20 cm \ Infrared spectra were measured for some of the 
crystals at room temperature with the Perkin-Elmer models 21 and 221 spectro¬ 
photometers. 

CRYSTAL STRUCTURES 


The crystal structure of the host lattice, and in particular the site 
symmetry at the transition metal ion position, are of central Importance to 

the understanding of the spectra. The crystals of ZnO, ZnS, and CdS used all 
12 4 

have the wurtzite structure , space group P 6^ m c - C fiv with the atoms occupying 
two sets of special positions (b) whose coordinates are (OOv) and which have 
the point symmetry 3 m = C^. It is convenient to think of the wurtzite 
structure as consisting of a hexagonal close-packed lattice of anions with 
cations occupying one half of the tetrahedral holes (and also forming a 
hexagonal close-packed lattice). Each atom is tetrahedrally surrounded by 
four atoms of the opposite type. There are also octahedral holes in the anion 
lattice which are normally unoccupied in the pure compound and Which lie on 


lines parallel to the hexagonal c axis. These normally empty sites form a 

B i 

rather open channel along the c axis. The ideal structure has c/a - (y-)2, 

v = 0 for M, and v • £ for X which makes the M-X bond length equal to 

(^-) 2**a. An atom at the center of an octahedral hole in the anion lattice has 


a slightly larger bond distance, equal to - 


This position is also surrounded 


by an octahedron of M atoms, three of which are at a distance (#)2 a and the other 
1 8 
17 — 

three are at a distance ( 2^)2 a. The actual value of c/a is slightly less than 


the ideal value for the three compounds being discussed. The M-X, tetrahedra 
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are therefore not regular; the M - X distance which is parallel to is 
slightly shorter than the other three equivalent distances. Taking 8 to be 
the difference between the two bond lengths (negative if the unique distance 
is the smaller one) and calling the three equal bond lengths d, one finds 

- 8/d * 0.0154, 0, and 0.008 for ZnO, ZnS, and CdS, respectively. 

13 14 

Both Cs^ZnCl^ and Cs 2 ZnBr^ are isostructural with K 2 S0^ which is 
orthorhombic, space group P n m a : ^ our m0 ^- ecules P er unit cell. 

The structure contains discrete tetrahedral ZnX^ ions, lying in mirror planes, 
and surrounded by eight Cs + ions. The tetrahedra can be described by taking a 
regular tetrahedron and tilting one bond in the mirror plane by about 5° away 
from the other bond which lies in the mirror plane. There is also a difference 
of about 0.2 X between the pairs of bond lengths lying in the mirror plane and 
the pair perpendicular to the mirror plane for the bromide, whereas the lengths 
are all nearly equal for the chloride. The true symmetry of the anion is m = C g . 

RESULTS 

The unpolarized absorption spectra of Ni* + in ZnO, ZnS, and CdS at 77°K 

o 4*4“ 

and 298 K are shown in Fig. 1. The unpolarized absorption spectra of Co in 

the same compounds and at the same temperatures are shown in Fig. 2. In each 

of these cases, the same crystal was used for the two different temperatures 

and the curves were drawn to the same scale, the base lines being arbitrarily 

shifted in order to minimize overlap of the two curves. It is apparent upon 

comparison of the spectra for these two temperatures that the integrated intensity 

of each of the bands is not very temperature dependent. Some of the spectra 

have been measured at temperatures up to 1000°K and it was found that the 

integrated intensity of the bands decreases slightly with increasing temperature. 

This clearly indicates that the transitions are not vibrationally induced^. 

The f numbers (oscillator strengths) of the principal room temperature bands 

for four cases are given in Table I, together with the wave numbers of the band 
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maxima and the assignment of the upper level of the transitions. The f numbers 
were calculated from the formula 

f = [0.431 x 10- 8 ] p( v > d ^ 

where p(v) is the optical density as a function of frequency v measured in cm 1 
c is the concentration of absorbing species in moles/liter and d is the thickness 
in cm. The integral was estimated by a graphical Integration of the area under 
the curve for each band, The f numbers are factors of 10 to 100 higher than 
typical values found for the same ions situated in a centrosymmetrlc ligand 
environment. This is one piece of evidence that the ions substitute for a host 
lattice cation at its regular lattice site and are therefore tetrahedrally 
coordinated. This is also supported by the assignment of the transitions which 
is discussed later. 

Spectra of Ni + + 

The room temperature spectra of Ni ++ all show three rather broad 
bands of increasing intensity in the direction of higher energy of absorption. 

The near infrared-visible band la the most Intense of the three and is split 
into two components in each case. With the exception of CdS:Ni, a steep 
absorption edge several thousand wave numbers less than the band gap energy of 
the pure host crystal is observed. The spectra of very dilute crystals have a 
very Intense absorption band to the violet end of the spectrum. This is shown 
in the lowest curve of Fig. lb for ZnS:Ni which is the 77°K spectrum of a 35 p 
thick crystal containing about 0.005% Ni ++ . The absorption band has its 
intensity maximum at 25,000 cm 1 and at higher absorption energies only the very 
steep absorption edge associated with the valence band to conduction band 
transition of ZnS is observed. The f number of this band is estimated to be 
nearly unity. This is clearly not a transition of the d-d type and is associated 






with an electronic process involving both the nickel ion and its ligand 


environment which will be referred to as a charge transfer process. There 
is some evidence, discussed in the Intensity section, that this transition 
plays a role in providing the observed intensity of the weaker, lower energy, 
predominately d-d transitions. 

A similar charge transfer process in dilute CdSsNi was not observed;as 
Fig. lc shows. However the lower band gap energy of pure CdS may be too 
small to permit one to observe the transition. The lower curve of Fig. Id 
shows the charge transfer transition, which has a peak maximum at 22,500 cm \ 
for ZnO:Ni. The band gap energy for ZnO is about 4100 cm" 1 less than that for 
ZnS whereas the charge transfer transition for these crystals containing Ni + + 
is only 2500 cm" 1 less. Similar high intensity, high energy transitions have 
been observed for Ni* + in crystals of Cs^ZnCl^ and CSgZnBr^ but are not shown 
here and will be the subject of a later paper. We only remark here that, unlika 
the case of Ni ++ in the wurtzite lattices, there is a whole series of interns, 
anisotropic absorption bands in the near uv for nickel in the cesium tetra- 
halozincates. Figure 3 shows the d-d transitions at room temperature for Ni ++ 
in the latter compounds. 

4 + 

Polarized spectra of the near infrared-visible bands of Ni in the wurtzite 
lattices measured at liquid helium temperatures are shown in Fig. 4. These 
curves were traced directly from the spectrophotometer recordings and are therefore 
plotted on a wavelength scale, unlike the previous figures referred to which are 
plotted on an energy scale. Each curve is labeled to show the orientation of the 
electric vector of the light beam with respect to the crystallographic c axis 
of the host crystal. The numbers on the figures are the wave numbers of the 
indicated peak maxima and more prominent shoulders and were obtained from the 
wave lengths by the use of Kayser's tables 1 **. The two 


polarizations are nearly 
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identical for ZnS:Ni and CdS:Ni. Since the splitting of the degenerate d n 
levels of the free-ion is caused primarily by the electric field of the nearest 
neighbor ligands in the crystal.it is concluded that either the environment of 
the nickel ion is cubic or local anisotropies are randomly oriented throughout 
the crystal. The observed anisotropy for ZnO:Ni, shown in Fig. 4a, has been 
interpreted as being due to a small trigonal field superimposed on a predominately 
tetrahedral field. The spectrum with J. is identical with that obtained by 
directing the light beam along the c^ axis. In the latter case, it has been 
pointed out that both JS and are J. c, but since this spectrum has been found 
to be Identical with the £ X spectrum, we conclude that the transitions are 
electric dipole ones. 

Spectra of Co** 

The spectra of Co* + in the wurtzite lattices, shown in Fig. 2, are somewhat 
similar to those of Ni ++ . There are three broad bands (the weakest one in 
the region of 3500 cm' 1 is not shown for ZnO:Co), the higher energy one being 
more intense, and the visible band is split into two main components at room 
temperature. The band at about 3500 cm' 1 la quite difficult to detect for the 
concentrations used and has only been observed at room temperature. The 77°K 
spectra show that the infrared band centered near 60C0 cm 1 is aplit into four 
main components Whereas the visible band is split Into a larger number of bands, 
some of which are rather narrow. 

A careful search using dilute crystals failed to reveal the presence of a 
charge transfer transition at energies less than the band gap energy. If the 
charge transfer process roughly corresponds to exciting an electron, from the 
valence band into the ground state of the transition metal ion, the larger value 
for the second ionization potential of nickel compared to that of cobalt indicates 
that the charge transfer process for cobalt would be at a higher energy than 
for nickel. This can be described as a greater trap depth for nickel compared to 
cobalt, and has been discussed by Klasens 1 ^ in a paper dealing with Impurity 


centers in ZnS. 
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The 4°K polarized spectra of Co** are shown in Fig. 5. The spectra are 

presented on a wavelength scale with the wave numbers of the individual 

4.4- 

components shown, just as for the corresponding Ni spectra. Many of the 
individual components are quite narrow, especially for ZnO:Co, which is not 
usually the case for spin allowed d-d transitions. The sharper components 
appear on the low energy side of the band whereas the higher energy components 
are usually broad. The spectra of the sulfides are found to be isotropic but 
the spectrum of ZnO:Co is highly anisotropic. The ^ 1 ^spectrum is 
identical with the axial spectrum for the latter case. The conclusions 
discussed in the last paragraph of the section on Ni** spectra apply here to Co"** 
as well. 

The uppermost curve in Fig. 2b, covering only part of the spectrum, was 
taken with a yery concentrated sample and clearly shows some of the weajc doublet 
statep which are not discernible in the dilute samples. Concentrated, thick 
•staples of ZnO:Co failed to reveal any weak transitions to the low energy side 
of main visible band, however a few extra details on the high energy side were seen. 

The Spectrum of ZnOsCu 

The spectrum taken with the light beam directed along the c^ axis is shown 
in Fig. 6a and the polarized spectrum measured at liquid helium temperature 

-1 9 

in the 6000 cm region is shown in Fig. 6b. The energy level diagram for d 

In a tetrahedral field including spin-orbit splitting is also shown in Fig. 6a. 

Mulliken's notation extended to the double group^ is used for the representations; 

is doubly degenerate and G has fourfold degeneracy. The broad weak absorption 
observed at -on temperature, which corresponds to the t 2 - e transition in the 
cubic field approximation, breaks up into a pattern of narrow absorption bands 
when the temperature is reduced to that of liquid helium. When spin-orbit 
splitting is included, two transitions are expected in T^ symmetry. The spin-orbit 






parameter, £ , for the Cu ++ free-lon is 829 cm *, therefore the lowest energy 
transition ( 2 T 2 > - G ( 2 T 2 > is expected to be found in the region 1000 - 1200 

There is a very strong absorption due to a combination band of the transverse 
and longitudinal optical mode lattice vibrations of pure ZnO in this energy 
range. Therefore, the crystal field transition cannot easily be observed and 
has not been observed. 

The G levels are split by fields of lower symmetry into two Kramers doublets. 
The G ( 2 e) level has only a second order splitting due to the combined effect 
of spin-orbit coupling and a low symmetry field. The pair of sharp lines to the 
low energy side of the spectrum have been interpreted as the two trigonal field 
components of G (*£). The lowest energy component is almost entirely absertt in 
the parallel spectrum. The pattern of lines between 5780 cm * and 6100 cm" 1 
is repeated at intervals of 430 cm’ 1 . The pattern broadens and becomes quite 
weak rapidly,appearing only asa weak structureless band after the first Interval. 
The interval of 430 cm" 1 is close to the reported value for the transverse optical 
mode frequency of 414 cm 1 for ZnO 18 . 

DISCUSSION 


The Hamiltonian appropriate to the problem being considered is given by 

* =£{--£- v* • f~ + ( &+• ( Vi} + X £ 


where the sum runs over all the n electrons of the central atom and V is the 

c 

crystal field potential acting on the ion. The first two terms uniformly shift 
all levels of the degenerate d n configuration, the next three terms, which 
partially remove the degeneracy, are arranged in increasing order of magnitude 
of their effect on the level scheme under consideration. The terms are all of 
such relative magnitudes that it has been necessary to diagonalize the matrix of 34* 
for the d n configuration in order to obtain a good fit between the calculated 
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energies and the observed spectra. This approach has been discussed by Liehr 
and Ballhausen 19 who give the matrices of ^4 for both weak field and strong 
2 

field basis functions for the d configuration. The potential given by the point 
charge model for a regular tetrahedron oriented so that the two-fold axes of 
the tetrahedron lie along the x a axes is 


V 

c 


. Sftif- cf > 3 cf*>f [ V 2 



( 2 ) 


-*M* 'r* 4 I LV+<fe )2 <V + V 4 >j 


where d is the distance from the central atom to the point charge e and the are 

20 

the spherical harmonics with relative phase specified by Condon and Shortley . 

In order to examine the effect of a field of trigonal symmetry, it is convenient 
to choose the axis as the *' axis. x n this case the cubic potential term 
for the tetrahedron, which is oriented so that one apex lies on the positive a* 
axis and a second lies in the x's' plane with direction cosines 



If the crystal field potential is neglected, the states of a d n configuration 
may be classified by the angular momenta functions |SLJM >, terms with different 
S, L, and J values having in general different energies. The functions may be 
expressed as linear combinations of Slater determinant product functions of the 

20 

individual d-electrons by standard techniques for the coupling of angular momenta . 
The effect of V c is to partially remove the (2J + 1) degeneracy of the SL terms 
and the resulting functions may be classified according to the irreducible 
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representations of the symmetry group of V c < The new functions jsUTT >» where r 
denotes the irreducible representation of the group and r *abels one of the set 
of basis functions, are degenerate with respect t:o r and can be expressed as 
linear combinations of functions with different H for each SLJ. Using these 
functions as a basis, the matrix of the Hamiltonian (1) may be factored into 
a number smaller of matrices, one for each r. The first contribution to the 
splitting of free-ion a n levels from (2) or (3) comes from the fourth power 
potential term V^. The effect of the odd power term is to mix in states 
of higher energy^ odd parity configurations with the Btates of the d n configuration. 
Low and Weger 6 have used second order perturbation theory in order to estimate 
the relative shift of the E and T 2 tetrahedral crystal field levels of (3d) 6 
because of mixing with (3d) 5 (4p) and (3d) 5 (4f) via They find, in a 
corrected estimate for divalent iron, a relative depression of the levels 
amounting to less than 10% of the energy separation given by the (3d) 6 
configuration alone. Calculations we have performed for (3d) 6 and (3d)^ 
indicate lass than a 1% admixture of the odd parity functions with (3d) n 
functions. This gives a negligible correction to the energiea, but is capable 
of giving a large contribution to the observed intensity of transitions between 
the energy levels. Discussion of the latter point la reserved for the aectlon 
concerning intensities. 

The term of Eq. (2) la Just - 4/9 times the corresponding term for an 
octahedron with the seme bond distance d. Taking the splitting of a single 
d-electron in a tetrahedral field to be 10 Dq, one expects to find Dq values 
for an ion in tetrahedral coordination to be about one half the value for the 
ion in octahedral coordination with the same ligands. Since the sign of the 
splitting parameter is opposite for the two cases, the level schemes will be 
Inverted for T^ symmetry compared to 0, and thus in the e orbitals of the 
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d-electron have lower energy than the orbitals. We will be concerned entirely 
with d n configurations of more than a half filled shelly which may be correlated 
with the simpler d^" n configuration if the particles are taken to be holes 
and therefore have a positive charge. Dq will then be positive for d 7 , 


The matrix elements of for states of d 7 , which have half-integral J 
values, were calculated and are given in Appendix 1. The basis functions 
|d 7 SUIT > needed to calculate the matrix elements of were obtained by 
applying the projection operators for the T^ double group to |JM > functions. 

This procedure gives the correct linear combination of functions with different 

V e 1 

M values for every set of SU values. The matrix of ) - is diagonal with 

U X iJ 

respect to S and L and is independent of Mg and and the matrix of ) 1^‘ is 
diagonal with respect to J and independent of M = Mg+ M^. The matrices of these 
quantities calculated for the free-ion functions were added to the matrix of 
to give the matrix of , which is factorable into two 9x9 and one 21 x 21 
matrices. For the d 8 configuration, the matrices given by Liehr and Ballhausen 


Using trial parameters obtained from the spectra by the application of 
perturbation theory, successive approximations were made by diagonalizing the 
matrices using an IBM 650 (which also gave the corresponding eigenfunctions) 
until a good fit was obtained between the observed spectra and the calculated 
energy levels. The final set of parameters obtained are given in Table II and 
the predicted energy levels are drawn on each of the spectra of Co** and Nl** , 
shown in Figs. 1, 2, and 3, together with the identification of the component 
states. The notation used for the irreducible representations is Mulliken's 


extended to Include the double group. The energy levels are designated by 

Tj(L), where L is the orbital symbol for the level in zero crystal 
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field, 2S+1 is the spin multiplicity of the level, denotes the irreducible 

representation for the orbitai. part, and stands for the irreducible representation 
2g+ i 3 

of a spin-orbit component of the Tj level. For example, T 3 - T^(F) 

stands for the spin-orbit component transforming like which arises from 

3 

coupling the spin triplet, orbital T^ symmetry states of the (free-ion) F 
level. Of course, states of the same representation are in general mixed, and 
such a designation for an eigenstate is meant to imply that the true eigenfunction 
is composed mostly of the designated function. The character table for 
which is given in Table III shows how certain elementary functions transform, 
and lists only the "unbarred" operations of the single group. The characters 
of the three extra classes E, 8C^, and 6§^ of "barfed" operations needed when 
the double group is considered are simply the negatives of the characters for 
the "unbarred" classes. The degeneracy associated with an Irreducible representation 
is given by the character for the identity operation E, The last column in Table 111 
gives the representations in C 3v which are correlated with the T d representations 
in the first column when the symmetry is reduced from T^ to C^. 

The appropriate first order descriptions of the energy levels for the d® 
and d^ configurations are shown in Figs. 7 and 8 respectively. The diagrams 
are drawn in the spirit of successive perturbations of decreasing magnitude 
on the free-ion levels which are given on the extreme left of the figures. The 

energies of the free-ion levels relative to the lowest one (set equal to zero) 

21 

are given in terms of the Racah parameters B and C, where B - p 2 - 5 and 

20 

C - 35 F^ and the F k are the Slater-Condon-Shortley integrals . Application 
of the crystal field perturbation splits the free-ion levels into levels 
characterized by Irreducible representations of the group whose energies in units 
of Dq relative to the unperturbed level are shown in parentheses. Spin-orbit 
Interaction gives further splittings into components whose representations are 
shown and whose energies relative to the crystal field levels are given in terms 
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of the free-lon spin-orbit parameter, Certain aspects of the two diagrams 
are common. There are only two states of the highest spin multiplicity, P and F, 
separated by 15B with F being lowest. The P state is not split in the crystal 
field and the splitting of the F state for d^ is inverted from that of d®. 

The spectra of N-» ++ and Co* + are thus expected to bear some resemblance to each 
other which is, in fact, the case. 

The observed splitting and relative intensities of some of the bands in the 

two polarizations of the ions in ZnO are explained by a small trigonal field 

in the presence of the predominant cubic field. It was pointed out in the section 

on Crystal Structures that the departure from T^ symmetry is due to the fact that 

the bond distance lying along the hexagonal ^ axis is slightly shorter than the 
other three equivalent distances. This gives rise to an axial field and the 
extra potential term which must be added to the cubic potential term is 


■■ft* 


where ^ is the measure of departure from T. symmetry defined earlier. The second 
a a 

and fourth power terms in Eq. (4) may cause additional splitting of the degenerate 
cubic field levels. 

In the discussion of the spectra of the individual ions which follows, it 
will be necessary to know the selection and polarization rules for electric 
dipole transitions. These are given in Table IV for the groups T^ and C^. 

There are no polarization rules for T d and a cross simply means that the 
transition is allowed. The electric dipole moment for T d transforms like Tg, 
whereas in C^ v the moment in the z direction transforms like and the moment 
In the x y plane transforms like E. 
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The Spectrum of ZnO:Cu 

We will discuss the simplest electronic system first; that of Cu ++ which 

g 

has the (3d) configuration. Since this configuration may be treated as one 
hole in a filled 3d shell, the problem is formally equivalent to that of one 
electron and one need only change the sign of Dq and of the spin-orbit matrix 
in order to obtain the correct cubic field matrices. The last column in Table III 
shows that the cubic G^) level splits into a pair of doubly degenerate levels 
having the representations Ej ^ 2 and ^ 3/2 ant * t *' e ground state becomes E ^ 2 
C^. According to the selection rules given in Table IV, the transition 
E 1/2 - E forbidden In the parallel spectrum, whereas Ejy 2 - E ^ 2 Is 

allowed in both polarizations. We therefore identify the line at 5784 cm ” 1 
as E 3 ^ 2 and the line at 5821 cm " 1 as E^. Therefore, the trigonal field 
splitting of the upper state is 37 cm 1 . In a cubic field, the energy of the 
transition is 10 Dq + £, and taking £ - 800 cm 1 one finds Dq z 500 cm 1 . 

22 

The energy matrices for a trigonal field given by Llehr were used in 

order to obtain the second order splitting, A - E [E ^ 2 ( 2 E) ] - E [E,^ ( 2 E)J. 

The result, expressed in terms of the potential given by Eq. (4),la 

A zsLi 11 (sL\ _ 120 _L_ 

a d \ 7 ' d 2 f 63 Njr/jM ; lODq , 


/ r V 

where the radial average value of the quantity r n /d n . We can estimate 

the magnitude of the second term by using the point charge definition of 

2 1). 

Dq » which is appropriate for this case, in order to eliminate 

(\) • Taking 5 : 800 cm 1 , and ( rj- ) = - 0.015, the value for pure ZnO, 

we find the contribution to A from the second term - 2.14 (~-)£ - 26 cm’ 1 . 

It is more difficult to estimate the first term since we do not know the 

relative values of and «> for the ion in the crystal. We have, however, 

d d 

estimated this term on the basis of three different assumptions, where in each 
case we have used the values (4-) : - 0.015, £ : 800 cm’ 1 , d : 3.7 a.u., and 
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Dq - 500 cm a) The Hartree-Fock wave functions for Cu^"*" recently tabulated 

by Watson 2 ^ were used to calculate the moments <r”> . We find <rS> ■ 1.028 

... 1 2 

and therefore obtain A : 8 cm . b) We have obtained <r > from the spectrum 
of ZnO:Ni ( vide infra) and have scaled it to the case of by use of the 

factor obtained from Watson's H-F functions and find < v^> ~ 0.685. The splitting 
in this case is calculated to be A = 14 cm" 1 , c) The H-F value for < r^>= 2.498 


and for d s 3.7 a.u. , s 5.60 . Substituting this value for 

d 2/ ‘ V' 

and eliminating <—r) as before, we obtain A s - 130 cm . With the exception 
d 

of the last case, the sign of A is given correctly, albeit the values are low, 

The magnitude of A is also rather high for the last estimate. In this connection, 
it is interesting to note that the H-F value of cr^ gives a good account of 
the strain-Induced splitting of the 2 E state of MgO:Cr*^, whereas the ratio 


together with the observed Dq leads to too large a splitting 2 ^. 
This is similar to the result we find here, and in addition the H-F ratio gives 
the sign incorrectly. Although one should not place too much reliance on such 

calculations of potential constants because of the inadequacy of the point charge 

25 

model, a fact which has been thoroughly discussed in some recent papers , we 
think they do lend support to the assignment and indicate that the trigonal 
field is reasonably well accounted for by the local distortion of the tetrahedron 
in ZnO. > 

The relative values of the intensities of the three observed transitions 
were calculated by assuming the mechanism to be configuration mixing of (4j>) with 
( 3 d) via the third power potential term in Eq. ( 3 ). Since the 3p functions 
transform like T ? in T d only the t 2 ( 3 d) orbitalB are mixed. The orbital 
functions in trigonal orientation given by Liehr 22 were used to calculate the 
tj — e moments. Liehr's notation has been used for the orbital functions, and 
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, 1/2 

< fc 2a W 6 b > = ' < C 2c Id £ a > : " ( 15 > ^ ~ l £ 

1 1/2 

< c 2a W e a > = < fc 2b ^ % > = ' ( 15 > ^ +1 i? 

, 1/2 

< c 2b U *. > = - < E 2 c l£l % > = C 15 > “>' & 


where co* = < pJ^l^iV^d " E 4 p^ ‘ We n0W UBe Che zero order functlons 
appropriate to E ^ 2 ( 2 T 2 ) and G fa) in order to obtain the relative intensities. 
It is unnecessary to use the first order functions or the odd power moment 
contributions from the terms in Eq. (4), since this would only give additional 
terms reduced by a factor of ( ^ ). We have furthermore assumed that the 
Kramers degeneracy is effectively lifted for the partner wave functions, and 
have taken the intensity to be proportional to the sum of the squares of the 
four individual moments for a given transition. The intensities are then 
found to be 

E 1/2 < 2 V ~ E l/2 (2e) : K 45 (a) ’ )2 ' I 

\, 2 fa) = K fj ( a )’) 2 1 

- E 3/2 < 2 E) = K (to') 2 1 

where the constant of proportonality between the intensity and the square of 
the moment is very nearly equal for the three transitions and has been called K, 
The observed intensity ratios are fa) ||: ( 2 E) J_ : E^ ( 2 E) J_ - 

1:2.5:5.3 compared with the calculated ratios of 1:2:3. The agreement is 
satisfactory and further supports the assignment. 

The spectrum therefore consists of the two trigonal field components of 

2 

G( E) followed by a vibrational progression characteristic of the host lattice 
whose frequency is perhaps slightly modified by the Cu ++ impurity. It may be seen 



19 . 


in Fig. 6b that the pair of sharp 0-0 bands appear with approximately the same 
splitting in 0-1 and although the lower energy band at 6223 cm * is not 
completely polarized, it is stronger in the perpendicular spectrum whereas the 
6257 cm 1 band is stronger in the parallel spectrum. This leaves unexplained 
the nature of the weak bands lying in the energy range 5900 - 6100 cm"* (and 
the repeated pattern at 430 cm" 1 higher energy). They may be due to local 
vibrational modes coupled with the electronic state; however we have no 
quantitative information concerning this point. 

Spectra of Ni* + 

The three main bands observed at room temperature for each of the Ni** 

3 3 

spectra are transitions from the ground state T^(F) to the excited states T^(F), 

3 3 

A^(F), and T^P). The symmetry selection rules for electric dipole tranaitions 
given in Table IV are that Tj^ - A 2 E, and T 2 are allowed. At low temperatures 
when all the electrons are in the lowest spin-orbit state, the allowed transitions 
to excited spin-orbit components are given by the selection rule A^ - T 2 only. 

Fig. 7 shows that each of the excited states mentioned above has a Tj level lowest 
3 3 

in spin-orbit coupling. Both the T 2 (F) and T^P) bands are observed to split 
at low temperatures with the lowest energy component of each band being the 
most intense and usually being the sharpest one. This is particularly evident 
for ZnO:Ni. It has not bean possible to pick out the four components of 3 T 2 (F); 
however, diagonalization of the matrix of shows that three components are 
grouped rather closely in energy with a total energy separation of approximately 
X./2 with the fourth component, A 2 lying about 3X/4 above the mean of the other 
three. The diagonalization also shows that the total splitting of the 3 T^(P) 
band is very nearly 2\, rather than the first-order result of 3A. More than 

3 

four components are found in the spectral region of T^(P) and the total energy 
separation is 1000 - 1500 cm 1 . The appearance of components with representations 
other than T 2 violates the selection rule A^ - T 2 only, which may be accomplished 
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by fields of lower symmetry and/or coupling with vibrational modes which cause 
the forbidden states to be mixed with allowed ones in higher order, Under these 
circumstances, it is difficult to obtain an accurate value of and the value 
X s - 250 cm * applies reasonably well to all the spectra. Table V gives 
the comparison between the calculated (using the parameters of Table II) and 
observed spectra of Ni ++ in the three wurtzite-type crystals at 4°K. The four 
observed bands which are correlated with the four calculated components of 
3 

Tj,(P) may be identified by reference to Fig. 4. 

In all three spectra, there are reasonably intense broad bands to the high 
3 

energy side of the T^(P) components which are not identified. Some of the 
extra bands which appear are probably due to transitions to different vibrational 
levels of the upper state. Detailed examination of the effect of coupling of 
degenerate electronic and vibrational (localized or impurity modified lattice 
mode) motions have not been undertaken for these cases. However, some calculations 

have been made of the effect on the energy levels in Tj symmetry of adding 

* 

potential fields having tetragonal and trigonal symmetry . It should be possible 
to use these results in treating the dynamical case by the addition of a restoring 
potential to the problem. The low symmetry field can be considered to arise 
in part from certain of the normal vibrational modes of the complex whose 
instantaneous spatial configuration contributes to the low symmetry potential field. 
The result for both types of distortion is that the ground state remains a non- 

3 

degenerate Aj state. The A 1 component of T^P) has an energy shift in the 
opposite sense to the other components of ^T^(P) and the ground state for both 

types of distortion. Similarly, both components of 1 D have quite different 
energy slopes than the ground state. The implication is that at the value of the 


* Performed in collaboration with P. K. Haltzer and D. S. McClure, RCA Princeton. 
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distortion coordinate where the ground level has a minimum, these excited levels 
will not have a minimum and the excited state potential surfaces may in fact 
have a steep energy slope. Therefore a transition from the lowest vibrational 
level to these excited states will cover a rather large number of excited state 
vibrational levels, giving rise either to a broad band or to a series of excited 
vibrational bands. This argument is identical to that given by Orgel^ who 
considered only the difference between the slopes of potential surfaces which are 
functions of a totally symmetric mode in order to account for the relative 

breadths of absorption bands. These considerations qualitatively account for the 

3 

fact that the higher energy bands in the T^P) spectral region and the bands 
identified as the E and components of are observed to be broad. A consideratioi 
of only a totally symmetric mode predicts in fact that the transitions to states 
of should be sharp, since the configuration of the ground state and S) are 
nearly the same and consequently the energy of these two levels plotted v£ Dq 
have nearly the same slope. 

The rather high intensity of the spin-forbidden transitions Aj - ^T^F) -* E, 

T 2 - ( l D) is the result of appreciable mixing of the latter states through spin- 
orbit interaction with triplet states. In the case of ZnS, which shows the most 
mixing, the (primarily) 1 D levels 11* below Tj, - 3 A 2 (F), The energies of all 
three levels are close to one another and the energies and eigenfunctions given 
by the matrix diagonalization procedure are: 

9150 cm _l ,t(r = T 2 ) 8 ,36( l G) - .75( l D) - .16 ^jCPJ+.SS 3 A z <F) - .06 3 T 2 (F) - .09 3 T 1 
9604 cm“ l f t(r = E) = - .57 ( X G) + .77 ( X D) + .24 ^(P) + .14 3 T 2 (F) - .02 3 Tj (F) 

9805 cm " 1 ,t(r = T 2 ) = -^S^G) +.38( 1 D) + .12 3 T 1 (P)+.80 3 A 2 <F) + .24 3 T 2 (F) + .06 3 Tj(f: 
A plot of energy vs Dq for values of B, C, and X appropriate to ZnS:Nl is shown 
in Fig. 9. At the value Dq = 475 cm * the levels are appreciably mixed and 

3 

- A^(F) has become highest. The same plot for Dq = 400 1 the value for CdS:Ni, 
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3 

shows there is slightly less mixing and Tj - A 2 (F) is lowest. However, the 

ratio G/B is larger for the latter case and the effect is to push up the 1 D levels 

by about 1000 cm* - . The larger energy separation between the singlet and triplet 
states is evident in Fig. lc. The singlet levels are still relatively intense, 
broad absorptions. In spite of the fact that there is a greater energy separation 

between the and levels in CdS:Ni, there is about the same amount of triplet 

1 3 

character to the D states as there Is for ZnS:Ni. The coefficients of the F 

functions decrease, but those of ^T^P) increase to 0.41 for E and 0.32 for T^. 

Attempts to observe the A^ - ^T^(F) - T^, - ^T^(F) transition were generally 
unsuccessful, although the transition is symmetry allowed and was expected to be 
reasonably Intense, based on the p - d mixing model to be discussed later. No 
absorption bands having energies less than about 4000 cm" 1 on down to energies 
comparable to the strong lattice vibrational absorption of the host lattice were 
observed, with the exception of one weak, narrow band centered around 1060 cm" 1 
in a 3 mm thick crystal of CdS: ( ~ 1 percent Ni). The energy of the transition 
is approximately(9/2 )\ in fields of T fl symmetry, and a measurement of the 
transition energy could provide a good means to determine \. This would be true 
only if the energy difference between the levels were unaffected by low symmetry 
fields or if none were present. The calculations discussed above showed that 
the energy la drastically affected in some cases and the arguments given above 
to explain the breadth of certain bands, apply here as well. If the T 2 - ^(F) 
level is split and broadened for small values of the distortion coordinates, as 
the calculations indicate, this would help to explain the failure to observe the 
transition except for the one uncertain case noted. 

The local anisotropy about the metal atom position in the wurtzite structure 
has been discussed and it was pointed out that the anisotropy, measured by 8/d, 
is largest for ZnO and is nearly zero for ZnS. The spectra for the two 
polarizations^ J_ ^c and j| ^ are found to be nearly identical for ZnS-:Ni and CdS:Ni. 
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Nickel is expected to substitute for zinc without much distortion based on 
ionic size considerations only, and thus the spectra may be expected to partially 
reflect the trigonal (C 3v ) field component of the tetrahedra in the pure compounds. 
This would explain the observed isotropy of ZnS:Ni, In the case of CdS:Ni, the 
nickel ion substitutes for the larger cadmium ion and the lattice probably 
relaxes toward the smaller substituted ion. A totally symmetric shortening of 
the four distances is not unlikely since the forces in CdS tending to produce 
the C^ v symmetry are not large^based on the value 5/d = - 0.008. Of course, 
in both these cases there may be a local distortion, which if randomly oriented 
throughout the crystal, would lead to an apparent isotropy of the spectra. These 
distortions would be evident only by the observation of splitting of some of the 
levels degenerate in cubic symmetry, 

ZnO provides an example where the polar nature of the pure host crystal 
is maintained upon the substitution of either cobalt, nickel, or copper for zinc. 
The pure crystal anisotropy of this compound is the largest of the three discussed 
and the similar ionic radii of these three transition ions and zinc apparently 

helps to preserve the anisotropy. Splitting of some of the spin-orbit components 
3 ++■ 

of T^P) fdr Ni in ZnO is clearly shown by the polarized spectrum presented 
in Fig. 4a. The T^ and T^ components are expected to split into A 2 + E and 
A^ + E, respectively. From Table IV, we have the selection rules A^ - Aj^ 
only in the parallel spectrum and A^ — E allowed only in the perpendicular 
spectrum. Since the ground state is A^ the band at 15159 cm"*' for ^ | | is 

A^ and the 15200 cm" 1 band for J_ £ is E and the observed splitting is 41 cm" 1 . 

3 

The splitting energies of Tj- Tj^P) calculated by the use of Eq. (4) are 
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and since ( —) is negative for ZnO we see that should be at a lower energy 

than E, which agrees with observation. There is no contribution from the fourth 

power term, since the states are derived from an orbital P state, Watson's 

H-F wave functions^ for Ni^ were used to obtain < r^ > a 1.133 and CrS- = 3.021. 

Taking d = 3.7 a.u. we find ^ z 5 **5 ^TJ!^ • If we use the observed 

d . 41 d 

value of Dq in order to obtain (—r > and the H-F ratio above, the trigonal 
d -1 

field splitting is calculated to be 350 cm . On the other hand, using the H-F 
2 

value ^^ 2 ^ 2 0.083 directly, we calculate a splitting of 66 cm Analogous 
d 2 

Co Che case of ZnO.sCu previously discussed, we find the H-F value of > gives 

2 1 |. 

good agreement with the observed splitting, whereas the H-F ratio <X > 1 

d d 

together with the observed value of Dq gives too large a splitting. The observed 

splitting was used to calculate 5 0.053 and therefore < r^> = 0.725. The 

-1 d 

splitting of 41 cm observed for Zn0:Ni, which is comparable to that obaerved 
for ZnOsCu, seems to be a reasonable value for the trigonal distortion characteristic 
of the host lattice. 

The T^ - ^T^(P) level has a calculated splitting of the same magnitude as 
a 

T 2 " however the pattern is inverted. The identification of the levels 

given in Table V together with the polarized spectra shown in Fig. 4a leads to 
the assignment that the 15488 cm" 1- E J. cj^band is the E component and the 
15500 cm" 1 ^ 1| band is the Aj component of ^ - ^(P). The observed 

splitting is less than that for the X 2 level but the order of the levels is in 
agreement with prediction. The selection rule A, - A. only for E I I c is not 
obeyed. The intensity of A^ - A^ is less than that for A^ — E however, and the 
selection rules seem to be obeyed for all observed transitions only in the sense 
that the "forbidden" transitions are usually weaker and broader than the "allowed" 


ones 
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The large values for the f numbers and the quite good agreement between the 
observed and calculated spectra based on symmetry for reasonable values of 
the crystal field parameters is good evidence that Ni** occupies a regular metal 
atom lattice site in the three wurtzite-type compounds. Although the octahedral 
holes are large enough to accomodate a foreign ion, we have no evidence from the 
spectra that the sites are occupied by N1** . The nearby octahedron of metal 
atoms gives a positive energy contribution which 1 b greater than the negative 
energy contribution of the anion near neighbors. If some of the neighboring 
cation sites were unoccupied, the octahedral crystal field stabilization energy 
of Ni*"* could help to make the occupancy of these sites by Ni** energetically 
favorable. A small fraction of the Ni ++ ions may occupy the octahedral sites, 
but in concentrations too small to be detected in the Bpectra. 

The spectra of Ni in Cs 2 ZnCl^ and Cs^nBr^ shown in Fig. 3 provide two 
further examples of Ni + * in fields of T^ symmetry. The CSjZnX^ structure has 
no octahedral sites and a substituting foreign atom would have to replace a 
zinc atom in one of the discrete ZnX^~ anions. In these structures it is 
expected that the resulting NiX^~ tetrahedra could distort more readily than the 
tetrahedra in the wurtzite lattice since the ligand bonding to the rest of the 
lattice ia weaker in the former case. Therefore, any tendency of a tetrahedral 

complex containing Ni**” to undergo a distortion to a lower symmetry because 
k 5 

of the Jahn-Teller effect ’ is expected to be more readily observable in the 
Cs 2 ZnX^ crystals than in the wurtzite-type crystals. The observed spectra are 
very similar to that of Ni ++ in crystals having the wurtzite structure. The 
energy levels predicted based on T^ symmetry appear to fit the observed spectra 
equally as well as the cases formerly discussed. The splitting of the visible 

3 

T^(P) band is again not entirely accounted for by reasonable values of the spin- 
orbit parameters. The predicted levels shown in Fig. 3 were calculated for 





26 . 


X z - 250 cm * in both caaes. Low temperature spectra of the visible band 
reveals that the lower energy component of the band is further split whereas 
the higher energy component remains broad, in analogy with the former spectra 
discussed. The values of Dq obtained are smaller for the halide ligands than 
for the oxides and sulfides. The Dq values given in TaLle II are in the order 
Br~< Cl~ < S a 0 which is in agreement with the spectrochemical series 
based largely on results for ions in octahedral coordination. 

The good agreement between the observed spectra and that calculated for 
a (3d)^ configuration in T^ symmetry indicates that perturbations of the energy 
levels by states of the odd parity (3d) 7 (hg) configuration are negligible. 

The experimental evidence is thus in accord with the calculated estimates of 
g-d mixing previously discussed, Even though the mixing coefficients are small, 
they still may account for the major part of the observed Intensities of the 
transitions. This point is discussed under Intensities. 

The calculated energies depend on four parameters B, Cj X,,and Dq. The 
effect of X and Dq on the energy levels has been discussed. The parameters B 

and C primarily affect the energies of the ^P, *D, *G, and levels. In the 
3 3 

free-ion the energy difference between F and P is given by 15 B, and B is 
found to be 1020 cm * for Ni + * . To a first approximation, the energy of the 
^P level is 15B+6Dq 17,500 cm 1 . Examination of the spectra of Ni ++ in the 

crystals discussed shows that the energy corresponding to the ground state to 
^P transition is a function of the type of ligand and is in the order s' < Br’ <Cf"<0". 
This means that the value of B in the crystal is reduced from the free-ion value 
and the amount of reduction is dependent on the ligand type. Defining the ratio 
B (crystal)/ B (free-ion) as p , we empirically find a linear relationship 

between p and the ligand polarisibility which Is shown in Fig. 10. The values 
27 

of a for the ligands studied are indicated at the top of the figure and the 
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^ 4 " 

experimentally determined values of p for Ni are shown as crosses, The 
straight line fitted to the experimental data extrapolated to a = 0 giveB 
p : 0.93, not the free ion value p = 1, Nevertheless, the empirical relationship 
is useful to give the approximate value of p to be expected in a crystal. 

Jorgensen ^ 3 has recently discussed the factors which influence the value 
of the F Integrals in crystals. He has referred to the tendency to depress p 
below one as a function of the central ion and ligand type as the nephelauxetic 
(■ cloud-expanding) series. The reduction of p below one can be qualitatively 
described as a decrease of the interelectronic repulsion energy because the 
charge clouds of the ligands penetrate in toward the ligand and screen the 
electronic charges of the inner metal ion from one another. The greater the 
ligand polarizlbility, the greater the penetration and hence a larger reduction 
in p. The nephelauxetic effect can also be described as an expansion of the 
charge cloud of the central atom which means that the average interelectronic 
separation is greater. Jorgensen considers that the reduction in p is direct 
evidence of covalent bonding and he discusses the effect for a number of Mn "'"salts 
An attempt has been made to express the observed values of p in terms of 

29 30 

the "covalency" parameter e introduced by Koide and Pryce and used by Stout 
to account for the relative depression of the free-ion levels in MnF^. For an 
ion in tetrahedral coordination, the amplitude of the t 2 radial functions is 
reduced by the factor^ 1-e'. Calculation of the coefficients of the direct and 
exchange integrals yields the following electrostatic energies: 

3 T l (P) = (1 - 4.8e) (28 A q - 3pB fl - 21C fl ) 

3 A 2 (F) = (1 - 6 e) <28A o - 50B o - 21C q ) 

3 T 2 (F) = (1 - 56) (28A o - 50B o - 21C o ) 

3 T t (F) = (1 - k.2e) (28A q - 50B q - 21 C q ) 
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3 3 

The energy difference between P and the weighted mean of the F levels is 
(1 - 4.8e) (15B q ). Equating this to I 5 B, one finds (3 a (1 - 4.8e). 

Using the observed values of |3, one obtains values of e for the oxide, chloride, 

and sulfide of O.O 5 I, 0.062, and 0.094 respectively. According to the idea of 
29 

Koide and Fryce , this means coefficients of the ligand a-combinations of 
0.23, O. 25 , and 0.31 for these three cases. 

Spectra of Co ++ 

The spectra of Co ++ in ZnO, ZnS, and CdS have been interpreted on the 
basis of tetrahedral coordination about the cobalt ions. The energy level 
diagram for the d^ configuration given in Fig. 8 shows the effect of successive 
perturbations on the free ion levels. The crystal field splitting of the F 
state is inverted from that of Ni** , which makes the orbital singlet lowest. 

The orbital selection rules are A 2 - T^ for electric dipole transitions. The 
lowest energy transition ^Ag (F) - S ? 2 ^ should be weak, which has been found 
to be the case. The energy of this transition is about 3500 cm " 1 and is shown 
to be a weak broad band in Fig. 2. The two principal bands seen in the spectra 
are the transitions ^A 2 (F) - (F) and ^A 2 (F) - (P). The spin-orbit 

splitting of (F) is readily apparent in the 77°K spectra shown. The first 
order value for the splitting is given by 6 \, where X a - 178 cm " 1 for the free-ion. 

In spin-orbit coupling, the ground state is four-fold degenerate and 
transforms according to the double group representation G. The degeneracy is 
lifted in fields of lower symmetry. The effect of fields of less than cubic 
symmetry is, in general, to split the four-fold G levels into two doubly degenerate 
levels. In the low symmetry fields, all the states have two-fold degeneracy and 
are call.sd Kramers doublets. This degeneracy is only removed by interaction 
with magnetic fields. The zero-field splitting of the ground state has been 
found to be about 1 cm 1 for both ZnS:Co and CdS:Co. The field acting on the 
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ground state appears to have axial symmetry but gj_ and g|| are nearly equal. 

The average values g : 2.35 and g - 2.05 have been measured for ZnS:Co and 

* 

CdS:Co respectively . The selection rules given in Table IV are that transitions 
from G are allowed to all three of the double group representations; E^ E 5/2 

and G. The small zero-field splitting of the ground state implies that the two 
Kramers doublets are equally populated at all temperatures used and the ground 
state is accurately represented as having four-fold degeneracy. The selection 
rules thus Indicate allowed transitions to all the spin-orbit components of 
the upper state of an orbitally allowed transition, unlike the case of Ni ++ 
where transitions from the ground state are allowed to only one spin-orbit 
component of the upper state. 

The parameters were obtained by fitting the observed Bpectra to the 
calculated energy levels obtained by diagonalizing the matrix of , A set 
of solutions for a range of parameters C/B, Dq/B, and \/B were used to construct 
plots o£ E/B vs. the parameters. The plots were used to obtain good graphical 
estimates of the parameters which best fit the spectra. The final set of 
parameters obtained are given in Table II and the calculated energies of the 
excited levels relative to the ground state Bet equal to zero are ahown in 
Fig. 2. The data for Co ++ in spinel shown in Table II was obtained by D. S. McClure 
and ia included for comparison with ZnO:Co. The former spectrum is not shown, but 
is quite similar to the latter. The good agreement between the predicted and 
observed spectra indicates that the configuration is predominately (3d) 7 (negligible 
3d - kg mixing), and that cobalt is in a nearly regular tetrahedral environment 
of ligands. 

* The zero-field splitting and the g values were obtained from paramagnetic 
resonance spectra measured by R. E. Michel, RCA Princeton. 
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The polarized 4°K spectra of the visible bands are shown in Fig. 5. In 
complete analogy with the corresponding Bpectra of Ni + * , the spectra of Co ++ 
in ZnS and CdS are isotropic whereas that of ZnO:Co is quite anisotropic. The 
anisotropy of the latter case can be explained by a C^ v component of the crystal 
field potential. The precise form of the splitting, its magnitude, and the 
identification of the components has not been worked out because of the 
complexity of the spectra. The discussion which follows is facilitated with 
the aid of Fig. 11 whieh is a plot of E/B va. Dq/B for a value of X./B * - 0.30 
which is approximately correct for all three host crystals. Dq/B °*0.54 for 
both ZnO:Co and CdS:Co and Dq/B s 0.50 for ZnO:Co. 

k 

In the free-lon, all the doublet states lie at energies higher than P, 
however for values of Dq/B greater than about 0.3, three doublet states lie at 
lower energies than ^P. The doublets are G - V and G - ^ 

where we have omitted writing the parent free-ion states since the crystal field 
2 2 2 2 

states are mixtures of G, H, D and P, primarily the former two. These three 
states have very nearly the same crystal field configuration as the ground state, 

which has almost a pure eo nfl . B -iirati on , The three doublets lie quite close 

k 

to the components of P and are strongly mixed with the latter for ZnO:Co. 

k 

As Dq/B is increased, the P components and the three doublets rapidly diverge 
and do hot mix appreciably. For Dq/B = 0.54 the doublets have very little 
quartet character and are observed as weak transitions in the region 12,000 - 13,00 cm* 
The broken upper curve in Fig. 2b is the spectrum of a very concentrated sample 

of ZnS;Co which clearly shows some of the doublets on either side of the 
k 

prominent P band. In the case of ZnO:Co, there is considerably more mixing 
between the quartets and doublets and the sharp Intense lines on the low energy 
side of the visible band shown in Fig. 5a are taken to be the doublet states. 
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With this identification, all the spin-orbit components of V lie. in the region 
of broad band absorption in the energy range 16300 - 17900 cm . Ab a further 
check on this assignment, the spectrum was measured using a sample in which the 
concentration times the optical path length was about a factor of 10 higher 
than for the crystal used to measure the spectra shown in Fig. 5a. No absorption 
lines were observed in the region 9000 - 15000 cm \ however several additional 
weak absorptions due to the higher energy doublets were observed in the region 
18000 - 22000 cm" 1 . The assignment also gives values of the. parameter 0 for 
the oxide and sulfides which are close to the values obtained for Ni**in the 
same crystals. The values for Co ++ are shown bb circles on the graph in Fig. 10 
and indicate nearly the same linear dependence of (3 on a as that found for Ni 
Since the relative energy, E/B, of the transition to the lowest three doublets 
does not markedly depend on Dq/B, the fact that they lie at about 2000 cm 1 
higher energy in the oxide than in the sulfides reflects the larger value of 
the parameter B for the oxide compared with that for the sulfides. 

The anisotropy in the spectra of ZnO:Co, which is more pronounced in the 
region of the sharp transitions on the low energy side of the visible band, can 
be qualitatively explained by consideration of the selection rules for the 
group C 3v . In C 3v the ground state becomes E^ 2 + E^ 2 and » aa explained earlier, 
these levels are equally populated for the temperatures used. We see therefore 
from Table IV that transitions from the ground state are formally allowed to both 
doublets Ejy 2 and E^. The orbital selection rules are 

(e 4 t 2 3 ; A 2 ) ^ (e 3 t 2 4 ; A z ) for B || c 

< e V‘ *2> 

where the parent cubic crystal field configuration has been given preceeding 
the C 3v representation. Transitions from the quartet ground state to the 
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doublet excited states are made allowed by virtue of mixing via spin-orbit 

or other interactions between the doublets and other quartets. Since the 

3 b 

allowed transitions are to states with the configuration e t 2 , we will consider 
only how quartet states of different syrnnetriesjbut with this configuration^are 
mixed with the doublets and the configuration label will be dropped. The states 
will bf> labeled by two symbols, the first giving the representation in followed 
parenthetically by the T^ representation, e.g. V 2 '^) stands for the spin- 
doublet whose orbital function transforms according to E in which comes 
from an orbital Tj. state in T d when the symmetry is reduced from Tj to C^. 

Second order L*S^ interaction couples 2 E fa) with Sl and ^ (N^) 

whereas 2 A 2 ^T^ + 2 E ( 2 Tj) is mixed with ^ (^Tj), \ (Stj), 
h h k h 

( T 2 ), and E < Tg). The trigonal potential can cause further mixing of 
states of the same symmetry in C^. From this analysis, one thus expects 
transitions from the ground state ^Aj (^Aj) to the excited states 2 E ( 2 E) 
and 2 E ( 2 Tj) for El j. <: and to 2 A 2 C 2 ^) for E^ | |^. This is qualitatively 
confirmed by examination of Fig. 5a which shows that in the energy region corresponding 
to the doublet transitions, the E J. £ spectrum is more intense and has more 
components than the E | j spectrum. This does not explain the large number 
of components having amall energy differences which are observed. We think the 
large number of closely spaced lines must be connected with the coupling of 
vibrational and electronic motions. 

Within the framework of cubic symmetry only, the digital computer 

calculations have shown that there is appreciable mixing between ^T^ (P) component* 

2L 2 2 

and T!, and T^. These calculations also show that the E l/2 ‘ A^ state in 

-1 k 

ZnOiCo whose calculated transition energy la 17706 cm has about 0.2 P character. 

The absorption band observed on the high energy side of the visible band having 
an energy of 17679 cm ^ in the E J_ c spectrum is identified as the transition 
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G - \ 2 (F) - E ^ 2 - 2 A L . This transition is observed at about I 5 OOO cm " 1 

in the sulfides and calculations show that the upper level has about 0.3 character. 

The four principal absorption bands observed at 77°K in the region 
corresponding to the V 2 (F) - (F) transition have been interpreted as the 

four spin-orbit components of the upper state. The value of X needed to account 
for the total splitting is about - 200 cm 1 in all three cases. This rather 
large value also gives reasonably good agreement between the observed and 
calculated components. Comparison between the calculated and observed spectra 
is presented in Table VI for Co ++ in ZnO, ZnS, and CdS. The calculated transition 
energies, together with the designation of the upper level, are also shown on 
the spectra in Fig. 2. In general, the agreement is good; however there are 
some absorption bands which cannot be specifically identified. These are 
probably due to transitions to different vibrational states of the electronic 
levels. There are, however, only a few instances of what look like obvious 
vibrational frequency intervals following a 0-0 electronic transition. ZnO:Co 
at 4°K has a set of sharp absorption bands at 6003, 6129, 6177, and 6277 cm" 1 . 

The 6003 band has a half-width of 10 cm" 1 , and the next three are successively 
broader. There appears to be progressions corresponding to frequency intervals 
of 89 cm " 1 and 115 cm " 1 in the set of sharp absorption bands to the law energy 
side of the visible absorption in ZnO:Co. The lowest energy component of 
Sj. (F) in CdS:Co is an intense sharp band at 5230 cm ' 1 with a half-width of 
about 7 cm 1 at 4°K. This band is followed by two weaker and broader absorptions 
centered about 5300 and 5355 cm 1 . A set of three relatively narrow absorption 
bands is also observed at energies of 5470, 5520, and 5530 cm" 1 . In this crystal, 
certain patterns of sharp bands seem to be repeated at intervals of approximately 
290 cm " 1 which is close to the reported value of the longitudinal optical mode 
frequency = 300 cm 1 of pur** GdS^ 1 . The 4°K spectrum of ZnS:Co also has a 
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Bet of sharp absorption bands at low frequency intervals in the infrared. 

The energies are 5760*, 5889, 5917, 5940*, 5986, 6020, 6053, 6070. 6135, 6161*, 
6218, 6240*, and 6277, where the four sharpest and most intense bands are marked 
with an asterisk. 

The low frequency intervals may be due to an internal vibrational mode 
of the complex. However, in crystals with the wurtzite structure, the complex 
is an artificial construction and is really strongly coupled to the rest of 
the lattice. One then must be concerned not only with the nuclear motion of 
the transition metal ion and its nearest neighbors, but may have to consider 
the coupling of impurity modified lattice modes with the electronic processes. 
While we have not considered this problem, we wish to point out that the low 
frequency intervals may be due to the excitation of acoustic mode vibrations 
of the lattice. 

4 + 

It has been shown that the values of B obtained for Co show a reduction 
from the free-ion value of the same magnitude as that observed for the B 
parameters of Ni** in the same crystals. In the latter case, qualitative 
arguments were given in order to rationalize the observed reduction. The some 
reasoning leads to the expectation of a reduction of the spin-orbit parameter X 
when the transition metal ion is put in a crystal. The g values obtained from 
paramagnetic resonance spectra and magnetic susceptibility measurements of 
crystals containing tetrahedrally coordinated Co** can be explained on the basis 
of a reduced value of ^. Nevertheless, interpretation of the optical spectra 
of Co** we have measured indicates that X is even larger than the free-lon value. 
It may be that X has been overestimated since it has been obtained primarily 
on the basis of taking the entire splitting of the (F) band as being due to 
spin-orbit splitting. Although the value X s - 200 cm"*' gives agreement with the 
total splitting,, the calculated intervals do not agree with the observed Intervals. 
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If we take the two lowest energy absorption bands observed In the 77 K spectrum 
3 -1 

to be the E ^ 2 anc * G components of T^ (F), the values X. = - 135 cm and 
X s - 130 cm * are obtained for ZnS and CdS respectively. The two higher 
energy components calculated from these values of X are = 6680 cm 

G = 6780 cm * for ZnS and E^ 2 = 5990 cm 1 , G s 6095 cm * for CdS. Comparison 
with the observed spectra shows that the two upper components lie close to the 
third main absorption band, and this identification leaves unexplained the 
highest energy band of the set which lies some WO cm"* higher than the third 
band. This alternative explanation is not completely unreasonable, since we 
have not taken into account the effect of fields of lower symmetry on the cubic 
field levels. In addition to the expected splitting of the G levels, low symmetry 
fields may alter the energy separations between the components. Similar effects 
on the cubic levels of Ni ++ have been previously discussed. Calculation of the 
effect of low symmetry fields are being planned and the results are expected to 
be valuable for determining a better value of X. Although one may be able to 
account for the observed splitting of the T^ (F) band by a smaller value of X 
plus an additional low symmetry field, it should be noted that the larger value 
of X gives reasonably good agreement between calculated and observed transition 
energies. The observed energies of the Individual bands at low temperatures in 
the visible region for the two sulfides agree well with the calculated energies 
and only a few extra absorption bands are observed. 

Another direct determination of the spin-orbit parameter of Co ++ in crystals 
has been obtained from the infrared spectra of CoO, CoBr 2 , CoCi 2 and CoF^^. 

In these crystals, cobalt is octahedrally coordinated and the (F) level Is lowest. 
The most intense infrared absorption band observed was identified as the Ejy 2 - S^(F) 
— G - ^T^ (F) transition whose energy is given in first order by 6 X. Values of X 

in the range - 190 cm 1 to - 210 Cm 1 were obtained. These values are in agreement 

3 

with those indicated by the total splitting of T^ (F) for cobalt in the wurtzlte 
type crystals. However, the spectra of the octahedral cobalt salts are also 
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complicated by the effects of low symmetry fields. At the present stage of the 
interpretation of the spectra, one can only conclude that the effective spin- 
orbit parameter obtained for Co ++ in crystals is anomalously high. 

The relation between p and e can be obtained for Co ++ in a similar manner 
to that discussed for Ni ++ and the result is p = (1 - 4.2 e). The values of e 
for the oxide and sulfide are 0.052 and 0.093 respectively. These agree quite 
well with the values obtained for Ni ++ and indicate that the amount of ligand 
mixing with the central atom d-eleetron functions is approximately the same for 
these two transition metal ions. 

INTENSITIES OF THE TRANSITIONS 
33 

Ballhausen and Liehr have given an account of the intensities of crystal 
field absorption bands of tetrahedral complexes by considering the effect of 
central ion 3d - 4g mixing and also by using the molecular orbital (MO) 
approach. Although they were not concerned with the intensity of transitions 
between individual spin-orbit components, their methods may easily be extended to 
this case. The basic idea of the method is that some odd parity functions are 
mixed in with the even parity d-functions so that the first order functions are 

^1 8 ^i c jp I p where is the unperturbed d-function 

P 

and 4> ° stands for an odd parity function. The terms which contribute to the 
electric dipole moment of the transition ^ ^ are of the form 

< >c n + < t° W *”• 

P is the dipole moment operator^ e^r^. Considering central atom d-g mixing, 
one finds that only the t 2 orbitils of the d-electron are mixed with g orbitals 
by the third power potential term of Eq. (2) ? which is proportional to x y z . 

The coefficient is equal to 
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in this case. Dipole moment matrix elements for the 


transition e -» and t 2 — can be expressed in terms of a common parameter a), 
if it is assumed that the radial parts of the e and t 2 orbitals are the same 


and the energy denominators E° - E are the same for all i and p. The moments 


for the cubic field basis functions 



The consideration of o-bonding by the MD method gives only slightly different 
results for the relative intensities. The ^ 1 8 are linear combinations of ligand 
functions ana, just as before, only the t 2 central atom orbitals are mixed. 
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The model leads to the same set of basic moments given In Table VXX, except that 

/ 1-6 \ 1/2 

the t„ — t moments are reduced by a factor y = I- j- "nm J compared 

22 ' 1-2 S,/e-e / 

with the e - t 2 moments and the parameter yA appears Instead of o>. The factor 
1/2 

(1 - e) is the coefficient of the t 2 function, S is the overlap integral 
between the t ^ function and the ligand g^, combination, and A is the ligand 
coupling integral. The parameters have been introduced in this way, since we 

wish to estimate the effect of ligand mixing by the formalism of Kolde and 

29 1/2 

Pryce , and will take S to be zero, thus making y : (1 - e) . 

When one includes the effect of n-bonding^ further corrections to the 

intensities are obtained. In this case the e orbitals combine with ligand g fl 

combinations and the t^ orbitals are coupled with both g^and ligand combinations. 

The moments for n-bonding only are again the same as above, however multiplied 

by a third parameter B. Using the full MO scheme, the calculated intensity 

contains four parameters, three of which are the coefficients c^ p of o- and 

^combinations which are added to the e and t^ orbital functions, and one which 

gives the ratio of B to yA. With such a large number of parameters, the model 

is not quantitatively useful unless some Independent means is available to obtain 

some of the parameters. This was not the case and estimates of the correction to 

be expected because of n-bonding on the basis of physically reasonable parameters 

were found to be small. This aspect will therefore not be discussed and we 

shall only discuss the d - g mixing and o-bonding models. 

The actual intensity calculations were performed in a slightly different 

manner from that indicated above. It was more convenient to use the spherical 

harmonic representation for the d-functions, since all the symmetry functions 

employed were those which transform under the operations of the group like 
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the spherical harmonic basis functions . The perturbed d-orbltal functions 
for d~g mixing are 
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8 2 -- d 2 + I” “*o 

8 L - d 1 - ioip* 

8 o = d o 

8j = d'+itqpj^ 


where ai is the same parameter referred to above. The dipole moments in this 
basis to first order in cn are given in Table VIII, After expressing the basis 
functions for » |d n SU;nr >, in terms of the linear combination of Slater 
determinant n-electron product functions, they were corrected to first order by 
determining the linear combination of product d nl g functions to be added. The 
moments in the direction for this basis denoted by P (Ur, L' JT') are given 
in Appendix II for the d 8 and d 7 configurations. The factors <14-Re) which 
appear arise from the diminishing of the tg - t^ moments compared with the 
e - t 2 moments. In the first approximation, the mixing of central ion d and 
g functions Is considered and e is set equal to zero. 

The eigenfunctibns obtained from the computer solutions of the matrix of 'H' 
are of the form 


V(E tj r) Z Y J C SLJ I-” SLJ 5 rr >• The intensity of the transition 
Y(E o D - r* ) is given by 

iK(E t - t Q ) ^ C SLJ c 's*l’j' < d n sujrrl^l d n sl'j'; r'r'>\ 

r.Y"' 


where we have assumed that the degeneracy, including Kramers degeneracy, is 
effectively removed for the wave functions, and we sum over the squares of the 
moments. The moments in the expression are given in Appendix II. The calculated 
results for CdSjNi and ZnS:Ni are given in Table IX, where the intensities are 
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all relative to the intensity of the transition - 3 T^(F) - - 3 T^ (P), 

and the observed values were obtained from the area under the absorption peaks 
observed at 77 K. The calculated intensities of the two lowest energy transitions 
is much higher than observed and the comparison between calculated and observed 
intensities is only slightly improved by considering the effect of covalency 

through the use of the values of e previously obtained. The calculated value for 

1 3 

the transition to 1’ 2 - T 2 (D) is comparable to that for the transition to T 2 - Ag (F) 

in both cases, which shows that the singlet state is appreciably mixed with 

triplet states. It Is interesting to note that the triplet character is mostly 

\ 2 (F) for ZnS:Ni and the intensity is about one third of that for 3 A 2 (F); 

whereas the triplet character is 0.31 3 ^ (P)+0.25 3 A 2 (F) + 0.10 3 T 2 (F) for 

CdS:Ni and the intensity of the transition to l D has a calculated value which 

is even larger than that for the transition to 3 A 2 (F). It can be seen in 
1 3 

Fig. 1 that D lies much closer to P for CdSjNi than for ZnS:Ni. It has already 
been pointed out that in the latter case both E and T 2 of 1 D lie at lower energies 
than T 2 - 3 A 2 (F). 

The calculated intensity pattern for ZnS:Co is shown in Fig. 12. Comparison 
of this pattern with the observed spectra shows that whereas the total calculated 
Intensity of the multlplets is in reasonable agreement with that observed, the 
observed relative strengths of the individual components within a multiplet do 
not agree well with the calculated values. The intensity ratio 
l {\(F) - k T 1 (P)} 

1 (\<F) - \ (F)} 

obtained by taking the sun) of the calculated components is 4.2 compared to the 
observed ratio of 4.6. The calculated ratio is more simply obtained by use of 
the cubic field configurations appropriate to the weak crystal field for the 
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three levels involved . The calculated resultB are only slightly altered by 




1 + 1 . 


Caking e t 0 and we have only indicated whether the intensity of the transition 
is increased or decreased by the sense of the arrow on each line in Fig. 12. 

The correction is almost always in the direction which gives better agreement 
with observed spectrum, however for e <0.1, the correction never exceeds 5 percent 
of the intensity calculated for e = 0. 

The calculated intensity pattern for ZnO:Co is quite close to that shown 
in Fig. 12, except that the values for the transitions to G - ^E, Ejy 2 - ^ 

2 

and G - T^ are about twice as large as those shown. The observed relative 
intensities of the sharp bands in the spectrum which have been interpreted as 
being due to these transitions is much greater than the calculated values. Even 
though the model used seems to be only qualitatively correct, thiB is good 
indication that the doublet states are even more strongly mixed with quartet 
states than has been found by the diagonalization procedure for the parameters 
given in Table II. A serious limitation on the assessment of the validity of 
the intensity mechanism is that we have not tried to take into account the 
effect of fields of lower symmetry and the calculations are all based on strictly 
T^ symmetry. The calculations should be a good first approximation however, 

Blnce the cubic field predominates. Since the total relative intensities of the 
multiplets is correctly given by the model, at least for the Co"*"* spectra, it is 
possible that lower symmetry fields may cause sufficient mixing of the spin-orbit 
components within a multiplet to account for the differences between calculated 
and observed intensities. 

The discrepancies between the calculated and observed intensities for the 

++ 

case of Ni cannot be easily explained. The most serious disagreement is that 
the intensities of the transitions Aj^ - 3 T^ (F) - T 2 - 3 ^ (F) and A^ - 3 T^ (F) T 2 ~ 

3 T 2 (F). are calculated to be much too large. Characteristically, the observed 
absorption bands are more intense, the higher the energy of the transition. 
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Qualitatively, one might expect that the relatively low energy charge 
transfer transition, which has a high oscillator strength, could give the most 
important contribution to the intensities of the lower energy, predominately 
d-d transitions. The terms which would contribute to the intensity are in 
thl. c... <+J V|^ |P,U°> <» c M <> << l£j» c > 

< - \ 

where and i|f° are the zero order approximation to the initial and final 
states with energies E° and E^ respectively, and 4* c is the charge transfer 
state with energy E c - Therefore, the contribution from the second term is 
greater the higher the energy of the transition (the smaller the value of 
E° - E £ ). It is not clear just how important the first term would be, but 
since it only differs for different transitions by the factor < ^» c |.Pj > and 
since the energy denominator is larger than that for the second term, we expect 
its contribution to be leas than that of the second term. 
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In a paper concerned with the intensities of octahedral complexes, Englman 
recently presented a model in which the odd vibrations of the complex serve 
as the perturbing potential, and he discussed the relative contributions of 
the two terms above for those complexes. If we assume that the V 3 potential 
of a tetrahedral complex causes the charge transfer state to mix with the 
d-electron states, and if the first term above is neglected, Englmari's results 
may be used to calculate the ratio of the intensities of the transitions. The 
results for ZnS:Nl are, denoting the intensity of the transition by the upper 
level symbol, 3 T, (P): 3 A 2 (F): 3 T 2 (F) = 1: 0.31:0.16. This model gives better 
agreement between calculated and observed intensities for the low energy transitions, 
however the model is not precisely defined as we have used it for the tetrahedral 
complexes, and is dependent almost entirely on the energy denominators E° - E c » 
which of course leads to lower calculated values for the lower energy transitions. 





43. 

The results merely serve to indicate how the charge transfer mixing can affect 
the relative intensities. 

We shall probably have to take into account the actual form of the orbital 
angular momentum for the charge transfer states, together with ligand a- and 
n-bonding in order to obtain a quantitative model which is capable of predicting 
the correct relative intensities. An important requirement for the model will 
be to obtain the potential function which acts on the electrons anywhere 
within the whole complex, e.g. central ion plus the ligand environment, since 
the odd power potential terms are the ones which cause mixing between states of 
opposite parity. 
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Figure 1 


Figure 2 


Figure 3 

Figure 4 


Figure 5 


Figure Captions 

Absorption spectra of Ni + * in crystals having the wurtzite type 
structure. Spectra measured at 77°K and 298°K; (a) ZnO:Ni, 
axial spectrum; (b) ZnS:Ni; (c) CdS:Ni, (d) ZnOiNi, polarized 
with IE J_ c^ and JE | | ^c. The origin of the optical density 
scale is arbitrary and the curves have been shifted relative 
to one another in order to minimize overlap. The same crystal 
was used for the two different temperatures. 

Absorption spectra of Co in crystals having the wurt 2 lte type 
structure. Spectra measured at 77°K and 298°K; (a) ZnO:Co; (b) 

ZnS:Co; (c) CdS:Co. The origin of the optical density scale is 
arbitrary and the curves have been shifted relative to one 
another in order to minimize overlap. The same crystal was 
used for the two different temperatures. 

Room temperature absorption spectra of Cs 2 (Zn^ -y Ni ) X 4 where 
y « 0.1. The curve labeled A is for X - = C|" and B is for X* * Br". 
Polarized absorption measured at 4°K in the region of the visible 
band of Ni ++ in crystals having the wurtzite type structure; 

(a) ZnO:Ni; (b) ZnS:Ni; (c) CdS:Ni. The origin of the optical 
density scale is arbitrary and the absorption curves for the two 
polarizations have been shifted relative to each other in order 
to minimize overlap. 

Polarized absorption measured at 4°K in the region of the visible 
band of Co + *in crystals having the wurtzite type structure; (a) 


ZnO:Co; (b) ZnS:Co; (c) CdS:Co. The origin of the optical density 
scale is arbitrary and the curves for the two polarizations have been 
shifted relative to each other in order to minimize overlap. 
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Figure 6 

Figure 7 

Figure 8 

Figure 9 

Figure 10 

Figure 11 

Figure 12 


Absorption spectrum of ZnO:Cu; (a) axial spectrum measured 
at 298 3 K and 4°K; (b) polarised spectrum measured at 4 °K. 
Energy level diagram for states of the d® configuration in a 
crystal field of T d symmetry. 

Energy level diagram for states of the d^ configuration in a 
crystal field of T d symmetry. 

Energy level diagram for the states of the d® configuration 
in a crystal field of T d symmetry obtained by diagonalizing 
the matrix of for the parameters C/B = 4.5, B = 560 cm ” 1 
and X s - 250 cm -1 . 

Plot showing the variation of the normalized electrostatic 
repulsion parameter p = B/B q with ligand polarielbillty, a. 
Energy level diagram for the states of the d^ configuration in 
a crystal field of T d symmetry obtained by diagonalizing the 
matrix of $r for the parameters C/B » 4.5, and A/B ■ - 0.30. 
Calculated intensity pattern showing the relative intensities 
of the spin-orbit components of ZnS:Co. 






Table 1. Oscillator strengths of the prominent room temperature bands 
for Ni and Co in ZnO and ZnS. The weight percent of impurity is 
given parenthetically after the host crystal and d is the sample thickness 
in microns. 

Wavenumber of band _£ Upper level designation 

maxima, cm 


ZnO (0,072 Mi) , d ■ 560 


15260,16180 

1.6 x 10~ 3 



8340 

8.« x 1Q“ 5 


V'» 

4600 

3.9 x 10“ 5 


3 t 2 (p> 


ZnS (0.070 Nl) , 

d - 90 


12500,13060 

1.9 x 10“ 2 


3 t x (p> 

8550,9390 

2.1 x 10" 3 


\<n, x t 2 (d> 

4940 

1.6 x 10~ 3 


3 t 2 (f> 


ZnO (0.048 Co) , d > 460 


15310,16360 

5.2 x 10" 3 


\<9) 

6170,7150,7690 

7.8 x lO -4 


\<V) 


ZnS (0.19 Co) 

, d o 32 


1373.3,14930 

1.3 x 10’ 2 


\(P) 

6640 

2.9 r lO -3 


\<v> 






Table II, Parameters obtained from the spectra measured at 77 K 



dd) 

DqCcnf 1 ) 

B( cm” 1 ) 

C/B 

X( cm" 1 ) 

E g (cm” 1 : 

ZnO 

1.95 

420 

770 

4,8 

(-250) 

25 800 

Zns 

2.34 

473 

560 

4.5 

(-230) 

29 900 

CdS 

2.52 

400 

570 

4.7 

(-230) 

19 400 

Cs„ZnBr A 

2.40 

330 

695 

4.5 


44 000 

cs 2 *nci 4 

2.24 

385 

723 

4.5 


50 000 

Free Ion 



1020 

4.1 

-340 



Co 


*2 


d(8> 

IXrtcnf 1 ) 

B(em“ l > 

C/B 

X(ea*" 1 ; 

MgAl 2 0 4 (a) 1,92 

400 

735 

4.4 

-210 

ZnO 1.95 

390 

773 

4.5 

-210 

ZnS 2,34 

375 

610 

4.4 

-193 

CdS 2.52 

330 

610 

4.4 

-193 

Free Ion 


997 

4.2 

-178 


ZnO 3 00 


(a) The data Was obtained from an unpublished spectrum measured by 

D. S. McClure and is included for purposes of comparison with ZnO:Co. 






Table III. Character table for the group Tj showing how some basis functions 
transform and also showing the correlation between the representations in Tj and C 


Single 



Representations of some Correlated 

basis functions 

representation 

in C 3V 

A 1 

xyz 

A 2 

(3. 2 -l). Jl <*V) 

i E 

\,\S 

Aj 4* E 

*,y»* 

A 1 -1- E 

4i> 

E l/2 


. *1/2 

<H>. <H>. <H>> 

<B> E l/2+ E 3/2 


Double 




Table IV. Selection and polarization rules for electric dipole transitions, 
j There are no polarization rules in T d and a cross means allowed. The axis of 

I the group C coincides with the hexagonal axis of the wurtzite unit cell and | 

means allowed for E ||c and ]_ means allowed for E J_ c. 

I 

l 



I 


I 



++ 

Table V . Comparison of calculated and observed spectra of Ni ._jThe 
parameters given in Table II, with the exceptions that Dq ** 460 cm -1 for 
ZnS;Ni and Dq * 385 cm * 1 for CdS;Ni, were used for the calculation. The 
slightly smaller values were used in order to obtain better agreement 


with the sharp lines at the origin of the bands. The observed spectra 
were measured at 4°K. 


Energy 

in cto"* 

Energy in cm 1 

Energy in cm"* 

Designation oj 
upper level 

ZnO:Ni 

ZnSi 

t Ni 

CdSiNi 


Calc. 

Obs, <£_[_£> 

Calc. 

Obs. 

Calc. 

Obs. 


4226 

4213 

4390 

4215 a 

4001 

3902 b 

Tj - 3 T 2 CF) 

4359 

4248 

4450 

* 

4070 

- 

E - 3 T 2 (F) 

4330 

4312 

4500 

- 

4161 

- 

T 2 - 3 T 2 CF) 

4552 

4400 

4650 

4500 

4332 

4413 

a 2 - 3 T 2 (F) 

8580 

8340 

9810 

9750 

7991 

7840 

T 2 ~ \ <P) 

13091 

12404 

9180 

8760 

9736 

9710 

T 2 - *T 2 CD) 

13564 

13365 

9705 

9180 

9975 

10254 

E - *E CD) 

15192 

15200 

12448 

12436 

12198 

12199 

T 2 - \ CP) 

15262 

15300 

12493 

12404 

12284 

12294 

E - 3 Tj CP) 

15543 

13486 

12785 

12779 

12500 

12509 

Tj - 3 Tj CP) 

15729 

15635 

12930 

12933 

12712 

12623 

A i " * T i (i>) 

19655 


15230 

15320 

14995 

15100 

t 2 - *t 2 CG) 


a) A series of about eight sharp absorptions are observed in the energy 
range 4215 - 4300 cm " 1 

b) a series of about fifteen individual sharp absorption lines are obaervad 
in the anergy range 3902 - 4413 cm " 1 









Table VI. Comparison of calculated and observed spectra of Co . The calculated 
values are based on the parameters given in Table XI and the observed 
spectra were measured at 4°K . The energies are all in cm" 1 , 


ZnOtCo 


ZnS: 

sCo 

CdSi 

iCo 

Designation of 






upper level 

Calc. 

Obs. 

Calc. 

Obs. 

Calc. 

Obs. 


6140 

6001 

5922 

5920 

5245 

5220 

*1/2.1 

6687 

6643 

6431 

6221 

5775 

5604 


7527 

7156 

7016 

6829 

6366 

6135 

*3/2 

7474 

7691 

7158 

7208 

6519 

6496 

G 

15155 

15145* 

12296 

12377 

11974 

- • ' 

. G-V 

15425 

- 

12530 

12437 

12185 


*1/2 “ \ 




12488 




15844 

- 

13052 

13188 

12742 

13043 

G- \ 




13761 




16236 

16262* 

13860 

13881 

13324 

13541 

G 

16292 

- 

13910 

13974 

13530 

13618 

*5/2 4 t 1 <p:| 

16689 

- • . 

14057 

14118 

13829 

13847 

G 

16930 

- 

14455 

14485 

14126 

14142 

*1/2 J 

17906 

17679 

15056 

13006 

14574 

14485 

*i/a- a V 0 > C 

18187 

18094 

15197 

15343 

14664 

14760 

G - 2 T 2 (G> j 


a) A series of sharp absorption bands are observed to higher energies. 8ee| Fig. 5a. 

b) The four spin-orbit components lie in the broad band absorption'region. See Fig. 5a. 

e) The level has about 25 percent 0 - 4 Tj(P) character. 




Table VII.Dipole moment matrix elements for the cubic field basis functions 


"V l-l V *- Jb “ i- 

* \z_ y 2 1 1 ! v ’ ■ Jr “ -L 

* e t 2 Ul * ■- Jb “ J_. 

* V.J 3 Ul *** * "" 5~ “ i. 

-•,2 |i| ‘x, * m Jb“Js. 

- 1 «, Ul V - * Jr “ NpZ J- 
‘V ul v - ■ n/> 

"V Ul S. > ■ Jr •f-* i. 

* n4 * *1 Ul u ’ ■ * u Lin ' 





Table VXI1.Dipole moment matrix elements for the spherical harmonic basis 


functions. The moments are given for circular polarization, where L B — 

* J~2 

and R ■ (i+tj). The entries in the table are multiplied by i«. 








Table IX Comparison of calculated and observed intensities for ZnS:Ni and CdasNi. 

The observed intensities were obtained from the 77°K spectra. The moments, M, are 
in the k direction and are multiplied by itok, for the j>~d mixing model and by ilfA Jt, 
for the ©“-bonding model, and the parameter"*: is explained in the text. The intensities 

are all relative to that of the transition A - *T. (f) -»- T - 3x (p) which lias 

been set equal to one. 21 


ZnStNi 


M 

U. l,r4 ■ 0 

"calc.^- °- 1 

*obs. 

Designation of 
upper level 

0,013 (1-0.77 #) 

0.004 

0,003 

0.007 

i 

S 

0.228 (1+0.1 ft ) 

1 

1 

1 

t 2 - 5 Tj (P) 

0.064 (1+0.47* ) 

0.062 

0.067 

0.079 

i, - \ «> 

0.040 (1+0.5 ft ) 

0.023 

0.025 

0.040 

r, ■ % m 

0.351 (1-0.28ft ) 

0.83 

0,77 

0.052 

T 2 - 3 T 2 (P) 

0.834 (1-0.36 «) 

1,18 

1.08 

- 

t 2 - 5 t 1 (F) 

0.019 (1-0.53 C) 

1 , tor£ ■ 0 
calc. 

0.010 

CdSlNi 

I for€ ■ 0 

calc. 

0.009 

X obs. 

0.006 

Designation of 
upper level 
t 2 - \ (O) 

0.206 (1+0.15#) 

1 

1 

1 

t 2 - 3 t x (P) 

0.076 (1+0.13ft) 

0.11 

0.11 

0.084 

x 2 - X T 2 (D) 

0.046 (1+0.87#) 

0.040 

0.046 

0.10 

t 2 - 3 A 2 (P) 

0.319 (1-0.06# } 

0.79 

0.75 

0.055 

X 2 - 3 T 2 (f) 

0.829 (1-0.36#) 

1.39 

1.26 

- 

T 2 - 3 Tj (P) 
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Figure lb 
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Appendix I. Evaluation of the matrix elements of electrostatic, 

Bpin-orblt, and cubic crystal field potentials for the configurations d* and dj . 
The basis functions used for the calculations are those 

appropriate to the SLJM scheme. The relative phases are those specified 

70 

by Condor and Bhortley , and are determined according to Kq. Al. 

| SUM > • £ < SUlgl^ | JM > |.SLM s ^ > <A1) 

K s * \ a M 

where < SttlgMj | JM > Is the Clebsch-Gordon coefficient. The phases of the 
| SLMgM^ > functions ere fixed relative to the | SLSL > function for each 
SL term according to the method of Gray and Wills by the use of the lowering 

operators - iL and - 18 , The phase of each state with * 8 and 
* 7 20 

• L is the same as that given by Condon and Short)ay who did not, 
however, specify the 2 P state, This is specified by taking 

\%l%l > u (hsj /2 ^< 2 “ 1+ " 2 *> * 6C2V-2*) - 12(2 V-l") - 5 <aV-l + ) 

<2 + 0%l + ) *4^6 <*V-r> -9<l + r-l) -9 <lVo“) 
whera the terms in parentheses stand for the Slater determinant product 
functions with quantum numbers and jn (i). Since there are two 
0 terms In the configuration, another classification is needed, which we will 
call R3. Both Spin quartets have R5 ■ lit 2 H, 2 G, 2 P, b 2 D, 2 P have ■ 21j 
and bhi has tj * 10, where the terms b 2 D and a 2 0 are those given by Condon 
end Shortly. 
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The phases of the states of d 7 relative to those of d 3 are 

m 

fixed by correlating the three electron function rt <i) with the 

seven electron product function whose three missing pairs of quantum 
numbers are the negatives of the three pairs in d 3 , where the standard 

order of listing according to Condon and Shortley is used, and the phase 

7 *V 3 

of the d state is (-1) times the phase of the correlated d state. 

A„ Matrix elements of electrostatic potential, 

The matrix elements for the states of the d 3 configuration 
20 

are given below in terma of the Racah parameters. 

< 2 H |(?| 2 H >■ ~ 3A M + ZC 

< 2 G jQ| 2 G > * 3A - 11 B + 3C 

«s 2 F | Q | 2 F > * 3A + 9D + SC 

< 2 0(10)j Q | 2 D(10> > » 3A ♦ 7B + 7C 

< 2 D(21) | Q| 2 D(21) > = 3A 4 3B «■ 3C 

< 2 D(10)| Q | 2 D(21) > * 3 pr b 

< 2 P | Q | 2 P > ■ 3A - 6B ♦ 3C 

< 4 P | Q | 4 P > = 3A 

< 4 P | Q | 4 F > “ 3A 15 B 

The matrix elements for the states of id 7 are obtained by adding 18A - 28B ♦ 14C 
to the diagonal matrix elements of the states of d 3 above. Since we are only 
interested in energy differences, it is unnecessary to use different values 

4 

for the two configurations, and taking the F term to be zero, one obtains 
the electrostatic energy differences for the states of both d 3 and d 7 shown 
in Fig. 8. 
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B. Matrix elements of L • S 

The matrix of L • S ' 1 ^ was obtained by direct 

diagonalIsatlon. The matrix of L « S was calculated for ,he basis consisting 

of the nineteen Slater determinant one-electron product functions which 

have M ■ Z ( m (i) + m. (i) } * 1/2. There are also nineteen SLJ terms 
i«l 9 1 

and each term has one and only one state with M = 1/2. The unitary 
transformation S which gives the |SLJ 1/2 > functions in terms of the Slater 
determinant, functions was used to transform the matrix of L • S> in the 


m scheme to the SLJM scheme. 

The matrix elementa < A | L • I B > are given below for the 
states of the d 3 configuration in units of fthe one-electron spin- 
orbit parameter. The matrix elements for states of the d 7 configuration 
are the negatives of those for the same states of d 5 . The table is arranged 
in five columns, each containing ten digitB. Columns onf. end two give the 
symbols for A and B respectively in the form 

“♦iliil.il a I °°- 


The third column gives the matrix element in decimal form with the sign 
following the number. This number is also given as the square of the ratio 
of the prime factors 2, 3, 5, 7, and 11, Column four is the numerator 
end column five is the denominator. The first two digits of each column 
give the exponent of the prime factor 2, the next two digits give the 
exponent of the prime factor 3, and so on. For example, the sixth entry 


I 





In the table 


0221050900 0221040900 1.6248076- 0101000001 0000020000 

is the matrix element 

< d 3 , 2 H 9/2 Jjj • j^J d 3 , 2 G 9/2 > = - < d 7 , 2 U 9/2 JjL_ • d ? , 2 G 9/2 > 

» - 1.6248076^ * - | 2 1 3 ~, Al j /7 £ . 

C. Matrix elements of the cubic crystal field potential V^. 

The matrix of is diagonal only with respect to S in the 
SLJM scheme, however it is diagonal with respect to T , the irreducible 
representations of the groups C> h or T d and is independent of Jf t the 
degeneracy label* The functions may be classified according to SLJtTF 
and the correct linear combinations of |SLJM> functions which form a 
basis transforming under the operations of the groups 0^ or T^ like the 
spherical harmonics were obtained by the use of projection operators. The 
set of operators for each f » When applied to the individual |JM > 
gives a set of functions which form a basis for the representation 
*(D 

P[ | JM > ■ Z 00 • fi IJM > (A2) 

where OO is the complex conjugate of the ij th element of the 

transformation matrix of the irreducible representation P for the operation 
R. The expression for the operators applied to the spherical harmonics is 
given by 31 

p ij i* * ■ * < l lV <rt (R) . -w"-' M ‘ » <*» 




mItrix elements of l-s in units of the one 

jjTHE OR ^CONFIGURATION 


| R5 L2J 
J*11/2 

4 221051100 

9/2 

411030900 

411030900 

1 411030900 
221050900 
221050900 
m 221040900 

4?/2 

411030700 
.. 411030700 

j 411030700 
II 221040700 
221040700 
221030700 

JI5/2 

411030500 

B 411030500 
411030500 
411030500 

9 411030500 
411010500 
411010500 
411010500 

0 411010500 

221030500 
221030500 
221030500 
210020500 
210020500 
,4 221020500 

«! 13 / 2 

1 411030300 

411030300 

1 411030300 
411030300 
411030300 
,, 411010300 

411010300 
411010300 
411010300 

8 210020300 
210020300 
210020300 
■ 221020300 

221020300 
221010300 




fl/2 


411010100 

411010100 

221010100 


R5 L2J 

221O51100 

411030900 

221050900 

221040900 

221050900 

221040900 

221040900 

411030700 

221040700 

221030700 

221040700 

221030700; 

221030700' 

4110305001 

411010500 ' 

221030500 

210020500 | ( 

221020500 

411010500 

221030500 

210020500 

221020500 

221030500 

210020500 

221020500 

210020500 

221020500 

221020500 

411030300 

411010300 

210020300 

221020300 

221010300 

411010300 

210020300 

221020300 

221010300 

210020300 

221020300 

221010300 

221020300 

221010300 

221010300 

411010100 

221010100 

221010100 


.5000000 

1.5000000 

1*5811363 

.6000000- 

1.6248076- 

.6000000 


.7905692 

.6123725“ 

,7500000“ 

.9682458 

.2500000- 

1.1666666- 

.5270461- 

,9660918 

1.0540927“ 

.5000000 

1.5491934 

.3333333 

I. 5275249- 
.3333333- 
,5000000 
.7637624- 
.1666666- 

2.0000000- 

\ 1.4491376 
J 1.5811383- 

, .3333333- 

[ .6324557 

j 1.9720266- 
\ .7500000- 

» 1.1456440 
\ .9354143 
\ .2500000 
'll. 8371172“ 
■ . V333333 

r “ 

* . 8333333 - 

II. 2472190- 

( . 6666666 * 

\ 


electron parameter 


NUMERATOR 
2 3 5 711 


2000000 

10000 

2000000 

101000001 

2000000 


10000 

1000000 

2000000 

1010000 


200 

10000 

100000100 

100010000 


201000000 

100 


100 


200000000 

1000100 

10000 


100000000 

10100 

200000Q 

1000100 

100 

3000000 


20000 

100000100 

200000000 


ZETA FOR STATES 


DENOMINATOR 
2 3 5 711 

200000000 

200000000 

100000000 

20000 

20000 

20000 


300000000 

300000000 

400000000 

400000000 

400000000 

202000000 

102000000 

1010000 

2000000 

200000000 

10000 

2000000 

1000000 

2000000 

200000000 

201000000 

202000000 


100010000 

100000000 

2000000 

10000 

2000000 

400000000 

400000000 

300000000 

400000000 

300000000 

2000000 

202000000 

202000000 

2000000 


01 
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where oCpIf are the Buleri.au angles for the rotation R ( cl p t ). 

For the sake of completeness, we list the symmetry functions obtained 
in this way for the half-integral values of J up to 11/2, the maximum 
needed for the configurations d 3 or d 7 „ We only list the function for 
positive t , where t = ♦ 1/2 for E^ V ■ + 1/2 for % 5/f2 and if “ * 1/2, 

+ 3/2 for 0. The functions for negative if are related to those for positive 
IT as follows: change the sign of the M values and multiply the function 
by (~l)^ 2J " 1 ^ 2 for r B and change the sign of the M values 

and multiply the function by (~1)^ 2J+1 ^ 2 for T a G. The axis system 
is oriented so that x, £, and * are parallel to the two-fold axes of 
the groups. 

Representation 


J 


Symmetry function 

Y 

in 0^ or 


SC M( J «^i*J| J ' M > 



J « 1/2 


| 1/2, 1/2 > 

1/2 

®l/2 

J ■ 3/2 


1 3/2, 3/2 > 

3/2 ] 

t 

) 

r 0 



|3/2, 1/2 > 

1/2 J 

1 

J <= 5/2 

4- 15/2 

, - 3/2 > + 15/2, 5/2 > ] 

1/2 

B 5/2 



- 5/2 > * |5/2, 3/2 > J 

3/2 

? . 



- |5/2, 1/2 > 

V2 

1 

J « 7/2 

- 7/2 > +f7| 7/2, 1/2 > J 

1/2 

®l/2 


i J- |7/2, - 3/2 > + {7 | 7/2, 5/2 > J 

1/2 

B 5/2 


i [ |7/2, - 5/2 > 4 >1*3 | 7/2, 3/2 > j 

3/2 

) 


rM* 717 "' 



f 0 


- 7/2 > - {5 | 7/2, 1/2 > "j 

1/2 

1 




Al-7 


Symmetry function 


= 9/2 — ^9/2, - 7/2 > * ^14 1 9/2, 1/2 > + 3 1 9/2, 9/2 > J 

ft* V ,_ 

£{^l 


9/2 3/2 > 

9/2, - 7/2 > * <12 \ 9/2, 1/2 > - \f7 \ 9/2, 9/2 


v R 

1/2 

3/2 
1/2 


3/2 

1/2 


1 


1/2 

G~1 


L- ( . 5 | 9/2, - 7/2 > + /l4 |9/2, 1/2 > - 3 | 9/2,9/2 > ? 

'JaT (_ 1 

■ 11/2 ^z=|->l35 ^ 11/2, - 7/2 > +^6 111/2, 1/2 > - f? |ll/2, 9/2 > ^ 1/2 ^1/2 

fli j 11/2, - 11/2 > ♦ {22 | 11/2, - 3/2 > + fTl 111/2,5/2 >} 1/2 E 5/2 

i—— 11 I 11/2, - 5/2 > ♦ 13 f2 I 11/2, 3/2 > + \fl65 \ 11/2, 11/2 > j 3/2 

nJmU ' ■ 

I— f-JTlll/2, - 7/2 > ♦ ^70 I 11/2, 1/2 > + 3fl5 l11/2, 9/2 > J 1/2 

fiori > ’ * 

^— JjTT 1 11/2, - 5/2 > + 'TTi \ 11/2, 11/2 > j* 3/2 | 

i —{7 sJT \ 11/2, - 7/2 > >Ta2 I 11/2, 1/2 > - 4 \ 11/2, 9/2 > ^ 1/2 J 

f7B C 

These symmetry functions were used to calculate the matrix elements 
< SLJot; f t V 4 | S’L’J* «L' 1 r t ' »j' s* where V 4 given by Eq. (2) is equal to 

*($) v 4 ° * JTi/ »/?V 4 } ■ 

Since < SLJ *. 5 T t Ifj! ■ J C M <SLJM | <A4) 

and using Eq. (Al), the matrix element may be expanded to give 






A.1-0 


6 r i r t '. * kj sj ‘“'Vs 1 J *. V J * r i *j> <v (J ' *’ <Y »j’> * 

< S , L*M s , M l ' I J»M» > •< SLK^l V 4 | S'L'Mj^' > 

where oL denotes the dP^ in case there is more than one combination tor 
some . Using the Wigner-Eckart theorem 31 , we obtain 


4 L *\ 

! SLM lI V 4 I SL'M^ > * (-1) L 

-M. u ML 


(SL || V, || SL') 


where n ■ 0 or + U and 


I is Wigner'a 3-J coefficient which is related to the 


Clebsch-Gordon coefficient ahd (SL || II SL') is the double bar matrix 
element, which is independent of the values of M^ and M^'. The double bar 
matrix element was obtained by calculating the matrix element for one 
particular <SLM^ | | SL^' > for each SL term and applying the formula 

above. The matrix elements in the SLJct ;T scheme were then calculated 
by the use of an IBM 650 program. The double bar matrix elements (twenty 
non-zero values) were stored in memory and the 3-J coefficients were 
calculated as needed by using a sub-routine written by 8. R. Polo. 3 '*' 

The coefficients (Ja^ Vj) were fed in on punched cards for the values 
G, 3/2} B . 1/2, and E 1/2. 
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The tables are presented in six groups, which separate the matrix 
elements Cor each of the three irreducible representations for the spin 
doublets and the Bpin quartets. There are five columns, each containing 
ten digits. The first column gives the symbols for < SLJ<l| as follows, 

2 S * i| R 5 1 _ 2 J | _ ci * Th ® 8econd column denotes |SL' J* «c * > 

in the same format. The third column gives the matrix elements in fixed 
decimal form in units of (positive) Dq, with the sign following the number. 
The next two columns give the matrix element as the square of the prime 
factors 2 to 47. The first two digits of column four are the exponent of 
2, and so on up to 11. Sach digit in column five is the exponent of the 

successive prime numbers 13, 17, 19, 23, 29, 31, 37, 41, 43, and 47. 

To these exponents must be added (algebraically) the number with modulus 

55 given at the head of the group. For example, the first elemement under 

the heading spin doublets of G representation is 

< 2 H, J « 11/2, 1 | Vj 2 H, J « 11/2, *•!>*- 0.2144504 Dq 




4.23 

3*11.13 


Dq 


The matrix elements may be used for the states of d 5 and d 7 in 0 h or T d 
symmetry by taking Dq positive for d 5 in 0^ and d 7 in T d and by taking 
Dq negative for d 5 in T^ and d 7 in 0^. 


Wf^SSE SSfisWfe P'&BFiiSHi!**T.« 


N quartets of E (/ REPRESENTATION 


30040060 

R5 L2J 
411030901 
411030901 
411030901 
411030701 
411030701 
411010101 


R5 L2J 
411030901 
411030701 
411010101 
411030701 
411010101 
411010101 


4« 6666665 
2*9814238 
3.6514835 
* 6666663“ 
1.6329923 


5753555555 


200 

20001.0000 

101010000 

101000000 


SPIN QUARTETS OF E. REPRESENTATION 
30040070 


R5 L2J 
411030701 
411030701 
411030701 
411030501 
411030501 
411010501 


R5 L2J 
411030701 
411030501 
411010501 
411030501 
411010501 
411010501 


,8571417 
3.8332592“ 
2.3904568 
3.1428566 
3.2071345- 


SPIN QUARTETS OF G REPRESENTATION 
30040080 


R5 L2J 
411030901 
411030901 
411030901 
411030901 
411030901 
411030901 
411030901 
411030902 
411030902 
411030902 
411030902 
411030902 
411030902 
411030701 
411030701 
411030701 
411030701 
411030701 
411030501 
411030501 
411030501 
411030501 
1030301 
411030301 
411030301 
411010501 
411010501 
411010301 


R5 L2J 

411030901 

411030902 

411030701 

411030501 

411030301 

411010501 

411010301 

411030902 

411030701 

411030501 

411030301 

411010501 

411010301 

411^30701 

411030501 

411030301 

411010501 

411010301 

411030501 

411030301 

411010501 

411010301 

411030301 

411010501 

411010301 

411010501 

411010301 

411010301 


.5000000 
2.7277236 
4.6656947 
.9091369 
1.7320508- 
.3779643 
2.6457313- 
2.8333332- 
3.5634831- 
4.1661903 
.3779643- 
1.7320508 
.5773500- 
.0952365- 
2.9692298 
4.4446712 
1.0516401- 
2.4688536- 
1.5714283- 
3.8491981- 
1.6035672 
1.6035672- 

2 • 6186 145' 


5755555355 


2000000 

202010000 

100010100 

2 

102000100 


5353555355 


2000200 

1040100 

901020000 

106000000 

203000200 

202000100 

202000300 

200 

700020100 

107000100 

202000100 

203000200 

201000200 

400000000 

605000000 

582000002 

503000100 

901000100 

202000002 

303000002 

304000100 

304000100 

603000100 


NUMBERS 

1112233444 

3793917137 

5555555555 


5555555555 


555555555# 


200000000 



Crystal Field Matrix Elements (cant.) 


SPIN DOUBLETS 
30020060 


R5 L2J 
221051101 
221051101 
221051101 
221051101 
221051101 
221051101 
221050901 
221050901 
221050901 
221050901 
221050901 
221040901 
221040901 
221040901 
221040901 
221040701 
221040701 
221040701 
221030701 
221030701 
221010101 

SPIN OOU0LETS 

30020070 

R5 L2J 
2.21051101 
221.051101 
221051101 
221051101 
221051101 
221051101 
221040701 
22.1040701 
221040701 
221040701 
221040701 
221030701 
221030701 
221030701 
221030701 
221030501 
221030501 
221030501 
210020501 
210020501 
221020501 


OF ^ REPRESENTATION 


5551555453 5555555555 


R5 L2J 

221051101 

2.5454544 

221050901 

.6146576- 

221040901 

3.7451267- 

221040701 

4.1871789 

221030701 

1.8018748- 

221010101 

221050901 

2.5454544- 

221040901 

3.4465616 

221040701 

3.8533730- 

221030701 

.4264012- 

221010101 

5.0142653 

221040901 

1.5555555- 

221040701 

.4969034- 

221030701 

1.1547004- 

221010101 

1.8516401- 

221040701 

1,4444442- 

221030701 

1.2909942 

221010101 

2.0701966 

22103C701 

1.0000000 

221010101 

5.8797472- 

221010101 


404000300 

604010000 

307010001 

107020001 

104030001 

404000300 

303000301 

103010301 

104000101 

404000003 

200000302 

200010102 

203000102 

305000002 

102 2000000000 
3010102 
105010002 
4000102 
104000004 


OF E^REPRESENTATION 


R5 L2J 
221051101 
221040701 
221030701 
221030501 
210020501 
221020501 
221040701 
221030701 
221030501 
210020501 
221020501 
221030701 
221030501 
210020501 
221020501 
221030501 
210020501 
221020501 
210020501 
221020501 
221020501 


1.3333330- 

1.6035672- 

.6900652 

.7968186- 

7.3029674- 

1.5936380 

1.8571427 

1.6598500- 

1.9166296 

2.9277001 

3.8332593- 

1.2857142- 

.8247860- 

3.7796447 

1.6495720 

1.0476187- 

4.3643578- 

1.9047616- 

2.0000000- 

2.1821787 

1.8095237 


5553555355 


400000200 

104000100 

101010100 

300010100 

501010200 

500010100 

2000000 

5010000 

204010000 

203010100 

404010000 

6000000 

201020000 

202020100 

401020000 

200000002 

401020100 

600020000 

202000200 

201020100 

20Q0Q0Q00 


5555555595 


2000000000 


20000000 





SPIN DOUBLETS 
30020060 


R5 L2J 
221051101 
221051101 
221051101 
221.051101 
221051101 
221051 101 
221051101 
221051101 
221051101 
221051101 
221051101 
221051101 
221051101 
221051101 
221051102 
221051102 
221051102 
22105H02 
221051102 
221051102 
221051102 
221051102 
221051102 
221051102 
221051102 
221051102 
221051102 
221050901 
221050901 
221050901 
221050901 
221050901 
221050901 
221050901 
221050901 
221050901 
221050901 
221050901 
221050901 
221050902 
221050902 
221050902 
221050902 
221050902 
221050902 
221050902 
221050902 
221050902 
221050902 


OF 0 REPRESENTATION 


R5 L2J 
221051101 
221051102 
221050901 
221050902 
221040901 
221040902 
221040701 
221030701 
221030501 
210020501 
210020301 
221020501 
221020301 
221010301 
221051102 
221050901 
221050902 
221040901 
221040902 
221040701 
221030701 
221030501 
210020501 
210020301 
221020501 
221020301 
221010301 
221050901 
221050902 
221040901 
221040902 
221040701 
221030701 
221030501 
210020501 
210020301 
221020501 
221020301 
221010301 
221050902 
221040901 
221040902 
221040701 
221030701 
221030501 

210020501 

210020301 

221020501 

221020301 


.2144504- 
2.0344721- 
1.0893539- 
.7131493- 
6.6374644- 
4.3452408- 
.2322628 
.0999500- 
.2383946- 
2.1849238- 
2.3414775 
.4767892 
.5109520- 
3.7547173 
.3916082- 
.2214520 
.5637934 
1.3493254 
3.4352145 
2.5706791- 
1.1062441 
.4947272 
4*5342518 
1.5866575 
.9894544- 
♦3462369- 
2.5443125 
.2727269- 
1.4878491- 
.3692736 
2.0145572 
6.0302260- 
.6672848- 
.3232862 
.7977236- 
5.8620498- 
.1740775 
1.2792040 
2.0889318 
1.5454543 
2.0145572 
2.0925551- 
4.6056618 
a 5096469 
1.4814874 
3 . 6556307 " 
1.2792040- 
.. I t 77276 
.2791451 


5351555353 


602000200 
7040'10200 
703020300 
704020200 
406020301 
407020201 
307000101 
304010101 
1102000101 
1303000201 
604000401 
1302000101 
603000301 
706000301 
804000400 
407010100 
404030200 
110010101 
107030201 
207010301 
204020301 
204010301 
405010401 
506010201 
404010301 
505010101 
600010101 
206000200 
205040100 
105000201 
104040101 
604020201 
605010001 
203000001 
204000301 
307000301 
203000201 
306000201 
605000201 
204000200 
104040101 
103000201 
403020301 
404010101 
204000101 
205000401 
306000201 
204000301 
305000101 


3555555555 


2000000 

2000000 

1000000000 

lOOQOOOOOO 

IOOQOOQOQO 

1000000000 

1000000000 

IOOOOGOOQO 

1QQ0Q000Q0 

1000000000 

1000000000 

1000000000 

1000000000 

1000000000 

1000000000 

1000000000 

1000000000 

lOOOC'OOOOO 

1000000000 

1000000000 

1000000000 

1000000000 

1000000000 

1000000000 

1000000000 

1000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

4000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2200000000 

2000000000 

2200000000 

2000000000 

2000000000 

4000000000 

2000000000 

2000000000 

2000000000 

2000000000 



Crystal Field Matrix Elements (c.ont.) 


221050902 

221040901 

221040901 

221040901 

221040901 

221040901 

221040901 

221040901 

221040901 

221040901 

221040901 

221040902 

221040902 

221040902 

221040902 

221040902 

221040902 

221040902 

221040902 

221040902 

221040701 

221040701 

221040701 

221040701 

221040701 

221040701 

221040701 

221040701 

221030701 

221030701 

221030701 

221030701 

221030701 

221030701 

221030701 

221030501 

221030501 

221030501 

221030501 

221030501 

221030501 

210020501 

210020501 

210020501 

210020001 

210020501 

210020301 

210020301 

210020301 

210020301 

221020501 

221020501 

221020501 

221020301 

221020301 

221010301 


Spin doublets of G representation (cont,) 


221010301 
221040901 
221040902 
221040701 
221030701 
221030501 
210020501 
210020301 
221020501 
221020301 
221010301 
221040902 
221040701 
221030701 
221030501 
210020501 
210020301 
221020501 
221020301 
221010301 
221040701 
221030701 
221030501 
210020501 
210020301 
221020501 
221020301 
221010301 
221030701 
221030501 
210020501 
210020301 
221020501 
221020301 
221010301 
221030501 
210020501 
210020301 
221020501 
221020301 
221010301 
210020501 
210020301 
221020501 
221020301 
221010301 
210020301 
221020501 
221020301 
221010301 
221020501 
221020301 
221010301 
221020301 
221010301 
221010301 


.4558421 
.1666643- 
.9092409- 
.7776156- 
1.8070157- 
.3030445- 
.9258196 
2.6457513- 
1.2121827- 
3.4641016 
2.1213203 
.9444443 
.5939133 
1.3801308 
1.3887300- 
4.2426407 
.5773500- 
5.5549201- 
.7559285 
.4629093 
.2063468- 
.1844261 
1.4846147- 
2.2677868- 
3*3946736 
2.9692298 
4.4446712- 
1.7320508 
.1428565 
.6388763 
2.9277001- 
3.5856857- 
1.2777531- 
1.5649214- 
3.5138210 
.5238091 
2.1821787 

2.6726124 
.9523807 
1.1664236 
4.7142856 
1.0000000 
2.4494897- 
1*0910893- 

2.6726124 


2.6726124 


.9047618- 

2.2162048 


604000101 

2000202 

1040102 

701020002 

704010002 

106000002 

305000102 

204000302 

506000002 

405000202 

106000202 

202 

500020102 

503010102 

107000102 

306000202 

203000202 

507000102 

404000102 

105000102 

200000002 

203010002 

407000002 

406000102 

303000104 

607000002 

504000004 

205000202 

204000002 

404010002 

405010102 

306010102 

604010002 

505010002 

204010004 

202000004 

403020102 

304020102 

602020002 

503020002 

206000004 

204000202 

305000202 

203020102 

304020102 


304020102 


202000002 

303000002 


2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2200000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

4000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 

2000000000 


2000000000 


2020000000 

2020000000 
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Appendix II. Electric dipole moments for the cubic field basis functions 
B 7 

of the configurations d and d . 


A. The d configuration 

The ground state of Ni * in T^ symmetry has the symmetry Aj and is 
primarily l^F^; A^ > with a small amount (0.1 - 0.15) of 1 3]E> 0 • >• 

The eight terms required to calculate the intensity of transitions are 
given below for the k direction. The common parameter Q equals w and 
6=0 for the £ - d mixing model, and is equal to if A with 
t ^ 0 for the (T-bonding model. 


< 'V MU T 2 » o > . . 1 | (i.) , (1 .6) k 

'■ V MU \> h *0 * * - 1 it o « - w*> i, 

, , . . 1/2 

< 3 f 4 ; A 1 I P| 3 f 3 ; T 2 , o > - . I §-<±) Q (1 - f e> k 

* \> \ 1 1 1 \ T J v ■ 1 r ! k >1/2 « (1 - fe‘>£ 


1/2 


< V A 1 h*o>= ( fo> Q (1 - 3« > k 

T * ft 1 1/2 

< V A 1 I'I V T 2 J o > “ ' 1 I ( I05 ) 9 (1 - 3€ > * 

< k x ip I 3 f 3 * t 2 r c > ■ o 

<V. Al iPi\,T al# >. -iia oL ) 172 


In a weak 
levels have nearly 


crystal field of T d symmetry, the four spin triplet 
the following configurations. 
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Vn-c!) VS {t J V-2t J V> 

% (F) = t/. 2 

3 T 2 (F) t 2 5 e 3 

\ *t,V) 

The moments for the spin-orbit components in this basis are 
3 .3 2 l 1/2 

< A x - (F> | P | T 2 - T x CP) > = - i | ) Q Cl * 6) k 

< A 1 " 3t 1 <F) I 11 T 2 “ % < F > > “ - * fj 0 £ 

< A x - \ CfO I P I T 2 - 3 T 2 <I**> > * i Cl - Jq6)Q £ 

* *i - ’ T i «» If I T j - \ <f> *■ - t f < 35 >* 2 « - 3*> « £ 

3 3 L 1 1/2 3 

< Aj - 3 t i (P) 1 1 1 t 2 - 3 Tj (P) > = if c ~ ) g Cl - ft) k 

< A x - \ (P) | P I T 2 - 3 A 2 (F) > =■ - i —• Q k 

< *1 - \ < r > l£l t 2 - % (F >»* -ih 0 « (1 - f 4> js 

< Aj - \ CP) I PIT, - \ <10 > - - ifc i- ) Q <1 ♦*) fc 

B. The d 7 configuration 

The moments for the spin quartets are given both in the |SLJ*t.flf> 
and |t 2 m e 11 <^r *)tPJT s> bases. The transformation matrices connecting 
these two bases are given below. 
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G 


t 2 e< 4 T 2 )-l 

t 2 e<4 ' r 2 ) " 2 

t 2 e(\)-l 

t 2 e( 4 Tl )-2 

t 2 e 2 ( 4 T 1 )-l 

t 2 .e 2 ( 4 T 1 )-2 

4 P -1 

3 

9__ 

13 

_7_ 

_ _1_ 

_7__ 

_ 1 

P 9/2 1 

2 iiT 

4 435 

4T3T 

4 j3? 

4^35 

2 or 

2 4 35“ 

A P -2 

i 

_9_ 

7 

1 

_ 3___ 

1 

- 3 

F 9/2 2 

2 

4 4TT 

4 4iT 

4HT 

4 or 

2OT 

2 05* 

4_ . 

2 


_ 5 

4 

3 

8 

6 

F 7/2 

"nTTT 


| 2 1Q. 

>| 210 “ 

Pair 

OTo” 

{ 210 ~ 


rn 

7 

_ _9 _ 

_3_ 

. 3 

3 

3 

F 5/2 


2'TtT 

2^70“ 

2 or 

2O0" 

470- 

4 70 

4_ 

1 

6 

3 


3 


6 

F 3/2 

F 

Jl05“ 

vTioT’ 


\[i05" 


4105 

4 P 





2 


1 

P 5/2 





or 


“ 4T~ 

A p 




2 


1 


P 3/2 




\"r 


“pr 



Only the moments In the k direction are given and we have 


I C G 
1 < G 
I < G 


-* * 1/2 5 
— * 3/2 5 
—► G) 


- (G.3/2 -* E 1/2( - 1/2)J 2 

■ [2M (G,-l/2 -*• E 5/2f - 1/2 )J 2 

= j2M (G,-l/2 -Hr G,3/2)J 2 + j^M <G,3/2 


The momenta in the | t 2 m e n (^P‘)5p^> basis are given below 
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Finally, we list the moments in the |SLJ<t ;TIf > basis. 

The same rules for calculation of the intensity apply. There are four 


columns of ten digits each. The first column gives < SLJot ;T*1 in the format 

_ i I ill.il 

where 

r t B 6 stands for ^/2 and ^ = 1,2 stands for 1/2, - 1/2 

pj =7 stands for and B 1,2 stands for 1/2, - 1/2 

« 8 stands for G and Vi = 1,2,3,4 stands for 3/2,1/2,-1/2, 

-3/2 

The second column gives |SL*J* ot' | p* y‘ > in the same format . The third 
and fourth columns give the moments M and M* respectively for the total 
moment expressed in the form 

< 3/2 L J <*. ; r r | P | 3/2 L* J* <*.’ ; f' y* > » (M + | M*) Q where 
Q has been defined previously. 
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d 7 CONFIGURATION DIPOLE MOMENTS 


G - E 1/2 


L2J a P r 

L2J a p x 

M 

M 1 

309010801 

309010602 

.3320079 

. 1593638- 

309010801 

30701C602 

.1128436 

.1247219- 

309010801 

101010602 

.2618614“ 

.1745743 

309020801 

309010602 

.4260064- 

.24343cl 

309020801 

307010602 

.3538151- 

.1905158 

309020801 

101010602 

.4000000 

*2666666- 

307010801 

309010602 

.3903600 

.3253000- 

307010801 

307010602 

.0872871 

.0581914- 

307010801 

101010602 

.1069044 

.0712696- 

305010801 

309010602 

.0488311- 

.0450748 

305010801 

307010602 

.4753953- 

.3678854 

305010801 

101010602 

.0658360 

*1975082- 

303010801 

309010602 

.0828078 

.0552052- 

303010801 

307010602 

.2839181 

.1728197- 

303010801 

101010602 

.1209486 

.3628458- 

105010801 

309010602 

.0983819- 

.2529822 

10501.0801 

307010602 

.3331258 

.3394112- 

10501.0801 

101010602 

.0923760- 

.2771281 

103010801 

309010602 

.1264911 

.0843274- 

103010801 

307010602 

.2262741- 

.3016988 

103010801 

101010602 

-.0461880- 

.1385640 


G - E 5/2 


L2J a P r 

L2J a P r 

M 

M ' 

309010803 

307010702 

.1999622 

.1825741- 

309010803 

305010702 

.4937860 

.3265986- 

309010803 

105010702 

.1891221 

.0872871- 

309020803 

307010702 

.0664015 

.0398409- 

309020803 

305010702 

.0772088 

.0712696- 

309020803 

105010702 

.1555555 

.4000000- 

307010803 

307010702 

.0141972- 

.0851835 

307010803 

305010702 

.3301587 

.2476190- 

307010803 

105010702 

.2256872- 

.2138089 

305010803 

307010702 

.0811482 

•0221312- 

305010803 

305010702 

.1088717- 

.1451623 

305010803 

105010702 

.0370327 

.0740656 

303010803 

307010702 

.4156136- 

.2529822 

303010803 

305010702 

.5576042 

.3394112- 

303010803 

105010702 

.3830039 

.2418972- 

105010803 

307010702 

.0276026- 

.0552052- 

105010803 

305010702 

.0370327 

.0740656 

105010803 

105010702 

.0923760 

.2771281- 

103010803 

307010702 

.3312314 

.4416419- 

103010803 

305010702 

.2756886 

.2715738- 

103010803 

105010702 

.0461880 

.1385640- 
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d 7 CONFIGURATION DIPOLE MOMENTS 
G ~ G 


L2J a P r 

L2J a P r 

M 

M ' 

309010803 

309010801 

.3381175 

.2397561- 

309010803 

309020801 

.3662335 

.3098898- 

309010803 

307010801 

« 1656431- 

.0752923 

309010803 

305010801 

.1495369 

.1199773- 

309010803 

303010801 

.0766651- 

.0511101 

309010803 

105010801 

.1886740 

.4293960- 

309010803 

103010801 

.1171080- 

.0780720 

309020803 

309010801 

.0469530- 

.1408590 

309020803 

309020801 

.2151657 

.2151657- 

309020803 

307010801 

.0690065 

.0230021- 

309020803 

305010801 

.3399761 

.2310775- 

309020803 

303010801 

.2342160 

.1561440- 

309020803 

105010801 

.1689473 

.1788854- 

309020803 

103010801 

.3577708 

.2385139- 

307010803 

309010801 

.6876699 

.4366954- 

307010803 

309020801 

.2990284 

.1610152- 

307010803 

307010801 


.0491807 

307010803 

305010801 

.1888235 

.1320344- 

307010803 

303010801 

.5529389- 

.3964468 

307010803 

105010801 

.1221789- 

.0956182 

307010803 

103010801 

.1912365 

.1274910 

305010803 

309010801 

.1043281 

.0782460- 

305010803 

309020801 

.2709184 

.1673320- 

305010803 

307010801 

.3875851- 

.1661079 

305010803 

305010801 

.3408225 

.2109853- 

305010803 

303010801 

.4445258- 

.3397189 

305010803 

105010801 

.3643546- 

.2649851 

305010803 

103010801 

.1269720- 

.0331231- 

303010803 

309010801 

.1916629 

.1277753- 

303010803 

309020801 

.0585540 

.0390339- 

303010803 

307010801 

.3442827 

.2295218- 

303010803 

305010801 

.0433683 

.1301051- 

303010803 

303010801 

.0354101 

.1062303- 

303010803 

105010801 

.0180299 

.0540898- 

303010803 

103010B01 

.0540898 

.1622696- 

105010803 

309010801 

.0975900 

.1951800- 

105010803 

309020801 

.0298142 

.1788854 

105010803 

' 307010801 

.1753001 

.0956182 

105010803 

305010801 

♦3643546- 

.2649851 

105010803 

303010801 

.2974942- 

.2163594 

105010803 

105010801 

.0206558 

.0619677- 

105010803 

103010801 

.0619677 

.1859032- 

103010803 

309010801 

.2927700 

.1951800- 

103010803 

309020801 

.0894427 

.0596284- 

103010803 

307010801 

.1912365- 

.1274910 

103010803 

305010801 

.1177711- 

.3533135 

103010803 

303010801 

.0540898 

.1622696- 

103010803 

105010801 

. 1445913- 

.4337741 

103010803 

103010801 

.0826236 

.2478709- 
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The Optical Absorption Spectrum of WiSiF^"611^0 
Donaid S. McClure 

R.C.A. Laboratories, Princeton, N. J. 

ABSTRACT 

A detailed study of the visible and near infrared absorption 
bands of NiSiF g *6H 2 0 has revealed many features commonly found in 
polyatomic molecular spectra. The most important result is the 
discovery that the excited states are not usually stable in the 
nearly cubic symmetry characterizing the ground state. This leads to 
the breadth of the absorption bands. The true electronic origins of 
many of these bards have been located approximately. The multiplet 
splittings which appear too large for the known value of the spin- 
orbit coupling energy, are understood when it is realized that some 
multiplet components are stable in the cubic configuration and others 


are not. 



I. INTRODUCTION 


The optical absorption of NiSiF g .6H 2 0 arises from the d 8 shell of 
the Ni +,f ion. There is one such ion per unit cell of the rhombohedral 
lattice. Each Ni ++ is surrounded by six H 2 0 lying at the corners of an 
octahedron which is compressed slightly along a 111 cube axis. The 
Sil' 6 groups and the Ni(H 2 0) g groups are packed as in the CsCl structure, 
but with a 111 trigonal distortion superimposed. The relevant structural 
data are shown in Table I. 

because of this simple structure the NiSiFg‘6H 2 0 crystal is 
especially suitable for a spectral study of the absorption spectrum of 
Ni ++ . The crystal field is almost cubic, but the small trigonal field 
provides a useful perturbation to help in identifying spectral features. 
The orientation of the trigonal axis at the Nl ++ ions coincides with the 
trigonal axis of the crystal. The substitution of DjO for provides 
another useful variable. 

The principal absorption bands in the spectrum have already been 
identified. Only a few of the singlet levels have not been observed. 

The purpose of this paper is to report on the fine structure observed 
throughout the spectrum and to give an interpretation of the band shapes. 
ThiB leads to a consideration of the effects of Jahn-Teller distortions 
and other instabilities in excited levels. It is shown in detail how 
the spectra are modified by the change in shape of the complex 
upon excitation. 




II. GENERAL DESCRIPTION OF SPECTRA 


a. Energy Levels in the Crystal Field 

The energy levels of the d R configuration in a cubic crystal 
field have been given in great detail by Liehr and Ballhausen^. The 
ground state in an octahedra 1 field is The calculated positions 

of the other levels are shown in Fig. 1, using parameters which best fit 
the NiSiF 6 *6H 2 0 spectrum at low temperatures. The calculated and 
experimental level positions are in good agreement, as Fig. 1 shows. 

The parameters used are: Dq = 910 cm -1 , X = 300, B = 8B0, 

C = 3960 (or F 2 ~ 1446, = 113.2), using the standard notation^. The 

positions and identification of individual components of a multiplet are 
not obvious, but the calculated triplet widths are in approximate 
agreement with observed band widths. Only in the broad band at 15,000 cm 1 , 
the 3 Tj (F) state is gross structure clearly discernable. The width of this 
structure is appreciably greater than the calculated width. 

b. Band Structure 

Some of the most interesting features of the spectrum are the 
regions of fine structure near the origins of the stronger bands. This 
structure is observable at 4.2°K in crystals 2 - 4 ran thick. It has 
been observed at the origins of 3 Tj (F), 1 T 2 (D) and ^ (p). Pappalardo^ 
has already reported the latter. Figures 2, 3 and 4 show the fine 
structure of these three regions. The spectra were taken with a Cary 14 
spectrophotometer and with a three-meter grating spectrograph. Tracings 
of the photographic plates showed no more structure than did the Cary 
records, but the resolution was somewhat better and the wavelength 
measurements were more accurate. These measurements are given on the figures. 



Roth normal and deuterated salts showed substantially the same fine 
structure. The widths of the narrowest lines measured from the plates 
were 4 cm” 1 for the 3 T^(P) and 3 T^(F) bands, and 1-2 cm -1 for the 
1 T 2 (D) band. The fine structure is partially polarized} in fact it is 
more strongly polarized than the main spectral bands. It is observed in 
all cases at least halfway up the side of the band after which it becomes 
less well resolved and merges with the coarser structure. 

The 5 T 2 (F) band origin is overlaid with overtone absorption 
of the H 2 0 molecules; therefore DjO was substituted. Some fine 
structure was observed, but it could all be attributed to H 2 p 
impurity or D^O overtones. However, the structure of the main band was 
brought out more clearly by deuteratlon. This structure is seen in 
Fig. 1. 

The 3 Tj(F) band has the complex and interesting structure at 
4.2°K shown in Fig. 5. The two broad peaks on the long wavelength side 
are more intense in the polarization perpendicular to the c-axis. 

The short wavelength peak at 15970 cm -1 is the sharpest of the strong 
features and appears to have the same frequency and intensity in both 
polarizations. The companion peak at 15780 cm” 1 on the other hand shows 
structure which is different in the two polarizations. This is the only 
intense feature in the entire spectrum to show such behavior. Associated 
with these two sharp peaks are progressions of what appears to be a totally 
symmetric Ni-0 frequency. An O-H frequency is also distinguishable, 
riper^ established the identity of the 0-H bands by noting their 
shift on deuteratlon. These vibration bands are reminiscent of those of 
molecular spectra. They appear only in this one spectral region. 



4. 

The 1 T 2 (D) band is the strongest of the intercombination bands 
so it was possible to resol\e the fine structure in a crystal 4 mm thick. 

In fact more lines were found here than in the other two bands. 

The (P) band shows no evidence of resolved multiplet. 
eomponent.8 and is featureless except for the fine structure at the 
origin. 

III. INTERPRETATION OF BAND STRUCTURE 

a. Sources of Spectral Structure 

The structure of these absorption bands is expected to be 
complicated for three reasons. First is the fact that they are forbidden 
and may only appear by way of odd vibrational perturbations which 
destroy the center of symmetry. Second, almost all of the excited levels 
are degenerate and may suffer Jahn-Teller distortions. Even states which 
are hot degenerate may become stabilized in unsymmetrical positions 
through the effects of coupling to nearby states. Third, the equilibrium 
configuration of the atoms in the excited state may differ by a totally 
symmetric displacement from that of the ground state. 

Due to the action of the Franck-Condon Principle, progressions 
of symmetric vibrations would be excited. Their intensity maxima would 
occur at various distances from the origin depending on the amount of 
displacement. In the case of the unsymmetrical upper state, progressions 
of non-totally symmetric vibrations will appear. The Jahn-Teller effect is 
a special case of this in which certain intensity anomalies may occur. 

For a given upper state one of these processes may be. the 
principal factor determining the band shape} to determine it theoretically 
one would have to determine the stable geometric arrangement in the 
excited state and then apply the Franck-Condon Principle and vibrational 
electronic-coupling theory. We have carried out a much simpler and more 
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qualitative procedure. A crystal field perturbation of low symmetry 
is applied by using the point charge model and its effects on each of the 
various states are calculated. If a state is found to depend strongly 
on an unsymmetrical displacement, it ought to have a broad absorption band, 
and the normal vibrational modes which can produce the displacement should 
appear in the spectrum. The four components of the ^T^(F) state should 
have approximately the same behavior toward the totally symmetric 
distortion but they may behave differently toward other symmetry types. 

The components of this state therefore provide a good example on which 
to test this procedure. 

3 

b. Origin of Structure of the T^(F) Multiplet 

In Fig. 6 is shown the energy of the ^(P) multiplet, and the 
nearby levels which Interact strongly with it, as a function of Dq, the 
cubic crystalline field strength. If we set Dq a Dq Q bS(a lg ) this diagram 
shows how the levels change with a totally symmetric displacement S(A lg ). 

The slopes of all the lines are nearly the same, indicating that all the 
multiplet components will suffer the same totally symmetric displacement 
in their excited states. According to the bandwidth criterion of Orgel (A) 
the components would then all have the same bandwidth (making the reasonable 
assumption that the force constants are the same in each state). If the 
identifications given in Fig. 5 are accepted, this is clearly incorrect. 

The conclusion must be that the displacements occur in some non-totally 
symmetric coordinate, and that it is by a different amount, or a 
different coordinate in each state. 

Before going to this let us estimate the intensity distribution 
of the Franck-Condon series which would arise from a totally symmetric displacement. 
The magnitude of the displacement can be estimated as follows. The contraction 




6 . 

of the ion-ligand separation observed for NiO, compared to the hypothetical 

o < 3 > 3 

crystal not having a crystal field stabilization is 0.07 A . In the T^(F) 
state the energy decreases as Dq decreases, hence there should be an increase 
of the ion ligand separation. Relative to the ground state this expansion 
will be y| x .07 X, using the ratio of the dependences of the two effects 
on Dq. Taking 0.1 X = x as a reasonable estimate, the intensity ratio 
of the 0-0 and 0-1 bands according to the Franck-Condon Principle (for 
harmonic vibrations) is 


*10 . ggl 

I 00 


= 1 


„18 


where a - 4x cvm/h and is taken to be 2 x 10 , corresponding to v 5 425. Thus 
we should expect progressions of totally symmetric vibrations which have 
their maxima near the origin. This is not what is observed, and it reinforces 
the conclusion already arrived at that the non-totally symmetric modes 
are important in determining the spectral structure. 

The next possibility to be considered is the tetragonal distortion 
represented by a partner of the e^ normal mode of an octahedron. Figure 7 
shows this and other normal displacements. In Figure 8, the energy is 
plotted as a function of It is seen that each energy level behaves 

differently under this displacement. The calculation on which this figure 
is based is a complete diagonalizatlon of the d® energy matrix under a 
tetragonal perturbation*. The latter is defined by a point charge model 
with the potential: 

.*r‘Y;+Qr“v« + 4‘ [yJ+JII «*+,*) ] 

and y is the charge on the ions. 


r 2 v 0 ^ 

v tetrag' t v cubic " 2^ 

where A = rO^J “ -§-) J Q = V (-^ - ■ 


* The calculation was performed by machine under the direction of P. K. Baltzer. 
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; is the Ni-0 distance in the direction of the tetragonal distortioi 


(a) and a is the distance in the perpendicular direc 


The displacement 


S 2a = 2 <z 3 ' Z 6 ) " (X 1 ' V ' (y 2 ' y 5> 
where the atoms are numbered as in Fig. 7. For infinitesimal displacements 


A positive value of S^ a means a lengthening of the c- axl s. The calcula tion s 

were carried out for the ratio * ^’ 2 where P * an< * ^ ” 105 r • 

r 4 2 

This ratio is obtained when -- = .25 R , where R is the bond length, 

2 

This is close to the value given by SCF functions. 

This mode of distortion, e g is the one commonly considered in 
calculations of the Jahn-Teller effect^. It has the correct symmetry to 
split the T^ and! Tjg states of a cubic Configuration. Thte pdrtnhr 
is thwsht to lead to the lowest part of the resultant energy surface 
although a combination of S 2a and S 2b may be the lowest energy configuration 
under some conditions (5). In the present connection we assume that the 
upper states are all subjected to a tetragonal (S 2a ) potential which 
may or may not lead to appreciable distortions into a tetragonal shape. 

We will examine the band widths in the spectrum to see if they can be 
explained as vertical transitions to Buch distorted upper states. 

Figure 5 shows the pattern of bands to be compared to Fig. 8. 

The two lowest components and T^ are here identified with the two 
broad peaks at 14100 and 14400 cm 1 . Fig. 8 shows that the T^ band is 
highly unstable under a tetragonal distortion. For a negative distortion 
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of one percent of the bond length (z-bond length increased two percent, 
x-bond length decreased by one percent) the energy of the E component 
decreases 1000 cm The zero point excursions are about two percent 
of the bond length and the zero point energy in the e^ coordinate is 
about 150 cm -1 . Therefore the destabilizing forces are considerably 
greater than the restoring forces, and a distorted equilibrium con¬ 
figuration is to be expected. 

Both the T 2 and E components show minor instabilities which 
are no more than the zero point energy in a displacement equal to the 
zero point displacement. 

The Aj component can have no first order energy change since 
it is non-degenerate, but coupling to other Aj states leads to 
powerful second order effects of the distortion. Within half the zero 
point excursion, its energy is lowered by twice the zero point energy, 
and this state must become unstable in the cubic configuration. 

The sharp peaks therefore correspond to levels which are 
not unstable under tetragonal distortions and the broad ones to levels 
which are unstable. The good correspondence between theory and experi¬ 
ment for the ”*Tj(F) multiplet of NiSiF^CH^O is encouraging and suggests 
that such correlations will be found for other bands and for other 
substances. 

c. Vibrational Structure of 3 Tj(F) 

In this section it will be shown that two true electronic 
origins in the 3 T 1 <F) multiplet may be identified. These give the 
multiplet width, and provide an independent path to the conclusion 
reached in the last section. 
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The sharp bands at 15780 and 15970 cm" 1 (called bands III 
and IV) are the origins of progressions of about 500 cm -1 vibrations. 

Three members beyond the origin are discernable. The two peaks are barely 
distinguishable in the 0-2 region. In order to explain the shapes of 
the blurred 0-2 and 0-3 features, we must say that the 15970 progression 
falls off more rapidly than the 15780 progression. There is no observable 
shift of the interval from 15970 upon deuteration; a minimum shift of 
40 cm“‘ would have been detected. A 40 cm - * reduction on deuteration 
was observed for the interval measured from 15780 cm”*; but since the 
15780 peak is multiple, three parts being distinguishable, the shift 
may only represent differences in the relative intensities of the three 
0-1 members. It seems probable therefore that the 500 cm” 1 interval 
represents a vibration associated mainly with the Ni-0 bonds. The 
fact that it appears in a regular progression shows that it is totally 
symmetric in the symmetry group common to upper and lower states. The 
previous arguments suggest that this group is the symmetry group of the 
ground state, and that the 500 cm 1 interval is essentially of an 
octahedral molecule. 

One or two OH frequencies of the water molecules appear in the 
spectrum. Added to an 011 interval is a series of 500 cm” 1 intervals. 

Our object here is to assign these bands and determine their origin. 

There is a total of four normal modes of the NiOl^O)^ complex 
for each normal mode of the water molecule in the group and they are 
, Ag, Eg. Raman assignments have been reported for ZnSiFg'dH^O (6) 
where the vibrations should be very similar to those of the Ni ++ compound. 
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Although no deuterated material was studied, these vibrations 
hydrates are reduced by a factor of 1.36 + .02 (7). The fre.qi 
the vibrations given in Ref. 6, and the corresponding values 1 
deuterated molecule assuming this factor are: 


3456 

3492 



symmetric OH stretching 


>ther 
Les of 
:he 


3515 2582 Ag J 

C antisymmetric OH stretching 
3519 2585 Eg! 

Since we know the positions of the absorption bands we can subtract the 

Oil or OD frequencies to find their origins: 

H 19068 - f 3456 * f 15612 

{3513 {15553 

D 18176 - P 540 = P 5636 

(2582 {15594 

Since these do not coincide with observed bands, one of them must be a 
hidden electronic origin. The electronic origins of the ndrmal and deliberated 
Species should not coincide exactly because of the zero point energy shift. 
Since we h^ve been able to locate the origins of several electronic 
transitions of NiSiPg«6H 2 0 we have some idea of the magnitude of this 
shift. Going from H. to D there is a blue shift of 73 cm -1 for the 
^(F) band, 45 cm" 1 for the 3 Tj(F) band and 110 cm" 1 for ^(D). 

These numbers were determined from the fine structure near the origins of 
the bands. Thus it is reasonable to expect anything from about 30 to 
150 cm" 1 . The possible origins above are therefore reasonable. 

Whichever pair is chosen, it is necessary to assume that the 0-H 
vibration observed does not add to either of the observed bands 15793 
or 15957, but instead adds to some different origin. We suppose this to 
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be a true electronic origin. If so, the OH vibrations must be u, not 

g, since they have appeared in a single quantum in a g -► g electronic 

transition. The Raman data do not give us the u modes, but from the 
energy difference between the Ag and Eg-raodes we can see that the. 
coupling of the OH vibrations of different water molecules iB weak; 
therefore the u-modes will not differ much from the g. This leaves an 
uncertainty of perhaps 40 cm -1 in the calculated position of the electronic 
origin. 

An additional uncertainty arises from not knowing whether the 
symmetric or antisymmetric stretching mode is the more active one. In the 
D^O compound, the first 0-D peak is double with a splitting of about 
BO cm" 1 . This may indicate that both modes appear, with the lower one, 
supposedly symmetric, the more intense. Although there may be other 
explanations of this doubling we take this as a reasonable possibility. 

Then the position of the origin, (averaging over both the A and E components, 
since both should appear in the spectrum), should be 15594 (H) or 
15623(D). The deuteratlon shift is reasonable in magnitude, and in the 
correct direction, 

The electronic origin at 15594 could be associated with either 
of the sharp peaks 15793 or 15957, The respective intervals, which 
must be single quanta of a u-vibration, are 199 and 363 cm" 1 . Because 
of the uncertainty in the location of the origin, and our lack of knowledge 
of the infrated frequencies, we cannot choose between these two possibilities, 
since each is a reasonable one. In fact it is now possible to ascribe 
both peaks to the same electronic transition. 
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In an attempt to distinguish these, choices the temperature 
dependence of the intensity in this region has been measured. The 
intensity should increase with temperature according to the function 

— 0 /T 

f/f = 1 + e where & - 1.44 u, with to the vibration frequency 

in cm” 1 . The other parts of the ^Tj(F) band should increase more rapidly 
since many lower frequency vibrations are effective. Bands III and IV 
actually increase much less rapidly with temperature than the bands I and 
II at the low energy end of the 5 T^(F) absorption region. Figure 9 
shows the 3 T^(F) band at a series of temperatures. The extinction co¬ 
efficient at the 6260 X peak increases by 27 percent from 4.2°K to 
300°IC while in the 6900 X region it increases by 100 percent in the same 
temperature range. The formula given above actually predicts an increase 
of 17 percent in the total absorption for the 363 cm" 1 choice and 38 
percent for the 199 cm 1 choice. Since we do not know the true area 
attributable to the bands of interest we again cannot distinguish between 
the two possible u-vibratlon frequencies. However, the temperature 
dependence experiments do tell us that the interpretation of the Bharp 
peaks of the 3 Tj(F) band as 0-1 bands is essentially correct. 

The vibrational structure at the red end of the ^T^F) band 
is probably associated with the lowest of the predicted multiplet components 
Aj. This structure appears to begin at 12690 cm” 1 (12735 in the DjO 
compound). There is some evidence for a progression of 380 cm 1 
vibrational modes. As Fig. 5 shows, there is a series of sharp bands 
having the same intensity spaced 50 - 100 cm 1 apart, then, at 380 cm 1 
a rise of intensity followed by a similar series of small intervals. 

After four of these 380 cm -1 intervals the first peak ut. 14150 cm -1 is 


reached, 


This idealization of the band structure is illustrated in 
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Fig. 10. 

The approximately 380 cm 1 interval must be a non-totally 
symmetric g-vibrational mode because of its uniformly increasing intensity 
in the progression. Since the lower two levels of the 3 T^(F) multiplet 
were shown to be highly unstable with respect to Eg displacements, the 
380 mode is probably an Eg mode. Now, since the upper state cannot 
be reached except via u-vibrational modes, where, are these modes in the 
spectrum? They muBt surely be the series of low frequency modes observed 
at the beginning of the spectrum end, apparently, superposed on each 380 
quantum as suggested by Fig. 10. There is no a'priori reason why the 
approximately 363cm” 1 mode postulated for the 15594 band should 
not also appear, and perhaps it is present but not resolved. 

The new and interesting feature of this band, according to 
the foregoing interpretation, is that the odd modes do not induce a 
large transition moment bv themselves, but do so only when accompanied by 
several quanta of the even Eg mode. This is not surprising, however. 

The simultaneous excitation of the two kinds of vibrations is necessary 
In order both to destroy the center of symmetry and to provide a large 
vibrational overlap integral. 

It is not known what the low frequency modes are. The 230 
could be one of the modes of the Ni(H 2 0)6 octahedron, but 35 and 88 
are too small for internal modes. They may correspond to motions mainly 
of the entire complex having some contributions from internal motions. 
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It has become obvious that a complete study of the infrared 
sppctrum of the crystal is needed in order to analyze the optical spectrum. 
But for the main purpose of this paper, it is enough to have established 
that there is an electronic origin of some multiplet component at 
15594, or possibly about 50 cm -1 lower; and another origin at about 
12fi90 cm -1 . 

The difference between the two established electronic origins 
of the 3 Tj(F) band is now known to be about 2850 - 2900 cm”*. This is 
probably the entire electronic width of the multiplet. It is more than 
twice as wide as one calculates from the "complete" theory of Ni ++ in 
a cubic field using X B - 275 cm” 1 . This value of X is consistent 
with measured g-factors and it ought to be about the same for excited 
states as for the ground state. The discrepancy therefore cannot be 
blamed on the value of X. The extra multiplet width undoubtedly a 
result of the instability of the lower multiplet component which causes 
its electronic origin to be red-shifted by the mechanism discussed in 
Sec. 11, and illustrated by Fig. 8. The additional multiplet width 
therefore represents the stabilizing energy gained by the lowest 
multiplet component by the tetragonal distortion. 

IV. CONCLUSIONS 

We have now shown by two experimental criteria, measurement of 
band widths, and measurement of electronic energy differences in a 
multiplet, that components of the 3 T^(F) multiplet have different 
equilibrium shapes. The other bands of Ni ++ have smaller spin-orbit 
coupling, so that the individual components are not resolved. But as far 
as can be learned, these bands also contain upper states which are unstable 


in the cubic configuration, 




I 


1 

I 
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Table I 

Structural Data for NiSiF »6H 0 
6 2 

a) one formula unit per unit rhombohedral cell, a = 6.21 a = 96°53' 

b) approximate atomic positions in rhombohedral coordinates 

Si 0 0 0 

Ml 1 I I 

2 2 2 

F + (xyz, zxy, yzx); x = 0.14 y B 0.31 z c 0.94 
0 * ( " " " ) x s 0.64 y ■ 0.01 z « 0.44 

c) nearest Ni-0 distance; 2.07 Z 

d) angle between C-axis and Ni-0 direction; 55.8° 

(for regular octahedron, 54.8°) 


i 

I 

i 

1 





FIGURE CAPTIONS 


1. Absorption spectrum of nickel fluesilicate. The portion from 8000 X 
to shorter wavelengths is taken from a crystal of NiSiFg^H^O at 77 K. 

The portion from 8000 $ to longer wavelengths is taken from a 
crystal of NiSiFg*6D o 0 at about 10°K. 

The theoretical positions of the absorption bands calculated using the 
following parameters are also shown: Dq = 910, B = 880, A = 300 
C/B = 4.5. 

2. Fine structure near origin of 3 T^ (F) band. 

3. Fine structure near origin of ^(D) band. 

4. Fine structure near origin of 3 T^ (P) band. 

5. The 3 T, (F) band of NiSiF *6D 0 for two thicknesses of crystal at 4.2°K. 
The apparent origins III and IV are marked as well as the progressions 
following them. 

6. Energy of multiplet components of 3 T. (F) as a function of a totally 
symmetric distortion (Dq). 

7. Normal displacements (or vibrations) of an octahedron showing standard 
notation for frequencies. 

8. Energy of multiplet components of 3 T, (F) as a function of tetragonal 

distortion, eg (the partner shown in rig, 7). The diagram has been 

calculated for the parameter values Dq = 900, B = 800, C/8 ■ 4.5, 

A * 320. The tetragonal distortion potential is based on a point 
charge model (see text). The horizontal scale width is about equal 
to the zero point displacement in the e. vibration. The negative 
qign corresponds to an elongation of th§ octahedron. 

9. The absorption in the 3 Tj (F) region at various temperatures (°K). 
Successive tracings are separated in the vertical scale. 

10. /in idealization of the structure near the origin of the 3 T, (F) 
band showing the progression of 380 cin which appears to Be 
present. This figure is to be compared to the original record 
shown in Fig. 2. 
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Sec. XV 

Exchange Coupled Mn Ions in ZnS Single Crystals 
Donald S. McClure 
RCA Laboratories, Princeton, N.J. 

ABSTRACT 

Evidence for exchange coupling between Mn ions in ZnS : MnS mixed 
crystals is manifested in the optical spectrum in two ways; the 
intensity of absorption is not proportional to concentration of Mn ++ 
and temperature dependent absorption lines appear. The latter 
give a measure of the exchange integral in the ground state of 
25°K, which is larger than the exchange integral of 12°K derived 
for pure MnS from its high temperature susceptibility curve. There is 
also spectral evidence for triples of Mn* + . These appear at concentrations 
above 8 percent and appear to be responsible for luminescence 
quenching. 





Exchange Coupled Mn ++ Ions in ZnS Single Crystals 
Donald S. McClure 
RCA Laboratories, Princeton, N. J. 

There is now both optical and microwave evidence for the 
existence of exchange-coupled ion pairs in crystal lattices. Owen (1) 
haB studied the Mn + pairs in MgO by spin-resonance, finding a value 
of the exchange integral J, of 30 cm -1 for the 110 pairs. This is 
about twice the value expected by comparison with the «J value in 
MnO where the lattice constant is only slightly smaller than in MgO. 

Statz and Koster (2) and Schawlow (3) have observed Cr +3 pairs in 
Al^Oj. The present note concerns optical evidence for Mn pairs 
in ZnS. The initial indication of interaction between Mn ++ ions 
came from a study of the 5300 X broad band absorption of the Mn ++ 
ion at various concentrations in ZnS single crystals. The intensity of 
absorption increased faster than the concentration. Later, experiments 
at 4,2°K and above on the absorption spectrum in the sharp line region at 
4600 X, revealed temperature dependent lines which were absent in dilute 
samples. These extra lines arise from the thermally excited pair 
states. 

Experimental! The specimens of ZnSsMnS were prepared by the vapor 

* 

transport method using iodine as the carrier. The deposition temperature 
varied from 750 to 9Q0°C, and the Mn ++ content varied from 1 to 11 
percent. Most of the samples were the cubic form of ZnS* They were 
from 1 to 5 mm in largest dimension. Pieces of crystal were polished 
with parallel faces and mounted in our optical dewar in the vacuum space. 
Thermal contact was provided by silver paste. A copper-constantan 
thermocouple was silver pasted to the sample. 

*Prepared by R. Nietsche and M. Lichtensteiger, RCA Zurich. 
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Spectra were obtained with a Cary 14 spectrophotometer. Temperatures 
between 4.2 and 77°K were obtained by letting the liquid helium evaporate 
and making the spectral runs during warm-up. The temperature change 
during such a run was never more than 5°, 

Results: Figure 1 shows spectra of crystals of several concentrations, 

including some pure MnS. The latter has the NaCl structure rather than 
the ZnS structure, so is not strictly comparable with the mixed crystals. 

The crystal thickness and approximate percentage of Mn are given in the 
figure and are also shown in Table I. The percentage Mn was determined 
by chemical analysis of crystals which deposited during preparation 
near the one chosen for study. There was a nearly uniform decrease of 
Mn concentration from the hottest region of the deposition zone to the 
coldest. Analysis of the actual samples must await completion of spin- 
resonance studies. The product of concentration and thickness is pro¬ 
portional to the number of Mn atoms in the light path, and in the 
absence of interaction between Mn ++ ions this product would be proportional 
to the optical density. Consequently we give in Table I the molar 
extinction coefficient,£ , at the peaks of the main absorption bands, where 

6 = P/Cl 

p » optical density, log I q /I, C * concentration in moles/liter 
1 » path length in crystal in cm. For a given band, this quantity should 
br a constant independent of concentration, but as Table 1 shows, it 
increases with concentration. In view of the uncertainty of the chemical 
analysis, and the obvious exception to the general trend of results in 
the case of sample 263B, the conclusions to be drawn from Table I are 
somewhat weakened. It is safe to say that the absorption coefficient 
increases with concentration, but we cannot say by exactly how much. 
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Figure 2 shows the results of an experiment in which the spectrum of 
a 5 percent crystal was obtained at various temperatures. Only the 
4600 8 region showing sharp line structure was covered by these 
experiments as the other regions did not show any dramatic temperature 
effects. The arrows in Fig. 2 point out the positions of absorption 
lines which are absent at 4.2°K but which grow with rising temperature. 

These evidently must arise from a level (or levels) near the ground state 
and within kT of it. 

The optical density of a temperature dependent line was corrected by 
subtracting the optical density of the background at the same wavelength. 

The background was taken to be the optical density at 4.2°K where there 
was no trace of the temperature dependent absorption line. This procedure 
would give the line strength if the width of all absorption lines in this 
region were to remain constant with rising temperature. Some broadening 
does occur but it is very slight. The rate of change with temperature 
is probably underestimated slightly. Table II gives the optical density 
as a function of temperature and the corrected optical density for a 
sample containing 5 percent Mn. If the logarithm of corrected optical 
density is plotted against 1 /T, a straight line is obtained whose slope 
is given by a Boltzmann factor. The 4653.4 A peak gave 29.6°K and the 
4590 £ peak gave 35.8°K, Considering the errors involved in the corrections, 
especially for the weaker 4590 £ peak, these two results are probably 
within experimental error. It will be seen that the two lines almost certainly 
arise from the same thermally excited level near the ground state. 

A second sample of the 5 percent crystal gave © = 28.5°K for the 
4653.4 £ peak. 
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A sample of 2 percent Mn in ZnS (1.27 mm thick) gave a value 
©= 37.8°K for the 4653.4 X peak (see fig. 3) . Although this is again 
somewhat high it is felt that the same type of thermally excited level 
as is present in the 5 percent sample is involved. 

Figure 4 shows the temperature dependence of the 4600 X region in 
an 8 percent sample. The structure is not as clearly marked, but the 
intensity increases strongly with temperature. 

The correct analysis of the data is actually somewhat more involved 
than the use of a simple Boltzmann factor. This will be discussed in 
the next sectibn. 

Discussion: The absorption spectrum of Mn ++ arises from transitions 

between a 6 S ground state and various quartet excited states. For 

purposes of orientation, the spectrum of the trans- Mn ( H 2 °^4 C1 2 molecul ® 

is shown in Fig. 5, All of the spectral assignments are shown in this 

figure, using the approximation that the crystal field is cubic. The 

resemblance of this spectrum to that of Mn ++ in ZnS is obvious. 

Assignment of transitions in the spectrum is based on the similarity of all 

Mn ++ spectra and on the theoretical term splitting diagram shown in 

Fig. 6. The parameter Dq measures the strength of the cubic field, and 

in order to fit the spectrum of Mn in ZnS, we must take Dq ■ 260 cm - *. 

The parameter, B, which measures the atomic term separations does not 

give all the level positions very well, but is about the same as for the 

Mn + * halides, 615 cm *• A better enalysis of the term positions of the 

ion is provided by the covalency parameter of Stout (4). The value of 

this parameter must be quite large for MnS in ZnS, and larger still for 

* The "2 percent” sample, ft 263B, is probably closer to 4 percent, 
according to Table I. 









MnSe in ZnSe. Fig. 7 shows the covalency parameter for a number of Mn + 
compounds. The single parameter does not give a perfect fit by any 
means, but it does arrange the ligands in an order which makes sense 
chemically. 

The transition in the region of 4600 X in MnSiZnS has two electronic 

l ^ 

levels for its upper states, the A^ and E levels, which are degenerate 

in a purely ionic model of the crystal field. The effect of covalency 

4 

on the two is somewhat different, however, and the E level is expected 

4 

to lie somewhat below . This is indicated in Fig. 7. It is 
difficult to separate these two levels in the spectrum, sb their 
absorption regions overlap. In the MnS;ZnS spectrum it appears that the 
origins of these transitions are: 4 E 4700 X <21270 cm -1 ); 4 A^ 4637 X 

(21560 cm" 1 ), a separation of 290 cm" 1 . The additional lines of about 
100 cm -1 separation from the origin are probably vibrational intervals. 

When the MnS concentration is above one percent, there is an appreciable 
statistical probability that for a given Mn ++ there will be a like ion 
in the neighboring cation lattice Bite. These coincidences will be referred 
to as pairs. There is only one geometrically distinct type of pair in the 
cubic ZnS lattice. At concentrations above 5 percent the number of 
triples becomes appreciable. There are many ways to place three foreign 
cations in the ZnS crystal such that only a single anion separates two 
of them. Fig. 8 shows the totals in these categories. The sum in all 
categories should be unity; at a composition of 10 percent only 89 percent 
of the atoms are accounted for by single, pairs and triples, so that 
quadruples and higher clusters are beginning to become important. 
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The exchange coupling in antiferromagnetic materials is usually found 
to be given to a good approximation by the isotropic spin-spin coupling 
hamiltonian: 


H = - 2J S 1 *S 2 (1) 

where J is the exchange integral, and are the total spins of the 

two ions. This same approximation should be equally valid for the 
ground state of ion pairs in ZnS. The resultant energy level system is 


given by: 


B * J <-X (X + 1) + 2S(S + 1 )) 

where X * Sj ♦ S 2 Sj + S 2 - 1 -0 t and S = j s i | “ 1 S 2 I • 

For J < 0, X ■ 0 is the lowest level, a non-degenerate level. 

The next higher level is X = 1 at 2J, a triply degenerate level. The 
partition function for the system is? 

X(X + 1) 0 /T 


X»0 


<2X + 1) e" 

where © “” J A. We need to know the population of the level X * 1 as a 
function of temperature. This is: 

»,«> - 20/T 

E T 


(4) 


N o (°) 

The intensity of the temperature dependent lines should be proportional 
to Nj(T), This must replace the Boltzmann factor used earlier. Q may be 
obtained from tables of the partition function of a rotator. 

For the 2 percent sample, the value of J/k obtained from this analysis 
is - 25 ♦ 2°K. Extremely accurate analyses of the high temperature 
susceptibility curve of p-cubic MnS (ZnS structure) by Doneban and Stevens (5) 
give J/k = - 12,45°K. One would not expect such a large difference between 


these two values because the metal-sulfur distance for p-MnS is 2,425 X, 
and that for ZnS is 2.35 X. However, OwenB results on Mn in MgO (1) were 
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similar to these. 

The failure of Beer’s law cannot be explained by the use of the 
hamiltonian (1). Physically, the phenomenon can be understood in the 
following way. The transition in single ions is forbidden by the selection 
rule AS = 0, since the ground state is a sextet, S = 5/2, and the excited 
states are quartets, S = 3/2. The transitions in single ions are 
observed only because spin-orbit coupling partially destroys the spin 
as a good quantum number. In the pairs, there are many values of the total 
spin, as we have seen in dealing with the ground state. The total spin 
in the ground state ranges from 0 to 5 while in the excited state, 
assuming only one excitation per pair we have S = 5/2 + 3/2, 5/2 + 3/2 - 1, 

- 5/2 - 3/2, or 4, 3, 2, 1. As we have just seen, the ground state 

having S = 1 is populated appreciably at' 20°K and above, and it may make 
transitions to an upper pair state having S = 1 without violating the 
spin selection rule. Even so, the transitions would not be any more 
intense than single ion transitions if the operator (1) were the only 
effect as an analysis of the eigenfunctions of the operator Cl) will 
show. However, the same spin states are generated by any other inter¬ 
action operator, since they arise because of symmetry properties, and it 
is only a matter of finding the operator which has the effect of increasing 

the intensity of the S - 1 -► S = 1 transitions above that due to the 

single ion transitions. 


One form of such an operator is (6) 

- - ^ 

where i,k sum over all spins of the two ions and 1,2 refer to the two ions. 


This operator is also isotropic and provides for the possibility of uncoupling 




the spins from the total spin S within single ions. The total spin 
of a single ion is therefore not preserved as it was by operator (1), 
and this provides the possibility of intensity increments in the transition. 
We do not expect that (5) will have much effect on the energy so that (1) 
is still valid in connection with the discussion of the ground state 
splitting. 

A theory of the intensity has not been worked out quantitatively, nor 
has the theory of the splitting of the excited states. The two lines at 
4590 and 4653 are probably the S *? 1 ltnea of excited states derived from 
4 6 4 6 

a Aj ♦ S and a E ♦ S pair respectively. 

According to the spin selection rules, we should not see absorption 
from the ground state of a pair to its excited states, since there are 
no excited pair states having S B 0. However, the absorption spectrum 
should actually be observable through the same mechanism which permits 
transitions in isolated ions. In this case then should we not see these 
transitions appearing like single ion transitions but blue shifted by 35/2 J 
or about 440 cm -1 ? They do not obviously appear in the spectrum. They 
would be weak because there is no mechanism to intensify them, but no 
weaker than absorption due to isolated ions. On a statistical basis 
in a 5 percent sample there are 54 percent of Mn ions having no near 
neighbor Mn , and 34 percent having orte ffi%. Mn + *. Thus the spectra 
due to unexcited pairs ought to be nearly as strong as the spectra of 
individuals. 

It would be desireable to unravel the spectra due to the two thermally 
accessible pair states (at low temperatures) from the spectra of individuals, 
but so far this has not been possible. There are definite differences. 
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however, between the spectra of 2 percent and 5 percent samples in the 
temperature independent regions as Figures 2 and 3 show. The spectrum of 
the one percent sample is so weak that only the two prominent lines at 
4653 and 4590 % are observable. 

The spectrum of the 8 percent sample has temperature dependent regions, 
(Fig. 4) but individual lines are not well resolved. It appears that 
at the lowest temperatures the absorption is stronger than in the 5 
percent sample. There are twice as many triples in the 8 percent sample, 
and they may make an important contribution. When three Mn ++ ions 
couple in three equivalent pairs, a X = 3/2 state is lowest (J < 0) and 
since X*3/2 Is possible in the excited states, thermal excitation is not 
required to match the multiplicities of the combining levels. This is 
probably why the extinction coefficient becomes rapidly larger at 8-11 
percent where the number of triples is comparable to the number of pairs 
and singles. There are several different kinds of triples, whereas in this 
lattice there is only one kind of pair; therefore the spectra of the 
more concentrated crystals become very complicated. In the 11 percent 
samples, no detailed structure may be observed for this reason. 

There la an Interesting connection between this work and luminescence 
efficiency. It la known that Mn phosphors begin to lose their efficiency 
when the Mn concentration is higher than one or two percent. Leverenz 
and North(7) showed that above one percent, the number of isolated 
ions in Zn 2 SiO^ (rbhdl) begins to decrease rapidly, following the curve 
of efficiency rather closely. Qualitative observations of the ZnS:MnS 
crystals showed that the 1 and 2 percent samples fluoresced brightly 
under 3660 X illumination, the 8 percent sample was appreciably dimmer and 
the 11 percent sample had no observable fluorescence. These results are 
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similar to those reported for ZtijSiO^rMn. The quenching of fluorescence 
in the case of ZnS:Mn appears, however, to be more closely correlated 
to the increasing number of triples rather than pairs, as can be seen by 
referring to Fig. 6. 
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Table I 


Molar extinction coefficient, £ , at 4.2°K for principal peaks in ZnS:MnS 
absorption spectrum as a function of concentration. C, in mole percent, 
is the analysis of a batch of crystals rather than the one for which the 
spectrum was measured, and is therefore only approximate. The so-called 
2 percent sample must be closer to 4 percent. For the 100 percent 
sample, pure MnS, the £ values are given for the bands which correspond to 
those of ZnSjMnS mixtures. Peak £ values for a b*nd are given in all cases. 





5300 8 


5100 8 


4600 8 

4200 8 

Sample C,°/ 0 

l,mm 

P 

£ 

P 

£ 

P 

L 

P 

L 

262A 

1 

1.05 

0.02 

0.41 

0.085 

1.80 

0.115 

2.40 



263B 

4 

1.27 

0.13 

1.30 

0.43 

3.65 

0.60 

5.10 

0.52 

4.40 

263C-2 

5 

1.18 

0.18 

0.66 

0.55 

2.04 

0.79 

2.89 

0.65 

2.38 

262D 

8 

0.55 

0.16 

0.78 

0.55 

2.72 

0.70 

3.40 

0.54 

2.63 

224 

11 

0.48 

0.67 

2.72 

2.03 

0.20 

2.52 

10.40 

1.46 

5.95 

246 

1Q0 

0.05 

0.66 

2.80 

0.61 

2.55 

1.62 

6.80 

- 

- 




Table II 


Strength of temperature dependent line at 4654 X as a function of temperature 
in two samples. p»optical density. The corrected optical density, e.g., 
p -.35 is divided by the maximum optical density in the last column for 
each sample. 


5°/o Sample, # 263C-1 
4654 X 


o ii 

"2 /o Bample U 265B 
4655 X 


T°K P 

4.2 0.35 

18 .52 

28.5 .65 

33.8 .71 

43.6 .78 

522 .84 

64.5 .86 

77.5 .91 

88.5 .91 



T i /{>• 23 \ 

\ .* 5 / 

4.2 .0 

12.5 .01 

17.5 .211 

22.5 .373 

26.0 .477 

29,0 .544 

34.0 .677 

39.0 .755 

46.0 .867 

58.0! 1.01 

68.0 1.00 



FIGURE CAPTIONS 
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1. a. Absorption spectra of mixed crystals of ZnS:MnS at 4.2°K for 
various concentrations. All are cubic zinc blende structures except 
for the pure MnS which has the NaCl structure. The concentrations, 
thicknesses and sample numbers are given in the figure. 

b. Spectrum of ZnS:MnS (approx 5 percent MnS) at 4.2 and 77°K, showing 
spectral details. Note that 21481 and 21780 appear only in the 77°K 
spectrum. The line 21262 at 4.2° and 21280 at 77° are probably 
the same transition. 

2. The 4650 X band of ZnS: 5 percent MnS at a series of temperatures showing 
the two temperature sensitive bands (marked by arrows). 

3. The 4650 X band of ZnS: 2 percent MnS at a aeries of temperatures. 

4. The 4650 X band of ZnS: 8 percent MnS at a scries of temperatures. 

The one distinct temperature sensitive band is marked by an arrow. 

5. The absorption spectrum of Mn (HjO^Cl, single crystal at room temperature 
showing the spectral assignments. 

6. The energy of the Mn* + states in a cubic field as a function of Dq. 

7. The covalency parameters for various compounds of Mn . 

B. Statistical probability of various species of aglomerates 
in ZnSjMnS mixed crystals. 
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